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Resumen 


Transmisión y conductancia en 


dispositivos electrónicos basados en 
grafeno 


El objetivo principal de esta tesis ha sido caracterizar analíticamente el 
transporte electrónico a través de una lámina de grafeno sometida a distin- 
tas puertas de potencial y/o constricciones. Hemos estudiado el coeficiente 
de transmisión en una lámina de grafeno sometida a diferentes barreras de 
potencial unidimensional. El hecho de que el coeficiente de transmisión sea 
unitario para ciertos ángulos, provoca que no pueda existir confinamiento de 
portadores de carga mediante una barrera de potencial eléctrico. Sin embargo, 
es posible obtener un pseudogap de energía en el coeficiente de transmisión, 
para ciertos valores de los parámetros, si se dirige de un modo adecuado el haz 
de incidencia de los electrones en la lámina. La introducción de una segunda 
barrera de potencial origina una serie de picos de resonancia en el pseudogap 
cuya posición ha sido identificada analíticamente. Finalmente, se ha conside- 
rado el método de la matriz de transferencia para construir el coeficiente de 
transmisión de una heteroestructura formada con un número indeterminado 
de barreras. 

A continuación, consideramos la aplicación de los resultados anteriores 
a un dispositivo físico real: un transistor de efecto de campo de grafeno 
(GFET). Construimos un modelo para el GFET y determinamos su con- 
ductancia dependiendo de la dirección de impacto de los portadores y su 
ponderación sobre la barrera de potencial. Describimos así una distribución 
de portadores incidentes alejada de la configuración uniforme característica 
de muestras homogéneas, y cuantificamos de este modo los efectos de inho- 
mogeneidad en la lámina de grafeno. 

Con el fin de mejorar y refinar nuestro modelo, hemos introducido el cam- 
po eléctrico asociado a la diferencia de potencial entre surtidor y drenador 
Vos, por medio de un potencial efectivo lineal en ciertas regiones de la lámi- 
na de grafeno bajo las top gate, que introducen una barrera de potencial 
proporcional a la diferencia de potencial entre la puerta superior e inferior 
Vra — Vga. 


En este proceso de someter a los portadores de corriente a diferentes im- 
pedimentos a su paso que controlen la corriente, sustituimos las láminas de 
anchura infinita por nanocintas de anchura fija, y las puertas de potencial 
por estrechamientos en la cinta. Con este objetivo, en el último capítulo es- 
tudiamos los efectos de una constricción rectangular en la conductancia de 
una nanocinta de grafeno. Obtenemos una expresión analítica de dicha con- 
ductancia, que pone de manifiesto cómo puede manipularse la relación de 
anchuras en las diferentes regiones de la cinta, para modular la corriente, y 
evidencia comportamientos muy diferentes dependiendo del carácter metáli- 
co o semimetálico de la zona del estrechamiento y del resto de la cinta. En 
general, veremos que la conductancia va a venir dada por la conducción de 
la parte más estrecha de la cinta. Finalmente, hemos comparado nuestros 
resultados para distintas cintas con parámetros adecuados, que exploran la 
diversidad de casos posibles, con los obtenidos por otros grupos que tratan 
el mismo problema con métodos de tight-binding convencionales. Aunque el 
acuerdo no siempre es del todo satisfactorio, creemos que este trabajo analíti- 
co, unido a una mayor potencia de cálculo, puede ser útil en la descripción 
del comportamiento de estos dispositivos y algunos otros relacionados. 


Palabras clave : Grafeno, GFET, nanocinta, nanoconstricción. 


Abstract 


Transmission and conductance in 
electronic devices based on graphene 


The main objective of this thesis has been to characterize analitically the 
electronic transport through a graphene sheet subjected to different potential 
gates and/or constrictions. We have studied the coefficient of transmission in 
a graphene sheet subjected to different barriers of one-dimensional potential. 
The fact that the transmission coefficient is unitary for certain angles, cau- 
ses that there can be no confinement of charge carriers through an electric 
potential barrier. However, it is possible to obtain a pseudogap of energy in 
the transmission coefficient, for certain values of the parameters, if directed 
in an appropiate way the incident beam of incidence of electrons in the sheet. 
The introduction of a second potential barrier originates a series of resonance 
peaks in the pseudogap whose position has been analytically identified. Fi- 
nally, it has been considered the method of the transfer matrix to build the 
transmission coefficient of a heterostructure formed with an indeterminated 
number of barriers. 

Next, we consider the application of the previous results to a real phy- 
sical device : a graphene field effect transistor (GFET). We build a model 
for the GFET and we determine its conductance depending on the direction 
of impact of the carriers and their weighting on the potential barrier. Thus, 
we describe a distribution of incident carriers away from the uniform charac- 
terization of homogeneous samples, and quantify in this way the effects of 
inhomogeneity on the sheet graphene. 

In order to improve and refine our model, we have introduced the electric 
field associated with the potential difference between source and drain Vps, 
by means of a linear effective potential in certain regions of the sheet of grap- 
hene under the top gate, which introduce a barrier of potential proportional 
to the difference between top gate and back gate potential Vra — Vae. 

In this process of subjecting the current carriers to different impediments 


in their path that control the current, we replace the sheets with infinite 
width by fixed width nanoribbons, and the potential gates with narrowings 
in the ribbon. With this objective, in the last chapter we studied the effects of 
a rectangular constriction on the conductance of a graphene nanoribbon. We 
obtain an analytical expression of this conductance, which reveals how the 
relation of widths in the different regions of the ribbon can be manipulated, to 
modulate the current, and evidence very different behaviors depending on the 
metallic or semi-metallic character of the narrowing area and the other areas. 
In general, we will see that the conductance will come from the conduction 
of the narrowest part of the ribbon. Finally, we have compared our results 
for different ribbons with appropiate parameters, that explore the diversity 
of possible cases, with those obtained by other groups that treat the same 
problem with conventional tight-binding methods. Although the agreement 
is not always completely satisfactory, we believe that this analytical work, 
together with a greater calculation power, can be useful in describing the 
behavior of these devices and some other related ones. 


Key words : Graphene, GFET, nanoribbon, nanoconstriction. 
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Capítulo 1 


Introducción 


1.1. Prólogo 


Esta tesis pretende combinar la parte teórica y la parte experimental en el 
conocimiento y comprensión del transporte electrónico en grafeno. Relaciona 
los campos de la Macánica Cuántica, la Electrónica y la Nanotecnología, en 
consonancia con el tiempo actual en el cual se demanda un conocimiento 
interdisciplinar. A lo largo de ella consideraremos el transporte cuántico de 
los portadores de carga en un material llamado grafeno. El término graphe- 
ne fue introducido por H. P. Boehm como una combinación de la palabra 
graphite refiriéndose al carbono en su forma cristalina ordenada, y el sufijo 
-ene asociado a los hidorcarburos aromáticos policíclicos en los cuales los 
átomos de carbono forman estructuras en anillo hexagonal [1]. Existen dos 
tipos extremos de transporte dependiendo de la extensión del medio y de la 
longitud del camino libre medio de los portadores. Cuando las dimensiones 
del medio (en nuestro caso las muestras de grafeno) a través de las cuales se 
propagan los portadores son superiores a su longitud del camino libre medio 
estamos en un régimen difusivo [2], mientras que si sucede lo contrario, el 
régimen de conducción es balístico [3]. La longitud del camino libre medio de 
los portadores depende mucho de la proporción de impurezas en la muestra, 
que actúan como centros de dispersión, y de la temperatura, que aumenta la 
probabilidad de colisión. En el grafeno de alta pureza, los electrones pueden 
viajar largas distancias del orden de micrómetros antes de colisionar, inclu- 
so a temperatura ambiente [4]. A bajas temperaturas (< 40 K) se llegan a 
alcanzar longitudes del camino libre medio tan elevadas como 28 ¡um [5, 6). 

Por todo ello, podemos afirmar que el transporte electrónico en grafeno, 
para muestras de dimensiones nanométricas y a bajas temperaturas posible- 
mente transcurre en régimen balístico. 
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En el título de esta tesis, con transmisión nos referimos al coeficiente 
de transmisión que calcularemos más adelante en los siguientes capítulos 
bajo diversas circunstancias. Sin embargo, el coeficiente de transmisión no 
es una cantidad medible macroscópicamente. No obstante, está relacionada 
(mediante el formalismo de Landauer) con la conductancia, para la cual sí se 
pueden realizar medidas experimentales. Recordamos que la conductancia 
relaciona de manera macroscópica la corriente total / con la caída de voltaje 
V, mediante l =G-V. 


La elección del grafeno viene motivada por ser un nanomaterial con una 
enorme potencialidad al permitir el transporte balístico a temperatura am- 
biente. Sin embargo, a raíz de su descubrimiento, se comprobó que ésta y 
otras propiedades inusuales no son exclusivas del grafeno. Éste es sólo el pri- 
mero de un nuevo grupo de materiales conocidos como materiales 2D a los 
que, por semejanza, se podría aplicar también nuestros métodos y resultados 
con sólo ligeras modificaciones. Todos ellos se aventuran, por sus propiedades 
inusuales, como el futuro en la Electrónica. 


1.2. Historia del grafeno 





Figura 1.1: Grafito. 


El carbono es la base de la vida, ya que forma parte de todos los compues- 
tos orgánicos. Aparece en la corteza terrestre en un porcentaje no superior al 
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0.1% . Se presenta comúnmente en su forma tridimensional en dos variedades 
alotrópicas, el grafito, representado en la figura Fig. 1.1 y el diamante. En 
las últimas décadas han ido apareciendo nuevas formas cristalinas como los 
nanotubos de carbono, los fullerenos y el grafeno. Aunque tienen en común 
con la mayor parte de los anteriores, la colocación de átomos de carbono 
en una disposición de red hexagonal (el diamante tiene una red cúbica), sin 
embargo presentan dimensiones muchísimo más reducidas, por lo que se se 
les denomina sólidos de baja dimensión o nanomateriales. Previo a la irrup- 
ción en el panorama científico del grafeno, los nanotubos de carbono fueron 
descubiertos en 1991 por Sumio lijima [7], el cual recibió el Premio Príncipe 
de Asturias de Investigación Ciéntifica y Técnica en 2008. Este estudio de los 
nanotubos de carbono sirvió para estudiar luego más fácilmente el grafeno, al 
tratarse también de un material de bajas dimensiones compuesto por átomos 
de carbono. 





Figura 1.2: Imágen de grafeno monocapa (zona con menor contraste) obte- 
nida en el laboratorio mediante un microscopio electrónico. 


El grafeno es un material obtenido experimentalmente por primera vez 
en un laboratorio en el año 2004 por André Geim y Konstantin Novoselov 
pertenecientes a la Universidad de Manchester [8] determinando sus asom- 
brosas propiedades (aunque existían trabajos experimentales previos como 
puntualiza Geim en [9]). En esa época se creía que debido a cálculos ter- 
modinámicos previos, no podían existir cristales aislados en dos dimensiones 
(2D) de manera estable. Frente a dicho prejuicio, los investigadores de Man- 
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chester realizaron experimentalmente con éxito el aislamiento de una lámina 
de grafeno, para ello usaron el método de "escisión micromecánica?. Éste es 
un método rudimentario pero eficaz, con el cual se consigue aislar una única 
capa de grafeno partiendo de grafito, que está constituido por un conjunto 
de ellas. El método consiste en utilizar cinta adhesiva para ir exfoliando el 
grafito sucesiva y repetidamente, pegando y separando, pegando y separan- 
do, hasta que únicamente quede una monocapa. Posteriormente es necesario 
depositarlo sobre un substrato, por ejemplo en una oblea de silicio. En el la- 
boratorio obtuve grafeno monocapa siguiendo el proceso anterior, tal y como 
se observa en la imagen obtenida mediante un microscopio electrónico en la 
figura Fig. 1.2, donde la tonalidad más clara corresponde a una única capa, 
y por tanto grafeno. La obtención experimental por primera vez del grafeno 
dió lugar a que sus autores recibiesen el premio Nobel de Física en el año 
2010 [10]. Cabe reseñar que aunque la realización experimental del grafeno 
ha tenido lugar en el siglo XXI; mucho antes, en 1947, ya se había tratado 
de forma teórica el grafeno [11] como derivado del grafito. Otro ejemplo de 
que el estudio teórico precedió a la obtención experimental del grafeno, es 
el estudio del comportamiento de una partícula relativista con masa no nula 
descrita por medio de la ecuación de Dirac, y la transmisión de esta partícula 
a través de una barrera de potencial por efecto túnel [12]. 


Existen muchas características sorprendentes que posee el grafeno frente 
a otros materiales, sin embargo nuestro interés se centrará en las propiedades 
electrónicas. Entre las diversas características que hacen que el grafeno sea 
un material muy diferente a los semiconductores convencionales están: 


= La estructura de bandas del grafeno ideal (sin impurezas) no presenta 
gap de energía entre la banda de valencia y la de conducción [13,14], 
es decir es un semimetal. 


= La relación de dispersión energía - momento es lineal en las bandas de 
conducción y valencia [15], frente a la relación cuadrática obtenida para 
los semiconductores convencionales. 


= Los electrones en el grafeno se comportan como partículas relativistas 
sin masa (ultrarrelativistas) [15], lo que confiere al grafeno una alta mo- 
vilidad electrónica. Experimentalmente se encuentra que la movilidad 
1. puede exceder 15.000 cm? V=! s7! a una concentración de portadores 
de 101% cm”? [16]. 


La movilidad extraordinariamente elevada a alta temperatura de los por- 
tadores de carga, tanto como 2-10? cm? V7! s7! [17], indica que los electrones 
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en el grafeno viajan en régimen balístico. Esta movilidad tan elevada [4] es 
una característica deseable en aplicaciones electrónicas de alta frecuencia [18]. 

La movilidad de los portadores y el nivel de dopado químico son figuras 
de mérito esenciales, y la caracterización de estas propiedades a gran esca- 
la junto con su uniformidad es un prerrequisito para la comercialización de 
dispositivos electrónicos constituidos con grafeno [19]. Por ahora, los disposi- 
tivos electrónicos basados en el grafeno no están implantados en el mercado 
debido a su coste, o porque existen soluciones de otro tipo, con mejores carac- 
terísticas, sin necesidad de usar grafeno. Se necesita poder fabricar grafeno 
más barato, en grandes cantidades y de la mejor calidad posible, para poder 
llegar a una comercialización de éste. Innovadores avances en la investigación 
y aplicaciones industriales del grafeno hacen que la producción en grandes 
cantidades y a un coste razonable sea una etapa necesaria antes de alcanzar 
aplicaciones en el mundo del mercado de consumo [20]. 


Existen aplicaciones en distintos campos, como aplicaciones derivadas de 
su impermeabilidad [21], sensores electrónicos para la detección de moléculas 
químicas y biológicas [22,23], la utilización de su capacidad como supercapa- 
citor [24], y como baterías [25]. También heteroestructuras activas de mate- 
riales 2D impresos con tintas de grafeno y nitruro de boro hexagonal, usadas 
para fabricar transistores de efecto de campo lavables, flexibles e impresos en 
tinta sobre material textil [26]. 


1.3. Características electrónicas del grafeno 


Las propiedades tan inusuales del grafeno, comparadas con áquellas de 
los semiconductores tradicionales, hacen esperar que las características del 
transporte en el grafeno, y en particular, el coeficiente de transmisión para el 
grafeno sea muy diferente al obtenido en los semiconductores tradicionales. 
Por la misma razón, parece natural que la mecánica cuántica tradicional, ba- 
sada en la ecuación de Schródinger no pueda explicar propiedades tan inusua- 
les. Los electrones en el grafeno sin impurezas se presentan como partículas 
ultrarrelativistas y deben ser descritos por la ecuación de Dirac [8,27]. 


El grafeno es un material que consiste en una lámina de átomos de cart- 
bono, de un único átomo de espesor, ocupando los átomos de carbono las 
posiciones de los vértices de hexágonos regulares que constituyen una estruc- 
tura de panal de abejas [16]. Esta estructura bidimensional es la causa de la 
notable flexibilidad estructural del grafeno, que se refleja en sus propiedades 
electrónicas. El átomo de carbono en hibridación sp? proporciona una es- 
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tructura trigonal plana soportada por enlaces C - C tipo dv, que cada átomo 
establece con los tres átomos vecinos, generando así una estructura hexagonal 
de enlaces localizados que proporciona robustez al material. Perpendicular a 
esta estructura se disponen los orbitales atómicos p, no hibridados (uno por 
cada átomo de carbono). Cuya superposición entre átomos vecinos da lugar a 
la formación de una banda de enlaces tipo r cuyos electrones se deslocalizan 
a través de la lámina de grafeno. 


Figura 1.3: Estructura de la red atómica del grafeno. Los círculos rojos (azu- 
les) representan a los átomos de la subred A (B), respectivamente. 


En la figura Fig. 1.3 observamos que la red atómica tiene una estructura 
triangular con una base de dos átomos por celda unidad [14]. Los átomos de 
la subred A(B) están representados por círculos rojos (azules). Un átomo de 
la subred A tiene como vecinos más próximos a tres átomos de la subred B, 
y viceversa. Tomamos como vectores de red 


do 53, V3), n= =D, (1.1) 


siendo a = 0.142 nm la distancia mínima entre dos átomos de carbono próxi- 
mos. Se obtienen los vectores de la red recíproca imponiendo 


d, bj = 210; (1.2) 
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Figura 1.4: Primera zona de Brillouin hexagonal para la red hexagonal. Los 
círculos verdes (morados) representan los puntos K ( K”), respectivamente. 


Y así obtenemos 


hs a, va, 20 co 13). (13) 


Al hacer la bisectriz de los vectores 41, ba, —d1, bh a+ bx, —(a] + b1) 
obtenemos la primera zona de Brillouin, tal y como aparece representada en 
la figura Fig. 1.4, donde los círculos verdes son los puntos K: 


2 IR 


y los puntos K” son los círculos morados: 


=P le): da 


La estructura de bandas del grafeno es conocida desde hace muchos años. 
Su estudio se remonta a 1947 con los trabajos de Wallace [11]. Él utilizó la 
aproximación de electrones fuertemente ligados, los cuales pueden saltar sólo 
entre átomos vecinos, situados a una distancia no superior a dos átomos. 
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Así consiguió describir por primera vez la estructura de bandas del grafeno, 
apareciendo descrita en [14] en su forma actual. 

Supongamos que tenemos N electrones p, del carbono, uno por cada 
átomo de la red, los cuales se extienden por el cristal con pseudomomento 
o vector de onda k, alojándose en dos bandas de energía dadas a partir del 
hamiltoniano de enlace fuerte, por las siguientes expresiones: 


a 
E 
! 


+Hty3+f(k) — t'f(k) (1.6) 





F(k) = 2cos (V3k,a) +4 cos Ch) cos (tea) (1.7) 


La fórmula Eq[1.6] describe la energía de los electrones en la banda de 
conducción (signo +), y en la banda de valencia (signo —) para el grafeno. 
Además t = 2.8 eV es la energía requerida para el salto entre el electrón de un 
átomo de carbono y el átomo vecino más próximo (es decir la realización de 
un salto entre distintas subredes A y B), y t' representa la energía requerida 
para el salto entre el electrón de un átomo de carbono y el átomo vecino más 
próximo en segundo grado (es decir la realización de un salto entre átomos 
de la misma subred bien sea A, o bien B). Vamos a tomar t' = 0,2 -t (en 
cualquier caso, siempre es una fracción propia de t), y a = 0.142 nm es la 
distancia C'— C, como dijimos antes. Observamos en la expresión Eq[1.6] que 
si el segundo término de la derecha fuese despreciable, entonces las bandas 
de conducción y valencia serían completamente simétricas, siendo habitual 
utilizar esta aproximación. 


Existen dos puntos de la red recíproca del grafeno con vectores de onda 
K y K' ( y sus equivalentes, como hemos calculado antes) con alta simetría 


=(Z 27 ) (3 271 ) 
3a 33 a?” $0. 3 30)> 
en los cuales f(K) = f(K”) =- 3, y por tanto según la ecuación Eq[1.6] la 
energía en dichos puntos es nula, es decir E(K) = E(K”) = 0. Dichos puntos 
K y K' se llaman puntos de Dirac. En la aproximación en que trabajamos, 
el espectro de energía resulta simétrico en torno a dichos puntos de energía 
nula. Apareciendo así, para pequeños valores de energía un cono de Dirac en 
torno a cada uno de los seis puntos de Dirac. 

Si consideramos únicamente la región del espacio de momentos próxima 
a uno de los puntos de Dirac, la energía de las bandas de conducción y/o 
valencia pueden ser aproximada entorno a dicho punto de Dirac como k = 
K + q 0 k = K' + q, resultando entonces que la energía se puede expresar 
como: 








(1.8) 
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Figura 1.5: Espectro de energía del grafeno frente a k,a y k,a. Aparecen seis 
conos de Dirac en torno a los seis puntos de Dirac, donde la energía es nula. 








O ql ql 
== ES 2 2 pUES AE 
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donde q es el momento del cristal relativo a un punto de Dirac, up es la 
velocidad de Fermi que viene dada por vp = 3at/2, y cuyo valor aproximado 
es Up > C/300 = 1 - 10% m/s, siendo c la velocidad de la luz en el vacío. La 
velocidad de Fermi es uno de los conceptos clave en el estudio del grafeno. A 
lo largo del tiempo se ha tratado de buscar un valor cada vez más exacto para 
la velocidad de Fermi, resultado de ello son diferentes artículos en los cuales 
se ha estudiado la variación de la velocidad de Fermi bajo ciertas condiciones, 
por ejemplo cuando se cambia la concentración de portadores n [28], tam- 
bién cuando el entorno en el cual se encuentra el grafeno se modifica, lo cual 
lleva a una modificación en la velocidad de Fermi [29], que es inversamente 
proporcional a la constante dieléctrica, alcanzando un valor de - 2,5-10% m/s. 


Introducción 





Despreciando términos de segundo orden en la expresión Eq[1.9], obtene- 
mos una relación de dispersión lineal entre la energía y el momento 


E (la)) = +horlal (1.10) 





La energía de dispersión en la expresión Eq[1.10] recuerda a la energía de 
partículas ultrarelativistas, las cuales son descritas en Mecánica Cuántica 
mediante la ecuación de Dirac sin masa. La relación de dispersión Eq[1.10] 
demuestra que los portadores de carga - electrones y huecos - en un estado 
de la banda se mueven a velocidad 


1 3 
v = +grad, E, = +-at al a 


h 2 ld al 

Observamos en la expresión Eq[1.11] que los electrones en la banda de 
conducción se mueven a la velocidad de Fermi en el sentido del momento, 
mientras que los huecos en la banda de valencia lo hacen a la velocidad 
de Fermi, pero en sentido contrario. En los semiconductores tradicionales 
definidos a partir de la ecuación de Schródinger, la dependencia de la energía 
con el momento era cuadrática, y por ello la velocidad dependía de la energía, 
mientras que en el grafeno la velocidad se mantiene constante, y cuyo valor 
es la velocidad de Fermi. Este comportamiento de los portadores recuerda al 
de una partícula ultrarrelativista descrita por un Hamiltoniano de Dirac [30], 
que puede utilizarse por consiguiente para describir el estado de los electrones 
en el grafeno bidimensional. 

La relación de dispersion Eq[1.10], ha sido observada experimentalmente 
en [31] a través de la dependencia de la masa ciclotrónica de los electrones 
del grafeno con la densidad electrónica. Dicha relación m* x y/n, expresa 
que la masa ciclotrónica es proporcional a la raíz cuadrada de la densidad de 
portadores, en lugar de constante como sucede con una relación de disper- 
sión de tipo parabólico, que es típica de los semiconductores tradicionales. 
La observación experimental de la dependencia de m* con yn evidencia la 
existencia de las partículas de Dirac sin masa en el grafeno. 


=+0p (1.11) 











Los electrones en una lámina bidimensional de grafeno se describen me- 
diante un hamiltoniano efectivo que consiste en dos copias del hamiltoniano 
bidimensional de Dirac, una para momentos cercanos a K y otra para mo- 
mentos cercanos a K”. Las funciones de onda para momentos cercanos a K y 
a K' están relacionadas mediante una simetría de inversión temporal [14,30]. 
Por ello, suele bastar con una única copia del hamiltoniano bidimensional 
de Dirac. En caso de incluir un determinado potencial externo V(x, y), su 
energía potencial U(x, y) será proporcional, y se puede expresar como 


Esopo PARA) UY (1.12) 
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siendo 0, y 9, las correspondientes matrices de Dirac 


dde a (1.13) 


Cada electrón en el grafeno con un momento próximo a uno de los puntos 
de Dirac se describe mediante una función de onda espinorial de dos compo- 


nentes 
E 2, y) 
Víx,y,t) =e *** Pe ) 1.14 


que cumple el sistema de ecuaciones acopladas 


10) 0 

Ep, = Up; — ES =- el (Lo) 
) 0 0 

Epa = Upa — ES + 5): (1.16) 


Para un electrón libre (U = 0) o sometido a una energía potencial cons- 
tante U = U, 4 0 la solución se expresa como una combinación lineal de 
ondas planas espinoriales. Las funciones de onda de la base se clasifican por 
su momento en el cristal k = (k.,, k,) y por su energía E. Los estados acce- 
sibles al electrón corresponden a las funciones de onda: 











ó Ey 1 
W,(x,y, t) EE IS e er Et) y ho +iky ) (L 17) 
[k] 


que están normalizadas en el rectángulo de dimensiones L x W. Dichos 
estados tienen las siguientes propiedades: 


= Un electrón en el estado W,(x,y,t) posee una energía dada por: 


E, =U, + hop |k| (1.18) 


y le corresponde un momento en el cristal k. Se mueve a la velocidad 
de Fermi en el sentido del momento < v, >= -gradx E, = UFI» y 
transporta una densidad de corriente eléctrica j, = —evpYVioV, = 
a en sentido opuesto al momento. Tiene además una helicidad 





definida según Eq[1.24] igual a 3 en el mismo sentido del momento. 
Estos estados constituyen la banda de conducción en el grafeno sometido 


a un potencial externo constante U.. 
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= Un electrón en el estado W_ (x, y,t) posee una energía 


Es == 1] (1.19) 


y le corresponde un momento en el cristal k. Se mueve a la velocidad 
de Fermi en sentido contrario al momento < v_ >= —U Ig» y trans- 


porta también una densidad de corriente eléctrica j_ = —evpW*oV_ = 
a en el mismo sentido que el momento. Tiene también una helici- 
dad definida igual a =>, en sentido contrario al momento. Estos estados 
constituyen la banda de valencia en el grafeno sometido a un potencial 


externo constante U.. 


La estructura periódica de la red en el grafeno impone unas condiciones 
de contorno que limitan el número de momentos k admisibles, y los convierte 
en un conjunto finito. Así, en cada banda hay un número finito de estados y 
un nivel de energía fundamental. 


Además, existe un conjunto de propiedades que son características del 
egrafeno como la fase de Berry, el zitterbewegung y la quiralidad. 


Fase de Berry 

Berry [32] estableció que un autoestado en un sistema cuántico, que se 
transporta lentamente a lo largo de un circuito C' variando los parámetros R 
de su hamiltoniano A(R), adquirirá un factor de fase geométrico expliy(C)]. 
La fase de Berry aplicada al grafeno es la fase que un autoestado expresado 
en la Eq[1.20] adquiere después de que k sea forzado a realizar un círculo 
cerrado a energía constante en torno a un punto de Dirac [33]. 


1 ik.r 1 
Yu = 30 ( 58 ) (1.20) 
La fase de Berry se define como [34]: 
0) 
0 = =$ (0) | an | W(k(t))) dt (1.21) 
E 
La integral tiene lugar a lo largo de un camino cerrado (círculo) a medi- 


da que evoluciona el tiempo. Aplicamos la ecuación Eq[1.21] al autoestado 
Eq[1.20]: 
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En la ecuación anterior Eq[1.22], el segundo y tercer término de la parte 
derecha de la igualdad son nulos al realizar la integración sobre un camino 
cerrado, por ello se reduce a: 


y 1/00 dl e 


Zitterbewegung 

El zitterbewegung (que podemos traducir por "movimiento tembloroso”) 
es un movimiento de vibración ultrarápido alrededor de la trayectoria clásica 
para un electrón que obedece la ecuación de Dirac. En 1930 E. Scródinger 
propuso la existencia de este efecto, al analizar el movimiento de un paquete 
de ondas de electrones en el espacio libre que son solución de la ecuación 
relativista de Dirac. 

De manera más formal, podemos decir que el zitterbewegung (jittery) es 
un movimiento teórico circular o helicoidal de electrones que obedecen la 
ecuación relativista de Dirac, el cual está provocado porque un electrón al 
moverse a velocidades relativistas puede engendrar su propia antipartícula y 
la interacción entre ambas provoca que en la trayectora del electrón aparezca 
un movimiento tembloroso. 

Esta interferencia produce que aparezca un movimiento de vibración ul- 
trarápido a una velocidad muy elevada (del orden de la velocidad de la luz) 
alrededor de la posición del electrón dada, perpendicular a su trayectoria, con 
una determinada frecuencia angular 2m.c*/h o aproximadamente 1.6-10% Hz, 
y cuya amplitud es muy pequeña, del orden de la longitud de onda Compton. 
Este movimiento tembloroso surge porque no es posible localizar la función 
de onda de una partícula relativista en una distancia inferior a la longitud 
de onda Compton. 

En el caso del grafeno, la antipartícula del electrón es un hueco, la longi- 
tud de onda Compton es del orden de 1 nm, y por tanto es posible detectar 
ese temblor utilizando un microscopio de alta resolución [35]. Otra perspec- 
tiva para el zitterbewegung en el caso del grafeno y desde el punto de vista 
de Física de la Materia Condensada, consiste en considerar que es un tipo 
especial de transición interbanda con la creación virtual de pares electrón- 
hueco [36]. También podemos considerarlo como un desorden intrínseco que 
da lugar a una conductividad mínima en el grafeno, incluso cuando la con- 
centración de portadores se aproxima a un valor nulo [36]. 


Helicidad 
La helicidad se define como la proyección del operador momento a lo largo 
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de la dirección del pseudospin [14]. El operador de Mecánica Cuántica para 
la helicidad viene representado por la expresión Eq[1.24] 


> 1 
A (1.24) 
2 |p] 
Los estados W(r) son autoestados del operador h, ya que: 
» 1 
hv(r) = +5 V(r) (1.25) 


donde los electrones tienen helicidad positiva, y los huecos helicidad nega- 
tiva. En la Eq[1.25] representamos los dos autovalores de h. El signo positivo 
o negativo depende de si a coincide con el sentido del movimiento, o es 
contrario a él. 

Un sistema u objeto es quiral si se puede distinguir de su imagen espe- 
cular, es decir al superponerlos no dan el mismo resultado. La quiralidad 
está relacionada con la rotación - dextrógira o levógira - de una partícula 
mientras avanza en sentido lineal. Para partículas sin masa, como sucede en 
el grafeno, helicidad y quiralidad coinciden. Por ello, los estados del sistema 
próximos a los puntos de Dirac K y K' tienen una helicidad o quiralidad bien 
definida. 


1.4. Otros materiales 2D 


Los materiales nanoestructurados poseen tamaños del orden de nanóme- 
tros en al menos alguna dimensión. Una buena clasificación de los materiales 
tiene en cuenta su dimensión (0D, 1D, 2D o 3D), y además su superficie 
libre [37]. Tradicionalmente los materiales eran 3D (tridimensionales), como 
la mayoría de los dispositivos de Electrónica. Sin embargo, en las dos últi- 
mas décadas han aparecido los materiales 2D (bidimensionales) como son los 
materiales laminados de un simple átomo de grosor, los 1D (unidimensiona- 
les) como las nanocintas (nanoribbons), nanotubos y nanohilos, e incluso los 
0D (cero dimensionales) como los puntos cuánticos. La característica prin- 
cipal de los materiales es su dimensionalidad, pues caracteriza la estructura 
atómica del material, y también las propiedades de éste [38]. Así sucede que 
el mismo elemento químico o compuesto puede tener distintas propiedades 
dependiendo de las dimensiones que tenga el material. De momento, a día 
de hoy, las nanoestructuras basadas en carbono están entre los materiales 
más importantes en el campo de la Nanotecnología (tecnología que manipula 
materiales con escala de nanómetros). 
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La realización experimental del grafeno, y sus asombrosas propiedades ha 
propiciado que otros investigadores se lanzasen a la búsqueda de otros mate- 
riales también en dos dimensiones, como el disulfuro de molibdeno (MosS3), 
el nitruro de boro hexagonal (hBN), el estaneno, el siliceno, ... Todo ello ha 
dado lugar al nacimiento de una nueva clase de materiales: los materiales 
2D [39], que se caracterizan por poseer un grosor de una única capa de áto- 
mos. La causa de las extrañas propiedades de los materiales 2D se debe al 
efecto de la dimensionalidad, que da lugar a propiedades diferentes cuando 
el material constituyente, aún siendo el mismo, se encuentra en estado 3D, o 
bien aparece en estado 2D. La obtención de estos materiales en dos dimen- 
siones, depende de la estructura molecular del material de partida en tres 
dimensiones, y de su estabilidad a temperatura ambiente. 


Los materiales 2D son básicamente de dos tipos: 


= Alótropos 2D de un elemento, como son los Xenes que son átomos 
que básicamente pertenecen al grupo IVA de elementos de la tabla 
periódica, y están organizados en una estructura de panal de abeja 
como el siliceno [40], germaneno [41], estaneno [42], borofeno [43] y 
fosforeno [44]. 


= Compuestos consistentes en dos o más elementos, como el nitruro de 
boro hexagonal y los TMD ( Transition Metal Dichalcogenides) que son 
atómicamente muy delgados, del tipo MX , donde un átomo metálico 
de transición M por ejemplo Mo, W es intercalado entre dos átomos 
de calcogenoides X por ejemplo S, Se, Te [45]. 


La familia de los cristales 2D incluye todo tipo de materiales conductores : 
metales ( ej. NbSe2), semiconductores ( ej. MoS» [46]) y aislantes ( ej. nitruro 
de boro hexagonal - hBN) [47]. 

El siliceno, el germaneno y el estaneno forman parte de los cristales de- 
nominados 2D-Xenes, que son similares al grafeno pero con distinto grado de 
pandeo [48]. El nitruro de boro hexagonal es un protector atómico capaz de 
soportar tanto condiciones de calor como de oxidación. Los TMD tienen co- 
mo característica que les diferencia del grafeno, poseer bandgap, y por tanto 
al ser semiconductores, son en principios buenos candidatos para aplicaciones 
en Electrónica y Optoelectrónica, incluso a temperatura ambiente. 


La característica crítica para los materiales 2D es la estabilidad [49]. Hay 
materiales 2D estables a temperatura ambiente en forma de monocapa por 


ejemplo el grafeno, el ya mencionado hBN ( también llamado grafeno blanco), 
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el óxido de grafeno, y MoS2, WS2, MoSez, WSez [50]. Otros son estables 
sólo en atmósferas inertes, o bien aún están por determinar. 


Un campo dentro de los materiales 2D está constituido por las heteroes- 
tructuras de van der Waals. Éstas consisten en planos atómicos aislados de 
diferentes materiales, que son colocados capa a capa de manera precisa según 
la elección realizada, para constituir heteroestructuras [51]. Esto es posible 
ya que existen materiales con enlaces fuertes en el plano, y enlaces débiles 
en la dirección perpendicular, es decir entre diferentes capas, lo cual provoca 
que sea posible aislar esas capas que constituyen ahora, cada una de ellas, 
un material 2D, ya que en la dirección perpendicular únicamente tienen un 
átomo de grosor. 


Una de las finalidades de estos materiales 2D es poder obtener transis- 
tores capaces de operar a temperatura ambiente [52], y así poder reducir el 
tamaño de los circuitos actuales. Los materiales 2D parecen candidatos ade- 
cuados para crear una nueva generación de dispositivos electrónicos, ya que 
permitirán un procesamiento relativamente fácil, acceso directo a los porta- 
dores de carga móviles, una movilidad muy elevada de los portadores de carga 
(a temperatura ambiente y a baja temperatura), alta conductividad termal, 
y un consumo de potencia mucho menor que el asociado a un transistor de 
silicio tradicional. 


A raíz de las investigaciones realizadas sobre el grafeno, surgieron técni- 
cas, como el crecimiento epitaxial [53], para preparar, detectar, transferir y 
caracterizar este material 2D, que posteriormente pudieron emplearse para 
los siguientes materiales 2D que fueron surgiendo. Así con una técnica ya 
desarrollada, el avance fue más rápido en estos últimos. 


La implementación práctica de materiales 2D en bienes de consumo re- 
quiere un buen conocimiento para su realización en dispositivos reales, así co- 
mo durabilidad en condiciones ambiente, y en condiciones extremas. Para ello 
es necesario realizar pruebas de esfuerzo eléctrico, termal, tiempo de vida, 
. ... Además para que el dispositivo funcione correctamente es necesario algún 
tipo de protección para estos materiales 2D. La clave para que la industria 
llegue a utilizar los materiales 2D es una producción a gran escala y a precios 
razonables. 


Entre las posibles aplicaciones están la utilización del MoS» para la reali- 


zación de FETs de radio frecuencias [54]. Además con el rápido crecimiento 
en el mercado de una electrónica de uso y una amplia conectividad a Internet, 
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la necesidad de una nanoelectrónica flexible con alta velocidad de operación 
y bajo consumo se ha convertido en uno de los temas más investigados. Los 
cristales atómicos 2D de van der Waals son prometedores candidatos por 
sus cualidades de miniaturización, elevada flexibilidad mecánica y realiza- 
ción electrónica [55, 56]. 


1.5. Organización de la tesis 


Esta tesis ha sido organizada en los siguientes capítulos: 


En el capítulo 1 hemos comenzado justificando el título de la tesis, comen- 
tado datos históricos sobre el grafeno, y las propiedades más características 
desde el punto de vista electrónico, del grafeno y de otros materiales 2D. 


En el capítulo 2 calculamos el coeficiente de transmisión para una lámina 
de grafeno de anchura indefinida sometida a una o dos barreras de poten- 
cial unidimensional obteniendo una expresión analítica cerrada. La relación 
entre el ángulo de incidencia y el de transmisión en una barrera de poten- 
cial depende de la energía del electrón incidente y de la altura de la barrera 
de potencial Up, apareciendo un ángulo límite bajo ciertas condiciones. El 
coeficiente de transmisión resulta unitario para incidencia perpendicular, y 
también para otros ángulos que verifiquen ciertas condiciones, provocando 
que no pueda existir confinamiento de portadores de carga mediante una ba- 
rrera de potencial eléctrico. El coeficiente de transmisión es mínimo, próximo 
a cero, para ciertos valores de los parámetros, de manera que podemos llegar 
a obtener un pseudogap de energía. El mismo modelo puede aplicarse al caso 
de dos barreras de potencial independientes. Al considerar dos barreras de 
potencial simétricas aparecen unos picos (tipo delta) de resonancia equies- 
paciados dentro del pseudogap. 


En el capítulo 3 hemos construido un modelo para el transistor de efecto 
de campo del grafeno (GFET) que relaciona los potenciales exteriores apli- 
cados en las puertas del dispositivo con los parámetros utilizados en nuestro 
modelo. Determinamos una expresión para la conductancia de los GFET ba- 
jo ciertas condiciones experimentales, en particular para temperaturas próxi- 
mas a cero kelvin, y baja polarización. En principio, consideramos incidencia 
no-perpendicular, para evitar la paradoja de Klein. Inicialmente los porta- 
dores pueden proceder de todas las direcciones posibles y además llegar a 
impactar en la barrera de potencial con la misma probabilidad. A continua- 
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ción, tratamos el caso de direcciones privilegiadas, es decir que tienen más 
probabilidad de proceder de esas direcciones que de otras, lo cual simula- 
mos mediante funciones de ponderación tipo campana de Gauss o tipo delta. 
Físicamente, esta dispersión que puede suceder al considerar direcciones pri- 
vilegiadas, la asociamos a muestras con pocas o bastantes inhomogeneidades. 


En el capítulo 4 realizamos una extensión del capítulo 2, aplicando múlti- 
ples barreras de potencial a una lámina de grafeno. En vez de obtener una 
expresión analítica cerrada, utilizamos el Método de la Matriz de Transferen- 
cia (MMT). Analizamos el caso de múltiples barreras de potencial de alturas 
iguales o diferentes. Como aplicación, también consideramos la barrera de 
potencial escalonada que permite simular cualquier potencial que varíe sua- 
vemente, aplicándolo a continuación a una función gaussiana. 


En el capítulo 5 realizamos otra extensión del capítulo 2, con el fin de 
mejorar y refinar nuestro modelo, para lo cual consideraremos que sobre la 
lámina de grafeno además de influir las barreras de potencial asociadas a 
las puertas superiores Vr, también lo hace el campo eléctrico asociado a la 
diferencia de potencial entre surtidor y drenador Vps, dando lugar a un po- 
tencial efectivo lineal en ciertas regiones de la lámina de grafeno. Lo anterior 
cambia el potencial y por tanto el Hamiltoniano, lo que nos lleva a construir 
de nuevo los estados electrónicos en los que es necesario introducir funciones 
de Hermite e hipergeométricas. Al someter dos barreras de potencial a campo 
eléctrico aparecen nuevos efectos. El orden que ocupan las distintas barrreras 
de potencial influye de tal forma que la segunda barrera de potencial puede 
llegar a quedar enmascarada. 


En el capítulo 6, consideramos una nanocinta de grafeno (GNR) de longi- 
tud infinita pero de anchura finita, que sufre a continuación un estrechamien- 
to. La conductancia básicamente viene dada a partir del número de modos 
que se propagan en la parte estrecha de la cinta, para lo cual comparamos la 
conductancia de la cinta que ha sufrido estrechamiento con la asociada a una 
cinta cuya anchura coincidiese con la parte más estrecha. En esta reducción 
progresiva de las dimensiones, consideramos una nanoconstricción rectangu- 
lar en la nanocinta de grafeno, con un doble propósito: comparar por una 
parte nuestros resultados analíticos con los que para una cinta idéntica ob- 
tienen mediante el método tight-binding nuestros colegas de la Universidad 
Complutense de Madrid y, por otro, con las medidas que realizan en el centro 
de I4D+i de la Universidad de Salamanca nuestros colegas experimentales 
sobre una nanoconstricción en una nanocinta construida expresamente para 
ello. 
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Las ventajas que presentaba el método tight-binding de ser escalable y 
poder definir muy fácilmente la forma de los bordes, junto con una inesperada 
lenta convergencia del coeficiente de transmisión analítico hizo aconsejable 
finalmente utilizar este último método para realizar la comparación con las 
medidas experimentales. 
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Capítulo 2 


Cálculo analítico del coeficiente 
de transmisión a través de una 
y dos barreras de potencial 


2.1. Introducción 


La transmisión electrónica a través de una o dos barreras de potencial, con 
distintas formas, en grafeno ha dado lugar a trabajos que tratan el fenómeno 
en modo diverso. El estudio de una barrera de potencial en una monocapa 
de grafeno para una partícula ultrarrelativista (sin masa) fue estudiado por 
Katsnelson y col. en [57], donde caracterizan el coeficiente de transmisión 
mediante una fórmula simplificada. Nguyen y col. en [58] utilizan el méto- 
do standard de la matriz-T' para estimar numéricamente la probabilidad de 
transmisión en una barrera unidimensional con formas diversas. Los fenóme- 
nos de transporte en el grafeno han sido también estudiados por diversos 
autores en [59-64] donde se justifica el efecto túnel de Klein y se obtiene el 
coeficiente de transmisión a través de una barrera de potencial. Adam y col. 
en [65] han considerado el papel de las impurezas cargadas y la concentra- 
ción en que éstas aparecen, que originan la dispersión de portadores de su 
dirección de propagación, y por tanto afecta de manera importante al trans- 
porte. La transmisión electrónica resonante a través de una doble barrera 
se ha considerado en [66,67], y en [68] por medio de las funciones de Green 
asociadas. 


Pretendemos estudiar el comportamiento de algunos dispositivos electróni- 
cos cuando se sustituyen los semiconductores convencionales por grafeno mo- 
nocapa. En particular, nos interesa estudiar el comportamiento de diversos 
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transistores de efecto de campo FET, que mediante un arreglo experimen- 
tal adecuado y haciendo uso de diferentes potenciales de puerta, someten 
al grafeno a potenciales externos que pueden ser simulados de distintos mo- 
dos. En la aproximación más sencilla, consideramos que dicho potencial es 
asimilable a una o dos barreras de potencial. Los electrones de conducción, 
caracterizados en la aproximación continua de baja energía como partículas 
ultrarrelativistas de Dirac, se transmitirían a través de los estados disponibles 
del grafeno sometido a estos potenciales. Desde los trabajos de Landauer [69] 
es conocido que el transporte coherente de corriente en estos dispositivos 
está determinada por la probabilidad de que un electrón pueda transmitir- 
se a través de él. Nuestro primer objetivo es obtener una expresión analítica 
exacta del coeficiente de transmisión T para estos potenciales, que servirá pa- 
ra discutir la influencia de las magnitudes que intervienen en el fenómeno: 
energía y ángulo de incidencia de los portadores, altura(s) de la(s) barrera(s), 
longitud y separación de éstas y construir después, en el régimen lineal, la 
conductancia del dispositivo. 





2.2. Coeficiente de transmisión para una ba- 
rrera de potencial 


Consideramos una lámina de grafeno monocapa idealmente infinita some- 
tida a una barrera de potencial finita. La lámina se considera lo suficiente- 
mente extensa tanto en la dirección x como en la dirección y, como para poder 
ignorar los efectos de borde. Esta aproximación desacopla completamente el 
movimiento electrónico, y convierte en idéntico el transporte de los electrones 
en los dos valles del grafeno, de modo que, para construir la conductancia del 
dispositivo o la intensidad de corriente, basta con considerar uno de ellos, el 
valle K por ejemplo, y multiplicar el resultado por 2. Una barrera de energía 
potencial unidimensional, representada en la figura Fig. 2.1, en el plano x— y 
está descrita por: 


0 si <a 
Ulx,y)=x% Un si a<x<a+L (2.1) 
O si a+L<xzx 


Ya que la energía potencial y el potencial eléctrico están relacionados me- 
diante U = q: V, siendo q la carga del portador que será —e para el electrón 
y +e para el hueco, donde e representa la carga elemental. Por simetría, 
podíamos haber puesto la barrera de potencial centrada en el origen. Sin 
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0 a a+L 0 


Figura 2.1: Barrera de potencial unidimensional de altura Uy y anchura L. 











Figura 2.2: Ángulos de incidencia Q y de transmisión 0, al impactar sobre la 
barrera de potencial. 
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embargo, se ha desplazado para tener un resultado más general, directamen- 
te aplicable a dos barreras de potencial. Naturalmente, nuestros resultados 
demuestran de manera expresa la independencia del fenómeno frente a la 
posición de la barrera. 
El hamiltoniano a resolver es entonces: 
a , £ A U(z, y) Pe 
H =0pl0P, +0yPy FU (a. y)1)= Ur ( P, Ula, y) ) (2.2) 


donde P_ = P, — Re y P, = P, + Es La solución del hamiltoniano 
anterior adopta una forma funcional diferente en cada una de las regiones 
determinadas por el potencial. Para construir los estados posibles conviene 
tener en cuenta algunos aspectos: 


1.- El hamiltoniano es independiente del tiempo, y por tanto la energía 
es una constante del movimiento. Así, los estados deben tener una misma 
energía fija E en todas las regiones. 


En las regiones extremas donde el electrón es libre, regiones / y 111, el 
módulo del momento debe ser | k |= k = Se con signo positivo para la 
banda de conducción, o negativo para la banda de valencia. Por otra parte, 
en la región de la barrera, región //, el momento debe ser diferente, puesto 
que debe tener en cuenta la energía asociada a la barrera de potencial. Su 
valor es | Y [E ed según sea el signo de E — Up. Esto significa que, 
para el caso considerado, en las regiones donde la energía potencial es nula 
(I y III) pueden propagarse indistintamente estados con energía positiva 
o negativa, es decir tanto de la banda de conducción como de la banda de 
valencia. En la región 11 en la cual tenemos barrera de potencial, en cambio 


hay que considerar varios casos: 





= Cuando E > Uy > 0, sólo existe un valor físicamente aceptable para 
el momento (en módulo) si se elige el signo positivo en la ecuación 
q = rd que corresponde a un estado de la banda de conducción. 
En este caso, los electrones se propagan en todas las regiones por la 


banda de conducción. 





= Cuando Uy > E > 0, debe elegirse en cambio el signo negativo en 
la ecuación q = o correspondiendo a un estado de la banda de 
valencia. En este caso, los electrones que se propagan por la banda de 
conducción en la región libre (zonas extremas) sólo pueden propagarse 
en la zona intermedia por la banda de valencia. Por tanto se comportan 


como huecos. Existe así un cambio en el comportamiento, ya que se 
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comportan como electrones en las zonas extremas (1 y 111), y como 
huecos en la zona intermedia (11). 


= Cuando E < 0, los electrones sólo pueden propagarse por la banda de 
valencia en todas las regiones. Es decir, se comportarán como huecos. 


2.- El potencial es independiente de y, por tanto la componente k, del 
vector de onda debe ser la misma en las diferentes regiones del espacio. 


Sin embargo, k.,. toma distintos valores dependiendo de la zona en que se 
2 


encuentre. En las zonas extremas / y IT es k, = E, / 537 — ko, mientras que 
F 


en la zona intermedia // es q4 = +, pea — ki. Donde el signo positivo 
F 


está asociado a ondas incidentes propagándose hacia la derecha, y el negativo 
a ondas reflejadas propagándose hacia la izquierda. Tanto k.,. como q, pueden 
ser reales o imaginarios, dando lugar, respectivamente a ondas propagantes 
o no propagantes, eventualmente evanescentes. 


Este comportamiento puede interpretarse en términos de las leyes de Snell 
de reflexión y refracción de ondas. Un electrón incidente con un ángulo 9d, en 
cada región deben cumplir: 


= El ángulo de incidencia (formado por el rayo incidente y la normal) d 
es igual al ángulo de reflexión d” (formado por el rayo reflejado y la 
normal). 


= El ángulo de incidencia Y y el ángulo de transmisión O (formado por 
el rayo transmitido y la normal) verifican: 


k, = ksin 9 = gsin O, | E | sind =| E— Uy | sinO (2.3) 


= Para determinadas energías, existe un ángulo de incidencia límite tal 
que la partícula se transmite a 90%. Para ángulos de incidencia supe- 
riores, la componente x del momento transmitido se torna imaginaria, 
lo que origina que las ondas planas de los electrones incidentes se con- 
viertan en ondas evanescentes, descritas por exponenciales reales. El 
ángulo límite está establecido por la condición 


0 
1l-= j 
- 2.4) 


$, = arcsin 
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De este modo, los electrones con energía inferior a 0.5U¿g no pueden 
sufrir reflexión total en la barrera y se transmiten con un ángulo de 
transmisión inferior al de incidencia, es decir, se acercan a la normal; 
los electrones con energía igual a 0.5 Uy se transmiten con un ángulo de 
transmisión igual al de incidencia, y los electrones con energía superior 
a 0.5 Up sólo pueden transmitirse por medio de una onda plana si el 
ángulo de incidencia $ está comprendido entre los valores —9, y Y; 
correspondientes al valor de su energía. Cuando lo hacen, se alejan de 
la normal. En particular, para una energía E = Up, el ángulo límite es 
nulo, por lo que sólo hay ondas planas transmitidas en el interior de la 
barrera para incidencia normal. Todo este comportamiento se recoge 
en las Fig. 2.3 - 2.4. 





positivo 


BIS NIN 





La 
o 





negativo 











E/Uo 


Figura 2.3: Representación del ángulo límite $, frente a E/Uy, tanto para 
ángulos positivos como negativos. Observamos que no existe ángulos de inci- 
dencia limite para valores inferiores a 0.5, y además para 1 alcanza su valor 
mínimo. 


En la figura Fig. 2.3 se representa el ángulo límite O, frente a E/Uo, tan- 
to para ángulos límite positivos como negativos. Se observa claramente 
en dicha figura que para valores inferiores a E/U¿ = 0.5 no existe ángu- 
lo de incidencia límite. Éste decae hasta que alcanza el valor E/Up = 1, 
a partir del cual comienza a crecer de nuevo. Del mismo modo, en las 
figuras Fig. 2.4 se observa la relación entre el ángulo de incidencia y el 
de transmisión para diferentes energías, en (a) para energías inferiores 
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a 0.5 Us, y en (b) para energías superiores a 0.5 U¿. Cuanto menor es 
E/Ugp, menor será el ángulo de transmisión con respecto al de incidencia 
y viceversa. 


De acuerdo con las consideraciones anteriores, la función de onda de 
energía E para el potencial en la ecuación Eq[2.1] puede ser expresada me- 
diante el espinor de Dirac: 


ehuy- RD ikyx 1 —ik.,x 1 
W;¡(x, y, t) = ABLW EAN 5H + re Ñi ek (2.5) 
k k 


ehuy= RD 1Quí 1 —1iqut 1 
A E 


ee 1 

Wip (x, y, t) = nada q TL ( sha ) (2.7) 
k 

(2.8) 


donde se ha denotado s = sign(E), sí = sign(E — Up), kz = kz + iky, 
+ = Qu Eldky, k= y k¿+k5 y q = yq + kz. Además se ha normalizado 
la función de modo que haya un electrón incidente en el sistema por unidad 
de superficie. Se ha considerado también que los electrones que inciden en el 
sistema desde la izquierda, pueden reflejarse o transmitirse en las zonas / y 
TI, pero en la región /1I, una vez atravesada la barrera, sólo quedan partículas 
propagándose hacia la derecha. Además, se ha supuesto también que el ángu- 
lo de incidencia es inferior al ángulo límite H < /. De no ser así, tendríamos 
en cuenta lo anterior. 








En las heteroestructuras semiconductoras los electrones obedecen la ecua- 
ción de Schródinger que es de segundo grado, y la continuidad de la densidad 
de corriente en la interfaz entre dos regiones contiguas del potencial, lo cual 
implica dos condiciones de contorno : la continuidad de la función de onda y 
la de su derivada en la mencionada interfaz [70]. En el grafeno, la ecuación de 
propagación es la de Dirac, que es de primer grado, por lo que será suficiente 
con imponer la continuidad de las funciones de onda, que ahora implica la 
continuidad de la corriente. En nuestro caso, las funciones de onda deben 
verificar la condición de continuidad de la función de onda en la frontera 
entre zonas adyacentes. De este modo, se debe cumplir en todo instante: 


Y ;¡(a, y,t) => Y ¡¡(a, y, t), (2.9) 
Vi¡(a + L, y, t) Vipi(a+L,y,t) (2.10) 


2 


Coeficiente de transmisión 








IS 









































Ia 
PIN + 
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Figura 2.4: Ángulo de transmisión O frente al ángulo de incidencia Y para 
diferentes valores de E/Uj. La línea gris representa O = Y. A medida que 
disminuye E/Up, menor será el ángulo de transmisión con respecto al de 
incidencia; mientras que a medida que aumente, sucederá lo contrario. 
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tal y como habíamos obtenido en nuestro artículo [71], lo cual proporciona 
los coeficientes 


e —e Pako kk (k— ss1q) sin(Lq;) (2.11) 
— K(kaqw cos(Lq) + ¿(k? — ss1kq) sin(Lqu)) 


ia(kx—Qqu) —4Lqx fo k 
E e e Y alkq_ + ss19ky) (2.19) 
2k(k qu cos(Lqs) + 1(k¿ — ss1kq) sin(Lqs)) 


10(ka +90) piLque 407 
B- e e Ko (lay 5519k+) (2.13) 
2k(k2 q; cos(Lq,) ST ¿(kz ES ss1kq) sin(Lq,)) 


ekolk, q, 
= (2.14) 
kz qu cos(Lq«) + ¿(k5 — ss1kq) sin(Lq.) 

















En el tratamiento de la corriente suele interesar obtener el coeficiente 
de reflexión y el de transmisión en una dirección determinada. En nuestro 
caso en la dirección x, de acuerdo con la definición del potencial dada en 
la ecuación Eq[2.1]. El coeficientes de reflexión R, que es la probabilidad 
de que una partícula sea reflejada, y el coeficiente de transmisión T', que 
es la probabilidad de que una partícula sea transmitida, vienen definidos 
por el cociente entre la densidad de corriente reflejada o transmitida, y la 




















A A 

densidad de corriente incidente: R = Gel e . LE O = 
Win. trans 

vw? 00 Vine 

En nuestro caso, tales coeficientes son 

k2(q — ss¡k Y? sin*( Lg. 
k2q2 cosA(Lq.) + (k? — ss1kq)? sin*(Lq,) 
2 2 
A ade (2.16) 





— K2g2 cos(Lq.) + (K2 — ss1kq)? sin" (Lqo) 


Los coeficientes de reflexión R y de transmisión T' deben obedecer la ley 
de conservación de las partículas; es decir R+T = 1. El coeficiente de trans- 
misión, que coincide con el descrito en [14], es simétrico respecto al ángulo de 
incidencia 7 (49) = T(—¿), lo cual significa que es independiente del sentido 
de incidencia de los electrones. En particular, para un ángulo de incidencia 
nulo (es decir, incidencia perpendicular a la barrera de potencial), el coefi- 
ciente de transmisión es la unidad. La barrera es transparente para todos 
los electrones que incidan perpendicularmente a ella independientemente del 
valor de su energía. Este fenómeno ya conocido para la ecuación de Dirac se 
le conoce como paradoja de Klein [15,57,59]. Esto no ocurre en los materiales 
semiconductores para electrones no relativistas, en los cuales el coeficiente de 
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transmisión para energías inferiores a la altura de la barrera de potencial es 
siempre menor que la unidad [72]. La barrera también es transparente para 





aquellos electrones con q. = FF, es decir cuya energía y ángulo de incidencia 
verifiquen 
nahup y” , es 
7 =(E-— Up) - E sinó (2.17) 


Debido al coeficiente de transmisión unitario tanto para incidencia nor- 
mal, como para los ángulos incidentes que verifiquen la condición expresada 
en la fórmula Eq[2.17), no permite el confinamiento de los portadores de car- 
ga mediante barreras de potencial eléctrico. Ya que al menos para ciertos 
ángulos, los portadores de carga están deslocalizados en el material. Esto es 
así, mientras no se abra un intervalo prohibido (gap) mediante algún pro- 
cedimiento [73]. Se ha considerado hasta aquí, al analizar las características 
de la barrera de potencial, que los electrones inciden con un ángulo inferior 
al ángulo límite. Si lo hacen con un ángulo superior, el análisis sigue siendo 
válido, y el único cambio que se ha de realizar en las fórmulas anteriores 
es cambiar q, por 1q», lo que transforma las ondas planas e+*'%" en ondas 
evanescentes eF1%*. Así, en el interior de la barrera, el transporte de elec- 
trones no contiene términos oscilatorios sino exponenciales reales. En cuanto 
a los coeficientes, se ha de cambiar cos(Lq.) por cosh(Lq,) y sin(Lq.) por 
¿sinh(Lq,). En particular, el coeficiente de transmisión resulta en este caso 








Edo 


x 


T pro 
k2g2 cosh*(Lq.) + (k7 — ssikq)? sinh*(Lq,) 





(2.18) 


Así, el coeficiente de transmisión dependiente del ángulo de incidencia, 
con respecto al valor del ángulo límite, se puede expresar explícitamente: 


"= Sió< 01 





E > (2.19) 














Ti= (2.20) 
(Y (7 

2 2 
(E) sin? dea) 
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Estos valores exactos, se reducen para valores elevados de la barrera de 
potencial a los valores aproximados para el coeficiente de transmisión con- 
siderado en [14]. Sin embargo, hay que mencionar que para valores no tan 
altos, la diferencia entre ambas expresiones puede ser significativa. 


El coeficiente de transmisión depende de cuatro variables, dos asociadas a 
la partícula incidente: su energía E y su ángulo de incidencia Q, y dos asocia- 
das a la barrera de potencial : su altura Uy y su anchura L. Según aparecen 
en las ecuaciones Eq[2.19, 2.20], las variables anteriores se pueden agrupar 

E 


en la energía relativa Do adimensional, el producto L - Uy que se expresa en 


eV - nm, y el ángulo de incidencia d. 


En la figura Fig. 2.5 representamos el coeficiente de transmisión frente al 
ángulo de incidencia para los parámetros L - Uy = 5, 10, 15 y 20 eV- nm. 
En (a) para E/U¿ = 0.3 y en (b) para E/Uy = 0.8. Se aprecia en ambos 
que, para una energía fija, el coeficiente de transmisión tiende a decrecer 
más o menos abruptamente a medida que aumenta el ángulo de incidencia. 
Cuando el producto £ - Uy es suficientemente grande empiezan a aparecer las 
resonancias descritas en la fórmula Eq[2.17]. Se observa en (a) que los valores 
de energía E /Uy < 0.5 se transmiten sin decaer abruptamente, cualquiera que 
sea el ángulo de incidencia. Cuando E/Uy > 0.5 en cambio, la existencia del 
angulo límite dependiente de la energía es evidente como observamos en (b). 

En la figura Fig. 2.6 representamos el coeficiente de transmisión frente 
a la energía relativa E/Uy para los parámetros L - Uy = 0.5, 1, 1.5 y 2 eV- 
nm. En (a) para fp = 7/12 y en (b) para $ = 7/4. En ambas figuras se 
aprecia que, para un ángulo de incidencia fijo, el coeficiente de transmisión 
es simétrico en torno al eje E = Uy cos”? o, lo que significa que estados de la 
banda de valencia con energía E = (cos”? p — €)Uo y estados de la banda de 
conducción con energía E = (cos”? H+€)U, tienen la misma probabilidad de 
transmitirse. 

Siempre que el ángulo de incidencia sea suficientemente grande, el coefi- 
ciente de transmisión empieza a decrecer desde T' = [14+[LU¿/(hvp)]? tan? p]]* 
en Ej = (1 +sin|o|)*Uo pasando por un mínimo de transmisión T' = 
cosh”*[LU¿tanó/(hvp)] en E = cos ?4Us y vuelve a recuperar su valor 
inicial en Ez = (1 — sin|H))-*Uj. En la práctica, aparece un intervalo de 
energía de AE = 2U,sin | d | cos”? 4 en torno al eje de simetría, para el 
cual la probabilidad de transmitir un electrón decrece prácticamente a cero, 
es decir, es una zona de baja transmitividad que de aquí en adelante deno- 
minaremos pseudogap. Puesto que sin |H| está comprendido entre 0 y 1, por 
tanto Ey puede variar entre Uy y Up/2 mientras que Ez puede variar entre Uy 
y co. Al representar la variación relativa de ese pseudogap de energía AE/U% 
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Figura 2.5: Coeficiente de transmisión en función del ángulo de incidencia 
siendo E/Uj = 0.3 (a) (0.8 (b) ) , para distintos valores de L - Uy (eV-nm). 
En (a), al ser E/Upy inferior a 0.5 el coeficiente de transmisión decae suave- 
mente para todos los ángulos, mientras que en (b), al ser E/Uj superior a 
0.5 decaerá abruptamente. 
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Figura 2.6: Coeficiente de transmisión frente a la energía relativa E/Uy para 
un ángulo de incidencia fijo y = 7/12 en (a) (7/4 en (b)), y distintos paráme- 
tros L- Up (eV- nm). Existe un eje de simetría para E = Uy cos”? f, en torno 
al cual para ángulos suficientemente elevados aparece un pseudogap. 
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frente al ángulo de incidencia ¿ en la figura Fig. 2.7, observamos para ángulos 
pequeños la variación del pseudogap es muy pequeña, mientras que en torno 
a 37m/8 se presenta una variación muy acusada, tanto mayor cuanto mayor 
sea el ángulo. 
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Figura 2.7: Anchura relativa del pseudogap AL /U, frente al ángulo de inci- 
dencia. Observamos una variación muy acusada para ángulos elevados. 


Puesto que en muchos experimentos se controla la altura de la barrera de 
potencial por medio de un potencial de puerta aplicado al dispositivo, resulta 
interesante ver cómo se comporta el sistema para valores diferentes de dicha 
energía potencial que, resulta fácilmente sintonizable al ser una variable que 
se controla externamente. En la figura Fig. 2.8 se representan las líneas de 
nivel del coeficiente de transmisión T para valores de 0.1 (amarillo), 0.3 (na- 
ranja), 0.5 (rojo), 0.7 (verde) y 0.9 (azul) frente a la anchura L y a la altura 
Uv de la barrera, teniendo como parámetros en el panel izquierdo E = 0.1 
eV y Y = 1/12 y en el panel derecho E = 0.3 eV y 6 = 7/6. Observamos 
que a partir de una cierta anchura L, que dependerá de los valores de los 
parámetros considerados E y /, aparece un pasillo asociado al pseudogap que 
está centrado en torno a un valor de Uy semejante a la energía incidente £, 
y que es tanto más ancho, cuanto mayor sea el ángulo de incidencia ¿. 


El caso de un pozo de potencial Uy < O también puede ser tratado de 
modo análogo a la barrera de potencial. Esencialmente se encuentran los 
mismos resultados. Los coeficientes r, A, B y t son formalmente los mismos 
que los obtenidos para la barrera. En particular, el coeficiente de transmisión 
T es idéntico al que hemos dado en la en la fórmula Eq[(2.18] para la barrera. 
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Figura 2.8: Líneas de nivel correspondientes a T= 0.1 (amarillo), 0.3 (na- 
ranja), 0.5 (rojo), 0.7 (verde) y 0.9 (azul). En el panel izq. para E = 0.1 
eV y $ = 1/12, y en el dcho. para E = 0.3 eV y y = 7/6. Observamos un 
pasillo asociado al pseudogap, centrado en torno a U¿y - E, y cuya anchura 
aumentará con el ángulo de incidencia d. 


La única diferencia que existe es el valor negativo de Up. Nótese que el coe- 
ficiente de transmisión cumple 7(E, Up) = T(- E, —Ub), que al aplicarlo a 
los electrones del grafeno, significa que la probabilidad de transmisión de un 
electrón de la banda de condución a través de una barrera de potencial de 
altura determinada, es la misma que la de un electrón de la banda de valen- 
cia con energía opuesta, a través de un pozo potencial cuya profundidad sea 
igual a la altura de la barrera anterior. Como consecuencia, las propiedades 
de conducción eléctrica del grafeno dopado n y el dopado p son idénticas, 
salvo en lo relativo al signo de los portadores. 





2.3. Coeficiente de transmisión a través de 
dos barreras de potencial 


Consideramos ahora dos barreras de potencial genéricas en la figura Fig. 2.9, 
con diferentes alturas U¡ y Uz y anchuras L¡ y Lo, separadas una distancia 
arbitraria D. 


39 


Coeficiente de transmisión 
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Figura 2.9: Dos barreras de potencial de alturas U, y Uz, y anchuras L¡ y La 
respectivamente, separadas una distancia D. 


El potencial de dichas barreras es ahora: 


O si <a 
U, si a<x<a+L; 
Ule.y)=x 0 si a+Li<x<a+Li+D (2.21) 
Us si a+ Li¡+D<zx<a+Li+D>wL) 
0 si a+Li+D+L<xzx 


Podemos construir la solución a nuestro potencial mediante la siguiente 
función de onda, troceada en las mismas regiones que el potencial: 


i(kyy—Ft) 
E e h ikyx 1 —ik7x 1 
Sa! (e Cos Jrs) | 
eu 01 y T 1 —1Q1 q T 1 
Wir(x, y,t) e e Aje ij q + Bye d q 
v2LW Si $174 


Ce Ukyy= Ft) ik 1 ¿e 1 
k k 
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a , e Ukyy- Ft) cd 1 e B sk 1 
vz, Y, ) —Ñ 2L 29€ So e 29€ So = 
(kyy t) 
e h ik 1 
= e 2.22 
Wy(x, y, t) 2LW e ( ste ) ( ) 


donde s = sign(E), s; = sign(E — U;), kz = ky E1ky, k = yk + kg, 
is = Qi, Elky Y Qi = 4/0í, + Ko: 








De manera semejante al caso de una sola barrera de potencial, se ha con- 
siderado que las partículas se propagan de izquierda a derecha. En cada zona 
pueden reflejarse y transmitirse, excepto en la región Y, donde sólo quedan 
partículas propagándose hacia la derecha, una vez atravesadas todas las ba- 
rreras de potencial anteriores. Al igual que en el caso anterior, se conservan 
la energía del electrón E y el momento perpendicular k, . La componente x 
del momento, se ajusta en cada región de modo que resulta +k.,, en el vacío 
y +q;, en cada una de las regiones de las barreras. 


Las propiedades de reflexión y transmisión (que vimos en el caso de una 
barrera), se cumplen también en cada una de las interfaces del sistema ahora 
considerado. Con la particularidad de que, ahora existen dos ángulos límite, 
asociados cada uno de ellos a una barrera de potencial. Para un ángulo de 
incidencia superior al ángulo límite en cualquiera de las dos barreras de 
potencial, el electrón se propaga en el interior de la misma por medio de 
ondas evanescentes. Pudiendo llegar a suceder en ninguna, en una o en ambas 
barreras de potencial. Los ángulos límite son: 


0; 
d;, =arcsin 1 — El 12 (2.23) 








Utilizando la condición de continuidad de la función de onda a través de 
los límites entre regiones contiguas, se debe cumplir: 


Y ¡¡(a, y, t), 
= Viri(a + Li, y,t), 
= Viy(a+L¡+0D,y,t), 
= Vy(a+ L¡+D+Loa,y,t), 


V¡(a, y, 

Vir(a+ Li, y, 

Vijr(a + Li +0D, y, 
Viy(a+ L¡ + D+ Lo, y, 





t) 
t) 
t) 
t) 
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lo que proporciona 


¿aka (A 
Z1+1i%3) k 


( (i(lgass, =- k) sin(Loq2,) — kz 02, cos(L2q2,)) (sk — s¡q1) sin(L1q1,) 
pete (a — kquss¡) sin(Liq1,) — ka41, cos(L1q1,)) 
(ls — 4252) sin(Lag»,) |. (2.28) 
elalka—q12)=iL1015 
(ZA +IZ) 2k 
( (kysiqu +skqu_) (aga, cos(La2q9,) + ¿(k; — ss2kq2) sin(Lag»,)) 
+ ek k, (ss kqu — k,q1_)(sk — s2q2) sin(L2q0, ) + (2.29) 
caleta Hilado Lo, 
TEREZN 2 
( (skqu, — k+51q1) (Tag, cos(Lo2q,) + ¿(k; — 582kq2) sin(Lag»,)) 





pe 


A = 





Bb; = 




















+ eee + ss1kq1)(s2q2 — sk) sin(Loq, ) ) (2.30) 
e Like 
O, = +53 kcMo 
A (21 +1Z3) is 
(Tag, cos(L2q2,) + ¡(ky — ss2kg2) sin(Lag»,)), (2.31) 
e Qa+2D+L1)ko po 
D¡ = 2 1 kk —sk)sin(L 00 
18 (Z, + 1Zo) y Ha yk+(s2q2 — sk) sin(Loq2,), ( ) 
77 e ¿(La L2)92. Hi(a DNkz 92, ) k2q1 L L 
= vita 2.33 
2 S WAETTA 2h (s22 pas q), ( ) 
(L1+L2)q2, + (a+D)(ko+92.7) 2 
e xq se 
Ba TESTA ma (92, sk — ky599), (2.34) 
; kg 42 
t = e MollitLo) Tatlo 12 : 2.35 
Li + 1Zo ( ) 


Zi = k2q1,42, cos(L:q1,) cos(L2q»,,) 
= (sisggo + 10 - skki(s1q1 + 8249)) sin(L1q1,) sin(L2q, ) 


+k5 (51520103 — sk(s1q1 + 82q2) + 12) cos(2Dk.;) sin(L1q1,) sin(L2q2,) 
(2.36) 
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Za = ki (5129100 — sk(s1q1 + 822) + e sin(2Dk.) sin(L1q1,) sin(L2q2, ) 
—k (02. (5s1/g - k5) cos(L2q2, ) sin(L1q1,) 


+41, (ssokq2 — ky) cos(L;q1,) sin(L2g2,)) (237) 


Y el coeficiente de transmisión resulta, explícitamente dependiente de las 
características de la barrera U¡, Uz, L¡, La y D, así como de las características 
del electrón incidente E y (;: 


T=1-R=tf =T|E,U,,Uz, Li, La, D q] = 


== [eos a] cos 7] + (0,0 tan? Hcos (E = 

















: Lje : Loe 
U, — Ecos? $ U, — ca [e | Sn [+2] 0 


cos Q cos Q €1 €) 








: Lie : Loae 
2DE cos id sn [e] en [52 | V 


(0, U) tan? ¿sin hor 


€1 €2 








5d Lje 1 
U,— Ecos*9 ] ES ses [ee] 














cos Y hup €1 
: Loco | 
Us — Ecos? $ [El SIN [es e 
> 2.38 
cos Y EE hup € ( ) 


siendo 





a = «[E,U, 0] =y(E-U)?-Etsidó,  ¿i=1,2 (2.39) 


Casos particulares de las dos barreras de potencial en la ecuación Eq[2.21], 
que según se demuestra mediante un cálculo directo, se reducen apropiada- 
mente a: 


= Cuando U;z; = 0 al caso de una sola barrera de potencial de longitud 
L,. 


= Cuando D =0, y U¡, = Uz a una sola barrera de potencial de longitud 
L;¡ + Lo. 


El coeficiente de transmisión es invariante frente al intercambio de las va- 
riables correspondientes a las dos barreras, y resulta por tanto independiente 
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del orden en que se coloquen las dos barreras de potencial. O lo que es lo 
mismo, el resultado es independiente de que la incidencia tenga lugar desde 
la izquierda o desde la derecha. 


También pone de manifiesto el fenómeno de acoplamiento entre ambas 
barreras de potencial, debido a la aparición no trivial del parámetro D. Si 
se considera que las barreras están muy alejadas, D > Lj, La entonces las 
oscilaciones sin(2Dk.,) y cos(2Dk.,) se hacen muy rápidas y pueden ser sus- 
tituidas por su valor medio que es nulo. En estas condiciones desaparece el 
acoplamiento entre barreras y el coeficiente de transmisión se reduce al pro- 
ducto de los coeficientes de dos barreras independientes T' = T| - Ta. 


También en este caso se cumple la simetría del coeficiente de transmisión 
respecto al ángulo de incidencia T(4) = T(—0) y es la unidad para un ángulo 
de incidencia nulo (paradoja de Klein). Es decir, las dos barreras de potencial 
son transparentes para todos los electrones que incidan perpendicularmente a 
ellas, independientemente de su valor de energía y de la altura de las barreras 
de potencial. La doble barrera de potencial también es transparente para 
aquellos electrones con q;, = re y qa, = 123 , es decir cuya energía y ángulo 
de incidencia cumplen la doble condición de la ecuación Eq[2.40] a la vez. Ya 
que si no, podría suceder que una barrera fuera transparente y la otra no. 





2 
(7) = (E-U)? — (Esing), ¡=1,2 (2.40) 

L; 

De forma semejante al caso de una sola barrera de potencial, consideramos 
que los electrones inciden con un ángulo inferior al menor de los dos ángu- 
los límite del sistema, asociados a cada barrera de potencial en la ecuación 
Eq[2.23]. Si lo hacen con un ángulo superior a alguno de ellos, o a ambos, 
el análisis sigue siendo válido, y el único cambio que se ha de realizar en 
las fórmulas anteriores es cambiar q;, por 1q;,, lo que transforma las ondas 
planas e*iz" en ondas evanescentes eFtis". En el interior de una o de las 
dos barreras de potencial, el transporte de electrones no contiene términos 
oscilatorios sino exponenciales reales. En cuanto a los coeficientes, se ha de 
cambiar cos(L;q;,) por cosh(£L;q;,) y sin(L;q;,) por ¿sinh(L;q;,.). 








En función de la energía del electrón incidente, del ángulo de incidencia 
y de los valores característicos de las dos barreras de potencial, así como 
de la separación entre ellas, podemos obtener una expresión explícita del 
coeficiente de transmisión para la doble barrera de potencial: 
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a = 00 (yy) os (y) 
= (20,0, tan? ósin? (+) + (E — UL) (E — Us) — E* sin 6?) 
sin (qe) sin (322) (2.42) 


€1 €2 


























: Lie : Loae2 
; 9D sin (4 ) sin (8 ) 
Za = UyUztan? ósin (5 =) du da 


UF €1 €2 











Es Lie 
U, — Ecos? $ (E (a) 
cos 











cos Q hup él 
: Lo€2 
E 2 L sin ( > ) 
Za a COS ( 2) dE (2.43) 
cos hup €2 


siendo 





e = Ecos, € = YE —U¡?-—E?sin?g  ¡=1,2 (2.44) 


En la figura Fig. 2.10 se representa el coeficiente de transmisión frente 
al ángulo de incidencia para dos barreras de potencial de alturas U¡ = 0.2 
eV y Uz = 0.4 eV, anchuras L; = 50 nm y Lz = 60 nm, separadas una 
distancia D = 0, 20, 40, 60 nm, siendo la energía incidente 0.1 eV en (a) 
y 0.3 eV en (b). Los parámetros utilizados son los típicamente considerados 
en las referencias [59, 67,74]. En dicha figura Fig. 2.10 se aprecia que, para 
una energía fija, la envolvente del coeficiente de transmisión tiende a decrecer 
más o menos abruptamente a medida que aumenta el ángulo de incidencia. 
La existencia del ángulo límite en la ecuación Eq[(2.23] que dependía de la 
energía y de la altura de la barrera de potencial es evidente al comparar 
ambas figuras. 


En la figura Fig. 2.11 se representa el coeficiente de transmisión fren- 
te a la energía de la partícula incidente para dos barreras de potencial de 
parámetros U, = 0.2 eV, Uz = 0.4 eV, L¡ = 50 nm, L2 = 60 nm, y dis- 
tancias de separación D = 0 , 20, 40 y 60 nm, siendo el ángulo de inci- 
dencia $ = 7/36 en (a), 7/12 en (b) y 7/6 en (c). En (b) se aprecia que, 
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Figura 2.10: Coeficiente de transmisión frente al ángulo de incidencia para dos 
barreras de potencial de parámetros U, = 0.2 eV, Us = 0.4 eV, L, = 50 mm, 
L3 = 60 nm, para diversos valores de D (nm), siendo la energía incidente 0.1 
eV en (a) (0.3 eV en (b)). En (b) la envolvente del coeficiente de transmisión 
cae abruptamente asociado a la existencia de un ángulo límite. 
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Figura 2.11: Coeficiente de transmisión frente a la energía de la partícula 
incidente para dos barreras de potencial de parámetros U, = 0.2 eV, U»z 
= 0.4 eV, L¡ = 50 nm, L3 = 60 nm, para diversos valores de D (nm). El 
ángulo de incidencia $ = 7/36 en (a), 7/12 en (b), $ = 7/6 en (c). Para 
ángulos de incidencia f suficientemente grandes aparecen dos pseudogaps, 
que al aumentar H pueden llegar a colapsar en un único pseudogap. 
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para un ángulo de incidencia fijo, el coeficiente de transmisión cae práctica- 
mente a cero en el intervalo energético [(1 + sin 49)7*U,, (1 — sing)"*U1], a 
continuación se recupera, volviendo a tomar valores próximos a la unidad, 
y vuelve a caer de nuevo a un valor casi nulo para el intervalo energético 
[1 + sin 6) 10», (1 — sin p)*U»]. En estos intervalos, el ángulo de incidencia 
es superior al angulo límite y la transmisión de electrones se hace por me- 
dio de ondas evanescentes en las zonas asociadas a las barreras de potencial 
donde suceda esto. Sin embargo, en (a) el ángulo de incidencia no es suficien- 
temente grande, como para que el coeficiente de transmisión caiga a cero en 
el primer mínimo, aunque sí en el segundo. Tanto en (a) como en (b), para 
6 >0, si (1—sing)*U, < (1+sin 9) *U, se pueden distinguir en el espectro 
de energías dos regiones claramente diferenciadas en las cuales el coeficiente 
de transmisión es mínimo. Si (1 —sing)*U, > (1 +sin (/)7*Uz ambas regio- 
nes donde el coeficiente de transmisión es mínimo se solapan y sólo puede 
distinguirse una única zona de valor mínimo. Esto es lo que sucede en (c). 
En este caso, a medida que aumenta el ángulo de incidencia f, aumentará la 
anchura del pseudogap. En los dos últimos casos subrayamos la existencia 
de pseudogaps efectivos de energía, dentro de los cuales, la probabilidad de 
transmisión del electrón cae prácticamente a cero. 


En los experimentos la energía potencial de las barreras U¡ y Uz2 son 
sintonizables por medio de los potenciales externos aplicados a las puertas 
superiores V, y V2. Puesto que su ajuste resulta sencillo, interesa ver los 
cambios que se producen en la transmisión para distintos valores de dichos 
potenciales. En la figura Fig. 2.12 represento algunas líneas de nivel del coefi- 
ciente de transmisión como función de los potenciales U; y Uz, para diferentes 
valores de energía E y ángulo de incidencia /. La figura Fig. 2.12 nos da las 
líneas de nivel para T = 0.1 (amarillo), 0.3 (naranja), 0.5 (rojo), 0.7 (verde) 
y 0.9 (azul), siendo L, = 50 nm, L2 = 60 nm, D = 20nm, $ = 7/12, para E 
= 0.05 (panel superior izq.) y 0.1 eV (panel superior dcho.). El panel inferior 
izq. y dcho. nos da las líneas de nivel con los mismos parámetros que las 
figuras anteriores, excepto y = 7/6, para E = 0.3 y 0.6 eV respectivamen- 
te. Al comparar estas cuatro figuras dos a dos, vemos que a un ángulo fijo 
cuando duplicamos la energía del portador incidente, aparece una especie de 
pasillo de pseudogap, centrado en U, = Uz = E, cuya anchura coincide con 
el valor de E que estemos utilizando. Es claramente visible en las figuras del 
panel superior, sin embargo en las figuras del panel inferior se aprecia menos, 
ya que en la figura del panel inferio derecho el pasillo es tan ancho que no 
aparece reflejado en su totalidad en la gráfica. 
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Figura 2.12: Líneas de nivel para para T = 0.1 (amarillo), 0.3 (naranja), 
0.5 (rojo), 0.7 (verde) y 0.9 (azul), siendo L, = 50 nm, La = 60 nm, D = 
20nm. Parámetros : E y 6 = 0.05 eV y 7/12 (panel superior izq.), 0.1 eV 
y 7/12 (panel superior dcho.), 0.3 eV y 7/6 (panel inferior izq.) y 0.6 eV 
y 7/6 (panel inferior dcho.). Aparece un pasillo de pseudogap centrado en 
U, = Uz = E, cuya anchura coincide con la energía que estemos utilizando. 
Las figuras de la derecha tienen una energía doble a las de la izquierda, por 
lo que la anchura del pasillo en tales casos se duplica. 
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2.4. Resonancias para una doble barrera de 
potencial 


Para dos barreras simétricas solamente existe un intervalo energético para 
el pseudogap, que coincide en valor con el que habíamos obtenido para una 
única barrera de potencial. Esto es debido a que únicamente dependía de la 
altura de la barrera de potencial y del ángulo de incidencia, que coinciden en 
ambos casos: 

[(1 + sin | 6 [)7*00, (1 — sin | y 1)7*00] (2.45) 


La aplicación de la transmisión perfecta T = 1 a la ecuación Eq[2.41] 
permite obtener una expresión cerrada de las energías que se transmiten 
integramente a través de las dos barreras simétricas. Obtenemos transmisión 
sin reflexión para el caso de barreras idénticas cuando se verifica: 











, Sin SEE, 0) y, Cos? $ — lsin[¿2F(E, 6)] DE 
4tan* o FE OY ( 200 FE 0) A cos q] 
LU, 2 
= Sia cos (] cos[, PUE, 6)1) =0 (2.46) 
siendo 
E E? 
F(E,Q) = ya = no UE sin? $ (2.47) 


que contiene el caso asociado al túnel de Klein. En el caso de una única 
barrera las resonancias se encuentran fuera del pseudogap. Sin embargo en 
el caso de dos barreras simétricas (y también para barreras desiguales), apa- 
recen además, superpuestas, resonancias sumamente estrechas (tipo delta) y 
equiespaciadas en el pseudogap, como se aprecia en la figura Fig. 2.13, donde 
se representa el coeficiente de transmisión para una barrera de potencial en 
rojo (Uy = 0.2 eV, L= 100 nm y 6 = 7/30), y para dos barreras de potencial 
con los mismos parámetros (U¡ = Uz = 0.2 eV, L, = La = D = 100 nm y 
6 = 1/30) en azul. Imponiendo que el coeficiente de transmisión sea unitario 
para dos barreras idénticas, obtendremos la expresión Eq[2.48], a partir de 
la cual obtendremos las resonancias: 





7 cos? p — 1 DE A 
2 tan cosp| = 
cos Y hop tanh ES ve sin? - 1)? 











(2.48) 
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Figura 2.13: En (a) T frente a la energía para una y dos barreras de potencial 
simétricas con Uy = 0.2 eV, L = D = 100 nm y 6 = 7/30. En (b), ampliación 
de la zona central de (a). Para dos barreras de potencial simétricas aparecen 
resonancias tipo delta equiespaciadas en la zona donde teníamos un pseudo- 
gap en el caso de una única barrera de potencial. 


Dichas resonancias, en nuestro ejemplo de la figura Fig. 2.13 (b) corres- 
ponden a los valores de energía E = 0.1836, 0.1987 y 0.2143 eV. 


Un aumento en el ángulo de incidencia H provoca un aumento en el inter- 
valo energético asociado al pseudogap en la expresión Eq[2.45]. Observamos 
que el aumento de £ (asociado a la anchura de las barreras de potencial), 
pero sobretodo de D (asociado a la distancia de separación entre las barre- 
ras simétricas), provoca un aumento significativo en el número de picos de 
resonancia dentro del pseudogap. Ello es debido a que la parte izquierda de 
la expresión Eq[2.48] depende de D implícitamente a través de una función 
tangente, mientras que la parte derecha de la expresión Eq[2.48] depende de 
L implícitamente a través de una función tangente hiperbólica que se en- 
cuentra en el denominador. 


Si nuestra doble barrera de potencial deja de ser simétrica, entonces la 
ecuación Eq[2.48] deja de ser aplicable para detectar picos en el pseudogap. 
A medida que se modifica la anchura, pero sobre todo la altura, de una de 
las barreras de potencial con respecto al caso de la doble barrera simétrica de 
potencial, los picos que aparecían en el pseudogap disminuirán en amplitud 
y se irán degradando progresivamente. 


Como conclusión, diremos que es interesante resaltar el hecho de que 
sobre el pseudogap puedan superponerse resonancias sumamente estrechas 


(tipo delta) a las energías definidas por la ecuación Eq[2.48)]. Creemos que 
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estas resonancias en energía, dentro de una zona donde básicamente no existe 
transmisión, puedieran ser utilizadas para la sintonización puntual a ciertos 
valores energéticos de la transmisión, con una tecnología en la cual se pu- 
diera controlar la dirección de transmisión del haz electrónico. Es decir, el 
dispositivo resultante se comportaría como filtro pasa-banda muy estrecho o 
filtro sintonizado. 


2.5. Resumen y conclusiones 


A lo largo de este capítulo hemos analizado los estados en una lámina 
de grafeno sometida a una o dos barreras de potencial. Para ello hay que 
tener en cuenta básicamente dos consideraciones: la energía es una constante 
del movimiento, y también que al considerar un potencial unidimensional 
U = U(x), la componente k, del vector de ondas mantendrá su valor en 
todas las zonas. Además, la relación entre el ángulo incidente / y el ángulo 
transmistido O en una barrera de potencial, depende de la relación entre la 
energía del portador incidente y la altura de la barrera de potencial E /Up, de 
tal forma que si es inferior a 0.5, el ángulo transmitido se acerca a la nomal; 
si es superior a 0.5 se aleja de la normal, en este contexto aparece el concepto 
de ángulo límite, que más adelante será clave; y si coincide con 0.5 el ángulo 
transmitido y el incidente coinciden. 


Obtenemos una expresión analítica exacta para el coeficiente de trans- 
misión de una lámina de grafeno sujeta a una barrera de potencial eléctri- 
co. Debido a que el coeficiente de transmisión es unitario para inciden- 
cia perpendicular (paradoja de Klein) o bien para áquellos parámetros que 
verifiquen la condición de la fórmula Eq|2.17], no puede existir confina- 
miento de los portadores de carga mediante barreras de potencial eléctri- 
co en el grafeno prístino con las características que hemos adoptado. La 
aparición del ángulo límite se observa claramente en la figura Fig. 2.5, ya 
que el coeficiente de transmisión decae abruptamente. Más interesante aún, 
es el hecho de que al representar el coeficiente de transmisión frente a la 
energía para valores dados de L - Uy y Q, observamos simetría en torno al 
eje E = Upcos”? 4, donde toma su valor mínimo T = cosh”*[LU¿/(hvp)]. 
Además para ángulos de incidencia suficientemente elevados en el intervalo 
energético (1+sin |p|)+Uy < E < (1—sin |6])*U0, el coeficiente de transmi- 
sión es prácticamente nulo apareciendo así un pseudogap de energía de valor 
200 sin |p| cos”? H en torno al eje de simetría. 


También hemos obtenido una expresión analítica exacta en el caso de 
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dos barreras de potencial. Observamos que el coeficiente de transmisión es 
independiente del orden en que se coloquen éstas. Aparece el fenómeno de 
acoplamiento entre ellas, siempre que la distancia de separación D no sea 
D > Li, Lo. También para dos barreras de potencial, al analizar el coeficien- 
te de transmisión - dependiendo del ángulo de incidencia - puede aparecer un 
pseudogap de energía asociado a cada barrera de potencial, que colapse en 
un único peudogap en energía si se verifica (1 — sin p)7*U, > (1 +sin $) *U». 


Por último, hemos considerado el caso de una doble barrera simétrica 
de potencial. Al analizar el coeficiente de transmisión observamos que tiene 
lugar la aparición de resonancias muy estrechas (tipo delta) dentro del pseu- 
dogap en energía, lo cual es sumamente interesante. Los valores energéticos a 
los cuales aparecen dichas resonancias vienen dados por la ecuación Eq[2.48]. 
El número de dichas resonancias extremadamente estrechas aumentará con- 
siderablemente con el aumento de la distancia de separación D entre barreras 
de potencial. Creemos que estas resonancias en energía, dentro de una zona 
donde básicamente no existe transmisión, puedieran llegar a ser utilizadas 
para la sintonización a ciertos valores discretos de energía en la transmisión, 
siempre que contasemos con una tecnología en la cual se pudiera controlar 
la dirección de transmisión del haz electrónico. Es decir, el dispositivo resul- 
tante se comportaría como filtro sintonizado. Subrayar que la mayor parte 
de los resultados de este capítulo aparecen en nuestro artículo [71]. 
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Capítulo 3 


Influencia del ángulo de 
incidencia en la conductancia 


3.1. Introducción 


Una vez calculado el coeficiente de transmisión para una lámina de gra- 
feno a través de una y dos barreras de potencial, vamos ahora a determinar 
la conductancia con objeto de comparar nuestros resultados con los expe- 
rimentos. Así, conseguiremos enlazar el cálculo teórico con una magnitud 
física experimental fácilmente medible. A partir de la densidad de corrien- 
te determinaremos una expresión general analítica con que caracterizar la 
conductancia efectiva, es decir la conductancia por unidad de anchura, del 
dispositivo. Construimos en este capítulo un modelo de funcionamiento para 
un transistor de efecto de campo de grafeno (GFET) y analizamos su con- 
ductancia por aplicación de la fórmula anterior, haciendo uso del coeficiente 
de transmisión determinado en el capítulo previo. Un aspecto relevante en 
el funcionamiento del GFET es el carácter direccional o no del ataque de 
los portadores sobre la barrera. ¿Tienen todos los electrones la misma di- 
rección de impacto o se distribuyen en diferentes direcciones con ángulos de 
incidencia diferentes? Analizamos este fenómeno asignando pesos estadísti- 
cos diferentes a las distintas direcciones de incidencia y encontramos en la 
conductancia del GFET una casuística variada que discutimos gráficamente. 


3.2. Conductancia efectiva 
El propósito de este apartado es aplicar los resultados obtenidos en el 
capítulo anterior que determinaban analíticamente el coeficiente de trans- 


misión a través de una barrera de potencial, en la descripción del funcio- 
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namiento de algunos dispositivos electrónicos, principalmente de un GFET, 
que debe definirse mediante su ecuación característica corriente-voltaje o, 
equivalentemente, mediante su conductancia. Para ello, debemos estudiar el 
transporte de portadores a través del dispositivo cuantificando la corriente 
eléctrica transportada por cada uno de los electrones que atraviesa el dispo- 
sitivo ocupando alguno de los estados electrónicos de alguno de los canales 
de conducción disponibles. Utilizaremos para ello el formalismo de Landauer- 
Búttiker [75,76] que conecta el coeficiente de transmisión con la densidad de 
corriente y la conductancia. 


Consideramos un GFET cuyo canal activo de conducción es una lámina 
de grafeno de dimensiones nanoscópicas, sometido a diferentes puertas de 
potencial que pueden establecer una barrera de potencial como la estudiada 
en [71] y seleccionar una cierta densidad de portadores. Dicho GFET se pola- 
riza por medio de una tensión externa drenador-fuente Vps que proporciona 
una corriente convencional en el sentido del drenador al surtidor, correspon- 
diendo a un transporte de electrones en sentido contrario. 


De acuerdo con Landauer, el papel de los contactos (S) y (D) es el de sen- 
dos reservorios que proporcionan un colectivo de electrones, con diferentes 
valores de energía y momento, en equilibrio termodinámico a un potencial 
químico uz en la fuente y ug en el drenador [69], cuya diferencia se esta- 
blece precisamente mediante el potencial de polarización 4, — Hr = eVps. 
Estos electrones compiten por llenar los estados cuánticos disponibles en la 
lámina de grafeno y generan sendas corrientes eléctricas cuya diferencia es la 
corriente total. 


La existencia de una polarización no nula en el dispositivo rompe el equili- 
brio establecido por los reservorios, sin embargo, se puede demostrar [69,77], 
que cuando tanto la polarización como la temperatura son suficientemente 
bajas, se pueden ignorar estos efectos y tratar el colectivo de electrones co- 
mo un conjunto en equilibrio a un potencial químico intermedio. En estas 
condiciones, cada electrón de energía E y momento k = (k.,,k,) = (k, 6), 
procedente de uno de los reservorios (por ejemplo del surtidor) que entra en 
el estado Y de uno de los canales de conducción de la lámina de grafeno de 
longitud L y anchura W, transporta una corriente eléctrica proporcional al 
flujo de probabilidad dado por 


Ji =—eop»* o, V (3.1) 


de modo que a la salida del dispositivo, la corriente del surtidor resulta 
proporcional a su coeficiente de transmisión. En efecto, 
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UF Ka 


se TE, Hp) se ET(E, $)cosp (3.2) 


Je evupV*o k] LW 


donde s = sign(E), y j? tiene dimensiones de intensidad de corriente por 
unidad de longitud. 


La corriente total del emisor es debida a todos los electrones con to- 
dos los posibles valores de k, y k,, que se introducen en el sistema y se 
propagan a través de los orbitales disponibles. Su número, con momentos 
comprendidos entre (k,, k,) y (kz + dk, ky + dk,), depende del número de 
estados disponibles y de su grado de ocupación. Se conoce que cada momen- 
to permitido ocupa una superficie en el espacio de momentos de un grafeno 
bidimensional de FE Xx > y por ello su número en un rectángulo dk.dk, es 
2x2 x dksdky (E x 57) = (LW/1?) dk dky si se considera la duplicidad de 


spin y la de valle [78]. 


En un sistema en equilibrio, la distribución de estos electrones está dada 
mediante la función de Fermi f,(£, yu) [79] 


1,1) = (14 0858)" (3.3) 


de modo que el número de electrones en este intervalo, en equilibrio a un 
potencial de Fermi y, procedentes del surtidor con un factor de ocupación 


JE, pL), es 


Melee =2 Xx 2 X ta dy 





LW 
fAE, 11) = Es kE, pudaicdiy (34) 
y cambiando al sistema de coordenadas polares dk.,dk, = |k| d[k|dó, e 
introduciendo la relación de dispersión E = shup|k|, tendremos 


Últ = ALO tE F(E, ur) (340) 


La densidad de corriente a que dan los estos electrones es entonces 
de 
h2vp 
así, la densidad de corriente total procedente del surtidor debida a todos los 

electrones con todas las posibles energías y ángulos de incidencia es 


Ide e 
dE hau J_ 


djs == —Sol E, H)T(E, 6) EdE cos ódó (3.6) 





cosado | T(E-S)LAE MEE (37) 
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Un cálculo análogo permite calcular la corriente procedente del drenador. 
Los portadores procedentes del drenador y surtidor compiten entre sí para 
transportar corriente a través del canal. Por ello, la corriente neta que fluye 
a través del canal es 


de [2 ji 
jo== 2 | cosódo | (TUE. OL(E qu) -Ta(E, 9) kE 10) EdE 


(3.8) 

De acuerdo con las propiedades de simetría de los coeficiente de transmi- 

sión T¿(E,9) y Tr(E, 6) son idénticos, y por tanto obtenemos la expresión 
general para la densidad de corriente que circula por el dispositivo: 


0 de 'ñ 
Ja = AN 


cosódo | T(EO (ELE qu) Jl Est) EE (39) 


En algunos casos y bajo ciertas condiciones, la expresión general anterior 
Eq[3.9] puede ser simplificada: 


= Si la polarización existente a bajas temperaturas es muy grande, en- 
tonces la contribución debida al reservorio drenador es despreciable ya 
que su potencial químico es despreciable frente al de la fuente. En dicho 
caso aproximaremos la corriente total a áquella debida únicamente al 
reservorio fuente 





4 3 
a ] cos o dó f re Of(E, uE dE. (3.10) 
h?up as 


= A baja temperatura (kgT < Ep — y), cuando el colectivo de electrones 
está altamente degenerado, la función de Fermi puede aproximarse por 
la función escalón f.(E, pu) — f(E, pa =0(E, ur) — 0(É, ), dando 


de 2 
le = = +5 9)EdE. ; 
q e [esoo fra d (3.11) 


= Si la polarización es muy pequeña, podemos aproximar la distribución 
de Fermi por su aproximación de Taylor en torno al valor medio de 
los potenciales químicos u = ES y considerar que los potenciales 
químicos se establecen por medio de la polarización uu, — pg = e (Ve — 
V.) = e Vos, es decir 


Of. o 


f(E, pr) — SE, ur) = “du Vos = OE 


ia (3.12) 
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y por tanto, 
2 E OLE, 
cosodo [ T(E, $) ( — HARD paz 
(3.13) 


= A temperaturas ultrabajas (a cero kelvin donde y = Ep») y polarización 
muy baja, podemos aproximar 


Ofo(L, yu) 
OE 


así a partir de las expresiones anteriores Eq[3.13] y Eq[3.14] nos que- 
dará, teniendo en cuenta la simetría del coeficiente de transmisión 


TE, 6) 3 T(E, 0), 


=0(E—u) =0(E-— Ep) (3.14) 


de? 7 
de =- Vos Er o T(Ep, 6) cosó do 
UF er 


NS ula 


2 
Ñ e Vos Er | T(Er, 6) cos p do (3.15) 
h Up 0 


En todas las fórmulas anteriores, el signo menos describe la dirección correcta 
de la corriente eléctrica transportada por electrones y huecos, en el sentido 
del potencial eléctrico decreciente. 

En ciertas medidas experimentales, el interés se enfoca no tanto en la 
densidad de corriente j.,, como en la conductancia de la muestra G. A partir 
de la ley de Ohm y la intensidad de corriente tenemos 


W 
Vos G = [ = / Jady = Ja W (3.16) 
0 


y podemos definir la conductancia efectiva o conductancia por unidad de 


longitud (anchura) Geoff 
G da 
Gerr= == 3.17 
EE (3.17) 
Aplicando la ecuación anterior Eq[3.17] a nuestra expresión general Eq[3.9] 
obtenemos 





Car f cosóds [TI (LLE 1) — JE 10) Ez 


h?vpVps 








ELA 
2 


(3.18) 
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En particular, para las condiciones experimentales expresadas en Eq[3.15] 
tenemos 





de? 7 
di zon | T(Er, 6) cos $ de (3.19) 


h 


e 
2 


La conductancia efectiva del dispositivo depende del potencial aplicado a la 
puerta superior (top gate) por medio del coeficiente de transmisión, y del 
potencial aplicado a la puerta inferior (back gate), asociado a la energía de 
Fermi Ep, pero no del potencial de polarización. Ello es debido a que la 
cuestión es que no importa la forma del potencial (de polarización) porque 
éste finalmente afecta al transporte a través del factor de ocupación [77]. 

Diferentes dispositivos con diferentes perfiles de potencial darán lugar a 
diferentes coeficientes de transmisión, que deberán ser determinados en cada 
caso, pero la fórmula Eq[(3.19] seguirá siendo válida mientras las condiciones 
experimentales bajo las cuales se obtuvo se mantengan. En lo sucesivo con- 
sideraremos las condiciones de trabajo que permiten utilizar la conductancia 
efectiva Eq[3.19], esencialente T = 0 K y Vps = 0. 


3.3. Modelo para el GFET 


La figura Fig. 3.1 muestra el esquema del dispositivo considerado: una 
lámina de grafeno rectangular (azul) encapsulada entre capas de h-BN [80], 
con contactos metálicos de Au/Cr (naranja) aunque en el dibujo únicamente 
hemos puesto Au por cuestión de espacio, que constituyen la fuente (source), 
el drenador (drain) y la puerta superior (top gate), y una puerta inferior 
(back gate) constituida por silicio dopado (amarillo), que se colocan bajo 
una capa de dieléctrico (SiO [81]). 

Las funciones que realizan las puertas inferior y superior en el esquema 
de la figura son diferentes. La puerta inferior en la figura Fig. 3.1 contro- 
la la densidad de portadores de carga y con ella la energía de Fermi de los 
electrones de transporte en el grafeno, y la puerta superior junto con la in- 
ferior modulan la corriente que fluye a través del dispositivo desde fuente 
a drenador, pasando a través del grafeno donde hemos aplicado potencial 
mediante la diferencia de tensión entre la puerta superior y la inferior Vr - 
Vea. Supondremos que la lámina de grafeno es suficientemente ancha para 
no ser considerada como una cinta, es decir la anchura debe ser mayor que la 
longitud de scattering W > Lscat1, la cual típicamente es del orden de 1 um, 
y por tanto no existirán efectos de borde. Además a temperatura próxima a 
O K, la densidad de portadores de carga n tiene una dependencia lineal con 
el voltaje aplicado a la puerta inferior (n x Va) [82] - como detallaremos 
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a continuación -, y una dependencia cuadrática con la energía de Fermi [7] 
(n x Ez.) como queda reflejado en el Apéndice A. 





Figura 3.1: Esquema del GFET. La lámina de grafeno representada en azul, 
las láminas de h-BN en gris, los contactos metálicos en naranja, el dieléctrico 
en rosa y el silicio dopado en amarillo. 


La energía de Fermi del grafeno puede modificarse aplicando voltaje a 
la back gate. Idealmente el pico de Dirac Vp, también denominado punto 
de neutralidad de carga debería suceder para un voltaje nulo aplicado a la 
back gate. Sin embargo, experimentalmente no es así, lo cual es debido a la 
aparición de dopado, dando lugar a una discrepancia respecto al caso ideal 
en la posición del pico de Dirac, que puede llegar a ser del orden de decenas 
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de voltios. 


Nuestro dispositivo se comporta como un condensador de placas plano- 
paralelas, constituyendo las partes metálicas, el grafeno y el silicio dopado, 
separadas una distancia d (típicamente de 300 nm), por medio de un dieléctri- 
co (con constante dieléctrica e) que en nuestro caso es SiO». Así, al aplicar 
un voltaje Vga dará lugar a una carga eléctrica, obteniendo una relación de 
proporcionalidad entre la densidad de electrones n (o de huecos p) y el voltaje 
aplicado a la back gate Vga. 


6 (Vaa — Vo) 
2 ed 

Por ejemplo en [83], « = 7.2 -10%% em”? V"! para un dieléctrico SiOz 
de 300 nm de grosor. En la expresión Eq[3.20] se ha introducido Vp porque 
se toma como origen para el potencial, ya que constituye un eje de simetría 
para la densidad de portadores. Cuando el voltaje aplicado a la back gate 
es próximo al punto de neutralidad de carga, la densidad de portadores es 
menor que cuando nos alejamos, ya que idealmente la densidad de portadores 
es nula en el punto de Dirac. Lo anterior nos permite modificar la naturale- 
za de los portadores ( tipo electrón o hueco) al trabajar con tensiones Vga 
superiores o inferiores al pico de Dirac Vp. 


= a(Vaa — Vo) (3.20) 


n—=p 


La resistividad p está asociada a dos componentes, una de largo alcance 
pL y otra de corto alcance ps: 


p=PL+Ps = == che (3.21) 
in—p|ep 

siendo u la movilidad de los electrones si Vga > Vp, o bien la movilidad 
de huecos en caso contrario. La resistividad aunque parece tender a infinito 
al aproximarnos a Vp, sin embargo alcanza idealmente un valor máximo de 
h/4e? = 6.45 kQ en dicho punto, el cual está asociado a una conductividad 
cuántica mínima. La conductividad o = 1/p tiene la característica inversa a 
la resistividad, y nunca alcanzará el valor nulo. Además, la pendiente de las 
curvas 0 = f(V,) lejos del punto de Dirac, caracterizan a la movilidad. 


3.4. Conductancia efectiva del GFET 


En el GFET de la figura Fig. 3.1 los electrones circulan por el canal como 
partículas ultrarrelativistas sin masa sujetas al formalismo de Dirac [84]. 
Atraviesan una barrera de potencial U, establecida por la puerta superior y 
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Figura 3.2: Energía de portadores versus posición en la lámina de grafeno 
(a). Angulos de incidencia y transmisión de portadores sobre la barrera de 
potencial (b). 


se transmiten por el dispositivo de acuerdo con el coeficiente de transmisión 
T(E, $) obtenido en el capítulo anterior [71]. El coeficiente de transmisión a 
través de la barrera de la figura Fig. 3.2, reulta ser: 
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sin? (EVE — Up)? — E? sin” 6) 
(E —U0)? — E? sin? ¿ 





T(E,U0,D,6) = | 1+Ujtan?o 


(3.22) 
depende de características propias del portador de carga como son la 
energía E, y el ángulo de incidencia sobre la barrera de potencial /, y de 
parámetros propios de la barrera de potencial como son su altura Uy y su an- 
chura D. No depende en absoluto de las dimensiones de la lámina de grafeno 
W y L (siempre que sus dimensiones permitan prescindir de los efectos de 
borde). El coeficiente de transmisión es simétrico respecto al ángulo de inci- 
dencia H y toma el valor 1 para incidencia normal, con independencia de los 
valores que puedan tomar el resto de los parámetros. Este es el caso de la pa- 
radoja de Klein para las partículas de Dirac [57,62,85], tal y como habíamos 
visto. Haciendo uso del coeficiente de transmisión en la ecuación Eq[3.22] 
construimos la conductancia efectiva de la muestra a baja temperatura y 
baja polarización en la expresión Eq[3.19]. Así, a temperaturas próximas a 
0 K, la conductancia efectiva de este GFET en unidades de e?/h (que es la 
conductancia básica cuantizada) es: 





2 +7 /2 
Geri = Eo EF | T(Ep, Uo, D, $) cos $ dó (3.23) 
Up T —1/2 
y explícitamente 
4|Ep| 
Gerr[ Ep, Un, D) = , 
pr(Er, Uo, D) E 


-1 





7/2 


D Duaril 
sor (Er — Up) — Ep sin 6) 


sin? ( 
1 + Ujtan?g 
(Ep — Up)? — E? sin? o 





-cos $ do (3.24) 


La conductancia depende del potencial aplicado a la puerta superior (a 
través del coeficiente de transmisión), del potencial aplicado a la puerta in- 
ferior (a través de Ep), y de la anchura de la lámina de grafeno Eql3.16], 
pero no de su longitud Eq[3.17]. El máximo valor de la conductancia efec- 
tiva Eq[3.23] se alcanza para el valor máximo del coeficiente de transmisión 
T(Ep,Uo, D,p) = 1, y por tanto 
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Figura 3.3: Conductancia efectiva frente a la energía de Fermi para una barre- 
ra de potencial de altura Uy = 200 meV y anchura D = 50(rojo), 100(verde) 
and 150(azul) nm. La conductancia tiene un máximo local en Ep =0.5-Up, un 
mínimo local en Ez = Up, y a continuación crece linealmente para Ep > U6. 


La figura Fig. 3.3 muestra la conductancia frente a la energía de Fer- 
mi para algunos valores de parámetros, especificados en el pie de la figura. 
Usamos como unidades para la energía me V, y para las longitudes nm (ade- 
cuadas para una lámina nanométrica de grafeno), tomamos hup = 658.2 
meV-nm y expresamos la conductancia en unidades proporcionales a e?/h. 
Vemos que la conductancia efectiva tiene un máximo local cuando la energía 
de Fermi alcanza el valor Ep = 0.5-Up. A energías superiores, la cantidad 
(Es — Up)? — Ed sin? d cambia el signo, pasando de positiva a negativa. Por 
ello, en la expresión Eq[3.24] el seno de esta cantidad se convierte en un seno 
hiperbólico, y la conductancia comienza a disminuir. Por tanto, la resistencia 
aumenta. Esta situación persiste hasta que se alcanza un mínimo local en 
Ep 2 Uy. Cuando Ep es mayor que Uy, debido a que el coeficiente de trans- 
misión tiende a ser unitario, la conductancia efectiva crece proporcionalmente 
a la energía de Fermi, de forma semejante a lo que habíamos obtenido en la 
ecuación Eq[3.25]. Observamos oscilaciones en la conductancia a energías de 
Fermi que verifican la condición Ep < Up. En particular, para Ep <0.5-0b, 
existen varias oscilaciones, mientras que para valores de la energía de Fermi 
que verifican 0.5-U4 < Ep < Uy solamente existe una oscilación. Para valores 


61 


Conductancia 





superiores, es decir, si Ep >> Up, la conductancia no depende del valor de 
la anchura de la barrera y varía linealmente con la energía de Fermi. 


En la figura Fig. 3.4, donde se representa la conductancia efectiva frente 
al potencial de barrera para una energía de Fermi Ep = 100 meV y anchura 
D = 50(rojo), 100(verde) y 150(azul) nm, se observa que a medida que Up 
aumenta, para valores de Uy que satisfacen Uy < Ep, la conductancia efectiva 
decrece linealmente y toma el mismo valor independientemende de la anchura 
de la barrera de potencial. Cuando Uy coincide con el valor de la energía de 
Fermi, aparece un mínimo local cuyo valor depende de la anchura de la 
barrera D. Para Uy > Ep, la conductancia aumenta y a la vez aparecen 
oscilaciones cuya amplitud se incrementa a medida que disminuye la anchura 
de la barrera. 
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Figura 3.4: Conductancia efectiva frente a la altura de la barrera de poten- 
cial para una energía de Fermi Ep = 100 meV y D = 50(rojo), 100(verde) 
y 150(azul) nm. Las tres curvas coinciden para Uy < Ep, y presentan un 
mínimo en Uy = Ep. La amplitud de las oscilaciones depende de la anchura 
de la barrera D. 


La dependencia de la conductancia efectiva con la anchura de la barrera 
D se muestra en la figura Fig. 3.5. La conductancia disminuye hasta que se 
estabiliza, mediante oscilaciones, en torno a un valor constante. La amplitud 
de dichas oscilaciones tiende a disminuir a medida que la energía de Fermi 
aumenta. 
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Figura 3.5: Conductancia efectiva frente a la anchura de la barrera, para Up 
= 100 meV y energías de Fermi Ej = 25(rojo), 50(verde), 75(azul), 100(ama- 
rillo) y 150(fucsia) meV. La conductancia disminuye hasta que se estabiliza a 
través de oscilaciones, cuya amplitud disminuye a medida que Ep aumenta. 


3.5. Conductancia efectiva ponderada 


Hasta ahora hemos considerado que, durante el transporte de corriente, 
la probabilidad de impacto de los portadores sobre la barrera de potencial 
es la misma para cada ángulo de incidencia. Otros escenarios deberían ser 
considerados, en los cuales, debido a las inhomogeneidades del material o al 
diseño del dispositivo, los portadores fluyeran en su mayoría en torno a una 
cierta dirección privilegiada dada fp o, con mayor generalidad, tuviéramos 
un conjunto de electrones cuyas direcciones de entrada estuvieran descritas 
por una cierta distribución de probabilidad P(p) representando el peso de 
una dirección dada en el colectivo de electrones. A continuación considerare- 
mos algunos de estos diferentes escenarios por medio de una distribución de 
probabilidad dada [86], y estudiaremos el efecto de la misma sobre la conduc- 
tancia efectiva del GFET. Naturalmente, esta función de probabilidad debe 
estar normalizada, y satisfacer por tanto la condición 


+7 /2 
. P(óldó= 1 (3.26) 


7/2 
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3.5.1. Distribución gaussiana 


Consideramos primero el caso de una función de distribución de proba- 
bilidad gaussiana para el ángulo de incidencia. Esta suposición representa la 
idea de que el transporte es de naturaleza balística, y que existe una cierta 
inhomogeneidad que provoca procesos de dispersión dentro de la muestra, 
los cuales desvían a los portadores de su trayectoria balística, llevándolos a 
impactar contra la barrera de potencial con diversos ángulos de incidencia. 
Ya que no todos los portadores sufrirán los mismos procesos de dispersión, 
podemos asumir que la mayoría de los portadores impactarán en la barre- 
ra de potencial con un ángulo similar a dp, aunque puede haber una cierta 
dispersión en la muestra que dirija algunos hacia direcciones diferentes con 
mayor o menor peso. Ásumimos que el ángulo de incidencia de los portadores 
obedecerá a una función de distribución de probabilidad expresada por: 


2 
P(6) =cexp (E (3:27) 
2a? 

El parámetro a de la función gaussiana determina la dispersión de la mues- 
tra abriendo o cerrando la distribución en torno a dy y produce la distinción 
entre una muestra bastante homogénea (una campana gaussiana estrecha, 
a pequeño) y una muestra con muchas inhomogeneidades (una distribución 
ancha, a grande). Se puede esperar transporte balístico en el primer caso, 
muy centrado en torno a una cierta dirección preferente Hp, y un colectivo 
disperso en el segundo. 

El parámetro c, que determina la altura de la gaussiana depende de a, 
ya que ambos deben verificar la condición de normalización Eq[3.28]. Como 
vemos en dicha expresión, a medida que c se hace cada vez menor, a se hace 
cada vez mayor y viceversa. 


T 
acV27 Er f (5) T (3.28) 
Cuando se considera este efecto, la conductancia efectiva del GFET Eq[3.23] 
debe ser ligeramente modificada multiplicando la contribución de cada direc- 
ción incidente por su peso estadístico. La conductancia efectiva para una 
distribución gaussiana G, queda entonces, considerando dy = 0? (cuando el 
coeficiente de transmision es la unidad): 


2c 





G¿(Ep, Up, D, a, c) 7 la 


(3.29) 


hvpr 


+7/2 ¿e 
Enf T(Ep,Uo, D, P) cos $ exp (-S) do 
—7/2 


64 


Conductancia 





La conductancia efectiva G¿ muestra ahora un comportamiento diferente 
según se aprecia en las figuras Fig. 3.6-3.7. 
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Figura 3.6: Conductancia frente a la energía de Fermi para una distribu- 
ción gaussiana en la dirección de incidencia centrada en torno a Hp = 0*, 
de anchura a = 0.64 (rojo), 0.21 (verde) y 0.08 (azul). Uy = 200 meV y 
D = 100 nm. La conductancia crece con oscilaciones de manera aproximada- 
mente lineal, hasta alcanzar un máximo local que depende de los parámetros 
utilizados. Posteriormente decae hasta Ep =— Up, para más adelante regresar 
a un crecimiento lineal. 


La figura Fig. 3.6 muestra el comportamiento de la conductancia cuando 
varía la energía de Fermi mientras que las otras variables permanecen fijas. 
Inicialmente la conductancia aumenta a la vez que oscila hasta que alcan- 
za un máximo local. El valor y la posición de este máximo depende de los 
parámetros de la función gaussiana utilizada. Posteriormente, la conductan- 
cia disminuye hasta que se alcanza un mínimo local cuando Ep “= Uy. Una 
vez superado dicho mínimo, la conductancia se incrementará monótonamen- 
te. En la figura Fig. 3.7, cuando utilizamos Uy como variable mientras que los 
otros parámetros permanecen fijos en valores constantes, podemos ver que 
la conductancia efectiva disminuye hasta que alcanza un mínimo en torno a 
Us -— Ep. También observamos que en funciones gaussianas más estrechas 
(lo cual es equivalente a una distribución con menos desorden), la variación 
de la conductancia es más pronunciada. Una vez que Up supera el valor de 
Ep, observamos oscilaciones en la conductancia, cuya frecuencia tiende hacia 
valores w =2D/Rhvp. 
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Figura 3.7: Conductancia frente a la altura de la barrera para una distribución 
gaussiana en la dirección de incidencia centrada en torno a dy = 0%, de 
anchura a = 0.64 (rojo), 0.21 (verde) y 0.08 (azul). Ep = 100 meV y D = 


100 nm. Presenta un mínimo en torno a Uy - Ep, que es tanto más profundo 
cuanto mayor sea el parámetro a. 




















E A ] 
7 0.3: x 50 ——0.64----- 150 
a PA 50 ——(0:21--=-- 150 
a 02 5 50 ——0.08----- 150 
y E "EA A AN A A A 

0.1f 

DOBLE > o, A ea oa 

0 50 100 150 200 
D (mm) 


Figura 3.8: Conductancia frente a la anchura de la barrera para una distri- 
bución gaussiana en la dirección de incidencia centrada en torno a dy = 0*, 
de anchura a = 0.64 (rojo), 0.21 (verde) y 0.08 (azul). Ep = 50 meV (línea 
continua) y 150 meV (línea discontinua). Ug = 100 meV. La conductan- 
cia presenta oscilaciones en torno a un valor constante, que es tanto menor 


cuanto mayor sea a. 
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La dependencia de la conductancia efectiva con la anchura de la barrera 
D para una distribución gaussiana centrada en 0% se muestra en la figura 
Fig. 3.8. Básicamente, la conductancia presenta oscilaciones en torno a un 
valor constante de la conductancia que se atenúan para barreras muy an- 
chas. Para un valor dado de la energía de Fermi, a medida que el parámetro 
a disminuye, el valor promedio que alcanza la conductancia aumenta, res- 
pondiendo al hecho de que cuanto más aguda sea la distribución incidente 
de electrones, mayor importancia tendrá en la conductividad el efecto Klein. 


Consideremos ahora el caso en el cual la distribución gaussiana está cen- 
trada en torno a dy 4 0”. Este ángulo Hp no nulo puede realizarse inclinando 
la orientación de la puerta superior en un dispositivo basado en grafeno co- 
mo se muestra por ejemplo en [87]. En este caso dy 4 0”, la condición de 
normalización está dada por la expresión Eq[3.30] 


A 


Como vemos en las figuras Fig. 3.9 y Fig. 3.10, los resultados son se- 
mejantes a áquellos obtenidos en el caso anterior dp = 0”, pero los mínimos 
observados son ahora más anchos. De hecho, cuando a tiende a ser muy próxi- 
ma a cero los mínimos llegarán a ser muy anchos. Se aprecia que a medida que 
la inhomogeneidad del medio disminuye junto con la dispersión de la distri- 
bución, el mínimo de conductancia se va haciendo más ancho, observándose 
así la apertura de un pseudogap de energía en que la conductancia decae 
abruptamente. 

La dependencia de la conductancia efectiva con la anchura de la barrera 
D para una distribución gaussiana centrada en 22.5” se muestra en la figura 
Fig. 3.11. La forma que adquiere la conductancia en este caso es bien diferente 
para los dos valores de energía tomados como parámetros. Para Hp = 50 meV 
presenta oscilaciones como en la figura Fig. 3.8, pero de amplitudes diferentes. 
Sin embargo, para Ep = 150 meV apenas presenta oscilaciones. 


3.5.2. Distribución delta 


Consideramos a continuación el escenario extremo, en el cual la función 
gaussiana es tan estrecha que puede ser representada mediante una función 
delta Eq[3.31]. Representa el caso en el cual el material es ideal y todos los 
portadores impactan en la barrera con el mismo ángulo de incidencia, do. 


P(6) = 4(6 — bo) (3.31) 
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Figura 3.9: Conductancia frente a la energía de Fermi para una distribución 
gaussiana en la dirección de incidencia centrada en torno a dp = 22.5”, de 
anchura a = 0.64 (rojo), 0.21 (verde) y 0.08 (azul). Uy = 200 meV, D = 
100 nm. Para ángulos considerables, el mínimo se vuelve tanto más ancho 
cuanto menor sea el parámetro a. 
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Figura 3.10: Conductancia frente a la altura de la barrera para una distri- 
bución gaussiana en la dirección de incidencia centrada en torno a (y = 
22.5”, de anchura a = 0.64 (rojo), 0.21 (verde) y 0.08 (azul). Ep = 100 meV, 
D = 100 nm. Para ángulos considerables, el mínimo en torno a Uy -— Ep se 
puede convertir en un pseudogap dependiendo del parámetro a utilizado. 
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Figura 3.11: Conductancia frente a la anchura de la barrera para una dis- 
tribución gaussiana en la dirección de incidencia centrada en torno a dy = 
22.5”, de anchura a = 0.64 (rojo), 0.21 (verde) y 0.08 (azul). Ep = 50 meV 
(línea continua) y 150 meV (línea discontinua), Uy = 100 meV. La conduc- 
tancia tiende a un valor constante mediante oscilaciones o no, dependiendo 
del parámetro Ep. 


En este caso, la conductancia efectiva está dada por: 


21E 
Gs.óo Ep, Uo, D, do) = (Ep, Un, D, 6) cos Po (3.32) 





Observamos que ahora la conductancia efectiva depende del ángulo de 
incidencia, así como del resto de los parámetros considerados anteriormente 
(energía de Fermi Ep, altura de la barrera Uy y anchura de la barrera D). 


En particular, si la corriente fluye perpendicularmente a la barrera, el coe- 
ficiente de transmisión es igual a la unidad con independencia de la energía 
de Fermi y cualquiera que sea el valor de los otros parámetros de la ba- 
rrera, como consecuencia de la paradoja de Klein [57,62,85). Por tanto, la 
conductancia a Hp = 0” es directamente proporcional a la energía de Fermi 


2 
Gao = [Er] (3.33) 
Ur 


Para incidencia no perpendicular, la conductancia presenta un pseudogap 
efectivo entre los valores de la energía de Fermi 
Vo Do 


KK < Ep < 3.34 
l+sindo — a 1 —sinódo ) 
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obteniéndose un mínimo efectivo para Ep = Uy/ cos? Hp (como mostramos 
en la figura Fig. 3.12), cuyo valor es 


hvup 


20 DU0t 
Gs,$0, min = e cosh”? (na) cos”* bp (3.35) 
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Figura 3.12: Conductancia efectiva frente a la energía de Fermi para una 
distribución tipo delta cuando la altura de la barrera es U¿ = 200 meV y 
su anchura D = 100 nm. Suponemos diferentes ángulos de incidencia : dp = 
0”(rojo), 2”(verde) y 5”(azul). Estas curvas son el caso límite de una gaussiana 
estrecha. A medida que aumenta o, el mínimo se hace más profundo, excepto 
por esto las gráficas son aproximadamente lineales. 


En las figuras Fig. 3.12-3.14 representamos la conductancia efectiva para 
una distribución tipo delta y ángulos de incidencia extremadamente pequeños 
dy = 0%, 22 y 5%. En la figura Fig. 3.12 frente a la energía de Fermi, en la 
figura Fig. 3.13 frente a la altura de la barrera de potencial y en la figu- 
ra Fig. 3.14 frente a la anchura de la barrera. El propósito de todo ello es, 
conociendo que T' = 1 para do = 0”, observar cómo al modificar el ángu- 
lo de entrada de los electrones en las figuras Fig. 3.12-3.13 va apareciendo 
un mínimo en la conductancia para una energía determinada. En la figura 
Fig. 3.14, la dependencia de la conductancia con la anchura de la barrera no 
muestra características destacables excepto por una pequeña oscilación cuya 
amplitud depende de la energía de Fermi. 


En este escenario, en el cual asumimos una función de distribución de pro- 
babilidad tipo delta, podemos observar como un hecho muy notable en las 
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Figura 3.13: Conductancia efectiva frente a la altura de la barrera asumien- 
do una distribución tipo delta cuando la energía de Fermi es Ep = 100 meV 
y la anchura de la barrera es D = 100 nm. Estas representaciones suponen 
diferentes ángulos de incidencia dy = 0”(rojo), 2” (verde) y 5”(azul). Las gráfi- 
cas presentan un valor aproximadamente constante excepto para Uy = By, 
donde aparece un mínimo, tanto más profundo cuanto mayor sea Qp. 
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Figura 3.14: Conductancia efectiva frente a la anchura de la barrera asumien- 
do una distribución de probabilidad tipo delta para energías de Fermi Ep = 
50 (línea continua) y 150 meV (línea discontinua). La altura de la barrera de 
potencial es Uy = 100 meV, además hemos supuesto diferentes ángulos de 
incidencia : fp = 0*(rojo), 2”(verde) y 5”(azul). Las gráficas toman valores 
aproximadamente constantes dependiendo básicamente de Ep. 
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Figura 3.15: Conductancia efectiva frente a la energía de Fermi cuando asu- 
mimos una función de distribución de probabilidad tipo delta y un ángulo 
de incidencia fp = 22.5”. La barrera de potencial se caracteriza por Up = 
200 meV y D = 100 nm. Aparece claramente un pseudogap, con valores 
prácticamente nulos en él. 
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Figura 3.16: Conductancia frente a la altura de la barrera Uy cuando asumi- 
mos una función de distribución de probabilidad tipo delta y un ángulo de 
incidencia fp = 22.5”, siendo los otros parámetros Ep = 100 meV y D = 100 
nm. El mínimo centrado en Uy - Ep toma valores aproximadamente nulos. 
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Figura 3.17: Conductancia frente a la anchura de la barrera asumiendo una 
función tipo delta centrada en torno a Hp = 22.5” para Uy = 100 meV, siendo 


Ep = 50(línea continua) y 150 meV (línea discontinua). Aparecen oscilaciones 
o no dependiendo de Er. 


eráficas lo que se acaba de describir analíticamente, la existencia de regiones 
de energía en las cuales la conductancia disminuye hasta valores próximos a 
cero. Se observa también que a medida que aumenta el ángulo de incidencia, 
este mínimo se hace cada vez más pronunciado y la región más extensa. En 
otras palabras, parece que el control de la direccionalidad de los electrones 
de entrada en el dispositivo puede utilizarse para provocar la apertura de 
un pseudogap efectivo en el GFET. En las figuras Fig. 3.15 y 3.16 para un 

ángulo de incidencia de dy = 22.5 "(o 7/8) se hace patente este efecto con la 
aparición del mínimo en la expresión Eq[3.35], con la conductancia práctica- 


mente nula entre los dos valores extremos de la energía determinados en la 
relación Eq[3.34]. 


Es interesante destacar que tanto la anchura del pseudogap como su pro- 
fundidad, pueden controlarse a voluntad con la altura de la barrera (relacio- 
nada con el voltaje aplicado en la puerta superior en el esquema del transis- 
tor mostrado en la figura Fig. 3.1), la anchura de la barrera y la dirección 
de entrada de los electrones en el dispositivo. El mínimo observado en la 
conductancia y la apertura del pseudogap como consecuencia, en la figura 
Fig. 3.15, evidencia que en ciertas condiciones el efecto túnel de Klein puede 
ser evitado, y que la construcción de un GFET efectivo basado en grafeno 
con una notable relación encendido/apagado podría ser factible. Creemos que 
una puerta superior no-perpendicular a la dirección fuente-drenador, puede 
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crear este pseudogap en una muestra de grafeno limpia y homogénea, sin 
impurezas [87]. También es posible controlar la dirección de incidencia de los 
portadores mediante guiado eléctrico como recoge Kim en [88]. 


3.6. Resumen y conclusiones 


A lo largo de este capítulo hemos construido una expresión cerrada para 
describir la densidad de corriente y la conductancia de los GFET conside- 
rados. Esta expresión relaciona el coeficiente de transmisión, la densidad de 
corriente y la conductancia del dispositivo con los parámetros físicos que go- 
biernan el modelo. Hemos considerado varias simplificaciones en la densidad 
de corriente, que reproducen ciertas condiciones experimentales. En parti- 
cular, baja polarización y baja temperatura para poder utilizarlo con un 
GFET. Las variables experimentales que controlan el coeficiente de transmi- 
sión, la corriente y la conductancia son los potenciales externos aplicados a 
sus terminales de: surtidor, drenador, puerta superior y puerta inferior. Éstos 
determinan el perfil de potencial y la energía de Fermi de los portadores de 
carga. 


Con objeto de mejorar el rendimiento del GFET descrito y para evitar 
las propiedades triviales que condicionan un impacto perpendicular sobre la 
barrera de los portadores de corriente (efecto de Klein), hemos considerado 
un impacto no-perpendicular. Esta direccionalidad en la transmisión nos lle- 
va a nuevas propiedades de la transmisión en el dispositivo. Hemos supuesto 
varios escenarios en la transmisión. Éstos se basan en un modelado de dife- 
rentes condiciones de dispersión en la muestra de portadores incidentes, que 
pueden ser atribuidas a diferentes grados de inhomogeneidad presentes en el 
medio (grafeno). 


Atribuimos diferentes pesos estadísticos en el transporte a los electrones 
que impactan en una cierta dirección, y modelamos el colectivo de electro- 
nes entrantes utilizando diferentes funciones de ponderación, una función 
gaussiana primero (describiendo un gran número de inhomogeneidades en el 
medio), y una función delta después, que selecciona una única dirección de 
entrada Hp (sin inhomogeneidades) y corresponde a un haz incidente perfec- 
tamente colimado. Al analizar la conductancia asociada, hemos encontrado 
en este escenario que la paradoja de Klein puede ser soslayada induciéndose 
la creación de un pseudogap efectivo. La posición y la anchura del pseu- 
dogap de energía dependen de la altura de la barrera de potencial y de la 
dirección fp de entrada de los electrones. Esta característica podría ser un 
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importante paso en el desarrollo de transistores basados en grafeno, siempre 
que podamos controlar la direccionalidad, ya que la ausencia de gap es una 
de las principales dificultades encontradas cuando tratamos de modular la 
corriente efectiva via la polarización. Creemos que una puerta superior no- 
perpendicular a la dirección fuente-drenador, puede crear este pseudogap en 
una muestra de grafeno limpia y homogénea, sin impurezas. Otra posibilidad 
para controlar la dirección de incidencia de portadores es mediante guiado 
eléctrico. 
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Capítulo 4 


Transmisión a través de 
múltiples barreras de potencial 


4.1. Introducción 


Una vez estudiado el coeficiente de transmisión a través de una y dos 
barreras de potencial. Consideramos a continuación el cálculo del coeficiente 
de transmisión para varias barreras de potencial mediante el formalismo de 
Dirac, utilizando el método de la matriz de transferencia. El método de la 
matriz de transferencia (MMT) fue ampliamente utilizado en el transporte 
electrónico a través de heteroestructuras constituidas por semiconductores 
tradicionales tales como GaAs/Al,(Gaj-.As, caracterizadas por la ecuación 
de Schródinger, tal y como recoge Liu en [89,90]. El método de la matriz de 
transferencia es un método sencillo para calcular el coeficiente de transmi- 
sión a través de múltiples barreras de potencial en una monocapa de grafeno. 
Se trata de un método numérico y analítico muy potente para justificar las 
propiedades de transporte en los sistemas, y además es computacionalmente 
muy fácil de implementar, ya que únicamente hay que considerar la multi- 
plicación de matrices de dimensiones 2 x 2. 


Pretendemos estudiar el comportamiento tanto de varias barreras de po- 
tencial, como de una barrera de potencial escalonada. Las primeras pueden 
llegar a fabricarse mediante un conjunto de electrodos aplicados al grafeno 
a través de los dieléctricos correspondientes [91], con un diseño constructi- 
vo semejante al utilizado para la fabricación de transistores de grafeno. La 
finalidad es utilizar varias puertas para obtener nanodispositivos versátiles, 
que puedan ser usados para distintas aplicaciones sin necesidad de cambiar 
su configuración. La segunda puede ser utilizada para simular el compor- 
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tamiento de potenciales que varíen suavemente [58,92]. En nuestro caso lo 
aplicaremos a la aproximación para una barrera gaussiana. 


Los resultados numéricos muestran que la probabilidad de transmisión 
depende de factores geométricos como las anchuras y las distancias de sepa- 
ración entre barreras así como del valor de potencial aplicado a las diferentes 
puertas (barreras de potencial), también de las características propias de los 
portadores como energía y ángulo de incidencia sobre la barrera de poten- 
cial, así como del número de barreras de potencial. Por ello cambiando la 
configuración de la estructura, las propiedades de transporte de los electro- 
nes pueden ser ajustadas adecuadamente a la aplicación que se pretenda. En 
particular, cabe destacar la aparición de pseudogaps con resonancias muy 
selectivas en energía cuando consideramos barreras de potencial de la misma 
altura. Así como una simetría especular y la posibilidad de ajustar indepen- 
dientemente los pseudogaps cuando las barreras sean de diferente altura. 


4.2. Múltiples barreras de potencial 


Consideremos múltiples barreras de energía potencial U;(x) en una di- 
mensión aplicadas a una lámina de grafeno monocapa, tal y como aparece 
en la figura Fig. 4.1. En dicha figura hemos detallado los coeficientes de ca- 
da una de las regiones de transmisión, así como su energía potencial. Una 
barrera de potencial genérica U;¡ estará comprendida entre Za;-1 y T2a¡, CON 
coeficientes asociados A9; y Boa;. La región a su izquierda vendrá caracteri- 
zada por coeficientes A3;-1 y Ba¡-1, y la región a su derecha por coeficientes 
Azi+r1 y Bois: 

Recordando que 4 = vpd -¿+U(a)Í, y Yi (a,y) = vila) - e*vY. La 
expresión general para una región libre, es decir donde no tengamos aplicada 
una barrera de potencial o bien sea nula, por ejemplo la zona 21-1 es: 


Viale) = Agiet? 4 Baje? (4.1) 


va (1) = Soj-1( Ag 10 Morte] — Bo; ¡e tor+tól) 


y para una región 2i donde tengamos aplicada una barrera de potencial 
U; no nula es: 


v4(x) = Aye Pie? y Bye Pio? (4.3) 


vila) = sal Ageitiortól — poe ile m+0a)) (4.4) 
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Aa. X3 X4 X2n-1 X2n 


Figura 4.1: Estructura general para múltiples barreras de potencial a lo largo 
de la dirección 2. 











donde 
E : 
k= ==. Ri= kcosto). y =ksinlo), (4.5) 
hup 
E=UNX" k 
Qi, = / e ) - k, 07, = arctan (E) : (4.6) 
s=sign(E), sa = sign(E — U;), (4.7) 


siendo U; la energía potencial de la barrera i-ésima. 
De manera general, al atravesar los portadores la barrera i-ésima impo- 
nemos: 


vá_,(2o;-1) => v4(29,-1), V%-1 (2-1) = vE(taj-1), (4.8) 
vál(ea) = Válua), válea) = VE,1 (025) (4.9) 


lo cual se traduce en las siguientes expresiones: 


Ao e k=T2i1 de Boj ¡e Per2i1 — Ao¡e BioT2i1 JE Baje “Riotaia 
s(Ag¡-e tona1+0) _ Baj_-1e erzi+el) 

pl Sa (Aye Mi 72i1+021] a Baje Wist2i1 +0) 
Aoje Pio T2 + Baje Pio? e Aoijje hera + Boga eN 
So (Aye ia 724001] _ Baje Mio t2a +02) 


= s(Agiy e teratél _ Boj, 0 oratel) (4.10) 
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obteniendo así las matrices: 


A e oTai—1 e kota1 

(P, ) a se lket2-1+0l —se Hhorai1+el) > (4.11) 

1924 q U2i—1 —10 2: y 1241 

(P; ) o ( $90 Mio 02i-1+02i] —sye Mio 2140211) ) > (4.12) 
os e iz Tai e ip Tai 

ee ) E ( soe Mis 22i+02:] —syje HU2io va +Óail) > (4.13) 

¡ko 22; —ikot25 

4 e' e 

(P; ) => ( setlkxt2+0] —se hoz tel) ) ) (4.14) 

resultando 


Ao; Asi 
Pi ( W ) ; 4.15 
E) E a 
P? ( es ) =P? ( a ) ? (4.16) 


por lo que al atravesar la barrera ¿ésima de potencial obtenemos: 





(Mo) = (POT PD PDT (eS (4.17) 


En general para n barreras de potencial, podemos relacionar los parámetros 
de la primera y de la última región mediante un conjunto de matrices 


Cs) 


La ecuación anterior relaciona los coeficientes del estado de salida con los 
del estado de entrada. Estamos interesados en describir la situación en que 
entra un electrón en un estado puro propagándose hacia la derecha, parte 
de él se refleja y otra parte se va transmitiendo por las regiones intermedias, 
donde sufre reflexion y transmisión en las distintas superficies de separación 
hasta que emerge transmitiéndose hacia la derecha. 

Entonces los valores que deben tomar los coeficientes en las regiones ex- 
tremas son: 


(Mi) - (Ma) - (Ma) --(M,,) ( os ) 


= Mi) (Gt) (4.18) 





Ar — 1, B; =", A9n+1 = t, Boan+1 =0 (4.19) 


donde hemos supuesto que incide completamente toda la onda en la primera 
región A; = 1, y en la última región no puede haber reflexión, ya que es una 
zona libre Ba, 1 = 0. Los coeficientes de reflexión y transmisión son r y t. 
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ONO 


De aquí obtenemos el coeficiente de transmisión 


ol 
Mi 





t (4.21) 
y puesto que en las regiones extremas no existen barreras de potencial, la 
probabilidad de transmisión es: 


T=1-f (4.22) 


4.2.1. Barreras con idéntica energía potencial 


Continuando con la notación, la anchura de la barrera de potencial i-ésima 
es L; y la distancia entre la barrera de potencial anterior y la siguiente es 
D;. En la figura Fig. 4.2 representamos el coeficiente de transmisión frente al 
ángulo de incidencia para 1 (a), 3 (b) y 5 (c) barreras de potencial, todas ellas 
con la misma energía potencial Uy = 0.4 eV, y misma anchura L = 20 nm, con 
distancia de separación entre las barreras D = 30 nm, para una energía del 
portador incidente E = 0.2 eV. Observamos que el coeficiente de transmisión 
es diferente dependiendo del número de barreras consideradas. Á medida que 
aumenta el número de barreras de potencial, la anchura angular de la zona 
cercana a $ = 0” donde T' es próximo a la unidad, es más estrecha, ya que la 
disminución del coeficiente de transmisión comienza antes. Alguna resonancia 
que observamos para el caso de una barrera de potencial, se convierte en 
una zona con múltiples resonancias muy estrechas al aumentar el número 
de barreras de potencial, así como la aparición de nuevas resonancias tipo 
delta [93]. 

En la figura Fig. 4.3 tenemos cinco barreras de potencial idénticas - con 
la misma altura de barrera Uy, y la misma anchura 20 nm y misma distan- 
cia de separación 30 nm. Fijamos la energía de los electrones incidentes en 
dichas barreras de potencial en 0.1 eV, y con ello realizamos la gráfica de 
contorno del coeficiente de transmisión en función de la energía de la barrera 
de potencial Uy y del ángulo de incidencia f. Para incidencia perpendicular 
el coeficiente de transmisión 7' es la unidad, lo cual sabíamos que sucedería 
debido a la paradoja de Klein, también es unitario en las proximidades de 
$ = 0, sin embargo en la zona angular opuesta y = T/2 es nulo. En ciertos 
intervalos del ángulo de incidencia (en torno a 37/16) aparecen alternativa- 
mente áreas (dependiendo de la altura de la barrera de potencial U¿) con 


81 


Múltiples barreras 








1 barrera 


0.2 - 








0.01, 





o 
013 + 
2/3 
oz! 
NIS 


P(rad) 





3 barreras 


0.21 
(b) 








0.01, 





[) 
oa? 
23 
cy! 
NE 


ó(rad) 





5 barreras 


(c) 











T 3x 


4 3 
ó(rad) 


St 
00/13 
Ia E 


Figura 4.2: Coeficiente de transmisión frente al ángulo de incidencia para 1 
(a), 3 (b) y 5 (c) barreras de potencial, siendo E = 0.2 eV, Up = 0.4 eV, L 
= 20 nm y D = 30 nm. A medida que aumenta el número de barreras de 
potencial, la zona angular próxima a la incidencia perpendicular donde T se 
aproxima a 1, disminuye. 
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Figura 4.3: Gráfica de contorno para el coeficiente de transmisión de 5 barre- 
ras de potencial idénticas en función de Uy y del ángulo de incidencia, siendo 
E=0.1eV, L =20 nm y D = 30 nm. Para múltiples barreras de potencial 
idénticas aparece la formación de minibandas cuando los electrones inciden 
de manera oblicua. 


fuertes oscilaciones en T' para un ángulo de incidencia df fijo. Para un ángu- 
lo de incidencia ligeramente mayor, cada zona de transmisión se convierte 
en un conjunto de pequeñas minibandas. Al desplazarnos hacia ángulos aún 
mayores, algunos picos en la transmisión permanecen hasta que desaparecen 
para ángulos muy elevados. Así, para múltiples barreras de potencial idénti- 
cas hemos observado la formación de minibandas de energía para electrones 
que inciden de manera oblicua sobre las barreras de potencial. 

En la figura Fig. 4.4 representamos la gráfica de contorno para el coefi- 
ciente de transmisión de un sistema constituido por 3 barreras de potencial 
de la misma altura Uy = 0.4 eV, misma anchura L = 10 nm y misma dis- 
tancia de separación D = 40 nm, en función de la energía y del ángulo del 
portador incidente. Para un ángulo fijo fp, claramente se observan zonas al- 
ternas de valor mínimo y máximo en el coeficiente de transmisión, a medida 
que aumenta la energía del portador incidente. Además observamos que el 
coeficiente de transmisión toma valores elevados para ángulos de incidencia 
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más grandes cuando la energía E es inferior a Us. 
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Figura 4.4: Gráfica de contorno para el coeficiente de transmisión de 3 barre- 
ras idénticas de potencial en función del ángulo de incidencia y de la energía 
del portador incidente, siendo Uy = 0.4 eV, L = 10 nm y D = 40 nm. Aparece 
la formación de zonas alternas de valor máximo y mínimo a un ángulo fijo, 
a medida que aumentamos la energía del portador incidente. Para E < Up 
existen ángulos grandes con coeficiente de transmisión elevado. 


Habíamos visto con anterioridad en la Sección 2.2 que cuando represen- 
tamos el coeficiente de transmisión frente a la energía del electrón incidente 
para una barrera de potencial aparecía un pseudogap. Cuando el número de 
barreras idénticas aumenta, el intervalo de este pseudogap es básicamente el 
mismo. Para una barrera de potencial, como habíamos visto, el intervalo de 
energías en el cual tenemos este pseudogap es: 


Ú U 
Erin = AE. PEA Bras = == (4.23) 
l1+sing 1 — sing 
por lo que la anchura energética del pseudogap es: 
Foff = Emaa— Emin = 2 Uy sinó cos 26 (4,24) 
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Así observamos que la anchura energética del pseudogap E,fy en la ecuación 
Eq[4.24] solamente depende de la altura de la barrera de potencial y, en ma- 
yor medida del ángulo de incidencia de los electrones. O sea, es independiente 
de la anchura de las barreras de potencial y de la distancia de separación en- 
tre ellas. Este intervalo energético E,f¿ depende más fuertemente del ángulo 
incidente /, que de la altura de la barrera de potencial Up. En efecto, puesto 
que sin varía entre 0 y 1, Emin varía entre U¿/2 y Up, mientras que Ema 
varía entre Up e infinito. Es decir, la variación de Ein está acotada, pero la 
variación de Emaz no. De ahí la variación tan acusada del intervalo energético 
E.¿¡ sobre todo con el ángulo de incidencia d. 


Puesto que la energía de barrera es idéntica en todas ellas, U,, el aumento 
en el número de barreras de potencial no influye sobre el intervalo energético 
del pseudogap, sin embargo produce un aumento en el número de oscilacio- 
nes y en el hecho de que sean más pronunciadas, en particular en la frontera 
entre el intervalo energético del pseudogap y las energías próximas a él. Lo 
anterior queda patente en la figura Fig. 4.5 (a) para 1 (línea punto-raya), 3 
(línea discontinua) y 5 (línea continua) barreras de potencial, considerando 
como parámetros Uy = 0.4 eV, p = 71/10, L = 20 nm y D = 30 nm. También 
cabe destacar que una diferencia importante entre el caso de una barrera de 
potencial y el de múltiples barreras de potencial es la aparición en este último 
caso de unas pocas resonancias (tipo delta) altamente selectivas en energía 
dentro del intervalo energético del pseudogap, lo cual también observamos 
en la figura. 


A continuación si se aumenta el ángulo de incidencia hasta alcanzar el 
valor 7/5 como representamos en la figura Fig. 4.5 (b), siendo el resto de 
parámetros los mismos que teníamos en (a). Observamos que ahora las re- 
sonancias tipo delta en el intervalo energético del pseudogap parece que han 
desaparecido. Así mismo, aumenta la anchura del intervalo energético del 
pseudogap (como habíamos establecido según la ecuación Eq[4.23]). Además 
las oscilaciones para el coeficiente de transmisión son más pronunciadas, de 
tal forma que algunos mínimos alcanzan en su nivel inferior el valor de cero, 
convirtiéndose así en nuevos pseudogaps. En particular se observa claramen- 
te para el caso de cinco barreras idénticas y ángulo de incidencia 7/5 en la 
figura Fig. 4.5 (b). 
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Figura 4.5: Coeficiente de transmisión frente a la energía del portador inci- 
dente para 1 (línea punto-raya), 3 (línea discontinua) y 5 (línea continua) 
barreras de potencial, siendo Uy = 0.4 eV, L=20 nm, D =30 nm, p =7/10 
(a) y 7/5 (b). El aumento del número de barreras de potencial no influye 
sobre el intervalo energético del pseudogap, aparecen resonancias tipo delta 
dentro de él (a). Al aumentar el ángulo respecto a (a), aumenta el intervalo 
energético del pseudogap, las oscilaciones se hacen más profundas, aparecien- 
do nuevas zonas de pseudogap (b). 
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Figura 4.6: Coeficiente de transmisión frente al ángulo de incidencia para 
barreras de potencial de alturas U,, Uz, Uz = 0.2, 0.4, 0.6 eV (línea punteada), 
0.6, 0.4, 0.2 eV (línea discontinua), 0.2, 0.6, 0.4 eV (línea continua), siendo 
E=0.1eV y L= D = 100 nm. Solamente cuando el ordenamiento de las 
barreras de potencial de diferente altura, visto desde izquierda o derecha 
coincidan, se obtiene el mismo coeficiente de transmisión. 


4.2.2. Barreras con diferente energía potencial 


Los resultados obtenidos en las expresiones Eq[4.11-4.22] pueden aplicarse 
también al caso de múltiples barreras con altura diferente. Creemos conve- 
niente aclarar que, en este caso, el orden de las barreras es significativo. En 
el caso particular de dos barreras de potencial el resultado es independien- 
te del orden en el cual las coloquemos, ya que portadores incidiendo desde 
sentidos contrarios, perciben igual el potencial. Sin embargo, si el número de 
barreras es tres o superior a tres, el coeficiente de transmisión dará el mismo 
resultado siempre que exista una simetría especular, es decir el orden en el 
cual colocamos las barreras de potencial mirando desde el extremo izquierdo 
coincida con el orden de colocación de dichas barreras cuando miramos desde 
el extremo derecho. 


Para tres barreras hay sólo tres disposiciones diferentes según se ilustra 
en la Fig. 4.6 donde representamos diferentes ordenes de colocación de tres 
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Figura 4.7: Coeficiente de transmisión frente a la energía del portador inci- 
dente para tres barreras de potencial de alturas U;,, Uz, Uz = 0.2, 0.4, 0.6 eV 
siendo L = D = 100 nm y 6 = 7/80 (a), 7/20 (b), 7/10 (c), 7/5 (d) radia- 
nes. En (a) observamos un mínimo asociado a cada barrera de potencial. Al 
aumentar el ángulo de incidencia /, los mínimos se convierten en pseudogaps 
(b). Al aumentar aún más (, los pseudogaps segundo y tercero, colapsan en 
un único pseudogap en (c). Al continuar con el aumento de df, sucede que 
todos los pseudogaps colapsan en un único pseudogap (d). 


barreras de potencial siendo E = 0.1 eV y L = D = 100 nm. Observamos 
que el caso de las barreras U¡, Uz, Uz = 0.2, 0.4, 0.6 eV (línea punteada) da 
el mismo resultado que el caso U;¡, Uz, Uz = 0.6, 0.4, 0.2 eV (línea disconti- 
nua). Ambos casos son diferentes al caso U,, Uz, Uz = 0.2, 0.6, 0.4 eV (línea 
continua) y también al caso U,, Uz, Uz = 0.4, 0.2, 0.6 eV. Por ello, para tres 
barreras de potencial de diferentes alturas solamente existen tres casos de 
ordenamiento diferente, y tres coeficientes de transmision distintos. 


Cuando tenemos barreras de potencial con diferentes alturas y represen- 
tamos el coeficiente de transmisión frente a la energía del electrón incidente, 
en principio obtenemos una zona de transmisión mínima asociada a cada una 
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de estas barreras de potencial, como podemos ver en la figura Fig. 4.7 (a). 
En dicha gráfica representamos tres barreras de potencial de valores U;, Uz, 
Uz = 0.2, 0.4, 0.6 eV coincidiendo la anchura de la barrera L con la distan- 
cia de separación entre barreras D, ambas de 100 nm y siendo el ángulo de 
incidencia $ = 7/80 radianes. Observamos que los valores de U,, Uz y $ que 
hemos asignado, no alcanzan un valor suficientemente elevado como para que 
obtengamos un pseudogap en el coeficiente de transmisión asociado a cada 
una de estas tres barreras de potencial. Es decir, los mínimos no llegan a caer 
en su parte inferior a un valor prácticamente nulo. En (b) conservamos los 
parámetros que teníamos, y únicamente modificamos el ángulo de incidencia, 
aumentándolo a un valor $ = 7/20 radianes. Ahora con un ángulo de inciden- 
cia mayor, observamos tres pseudogaps perfectamente formados, cada uno de 
los cuales está asociado a una barrera de potencial. Los valores de energía 
inicial y final de dichos pseudogaps están dados por la expresión Eq|[4.23], 
siendo ahora U, la altura de cada barrera U;. El límite superior del pseudogap 
varía mucho más que el límite inferior del pseudogap, con respecto al ángulo 
de incidencia. Debido a esto, puede llegar a suceder que el límite superior 
de un pseudogap alcance, o supere, el límite inferior en energía del siguiente 
pseudogap. Entonces, ambos pseudogaps colapsarán en un único pseudogap. 
Si mantenemos fijas las diferentes alturas de las barreras, esto sucederá para 
un cierto ángulo de incidencia. En (b) esto sucederá si aumentamos el ángulo 
de incidencia hasta y = 7/15.6 radianes, y para ángulos superiores a éste. 
El nuevo pseudogap es más ancho energéticamente, tal y como vemos en la 
figura (c). En dicho caso hemos representado el coeficiente de transmisión 
para $ = T/10 radianes, manteniendo el resto de los parámetros anteriores. 
Observamos que el segundo y tercer pseudogap han colapsado en un único 
pseudogap. Así, si seguimos con este proceso de ir incrementando el ángulo 
de incidencia, llegaremos a tener un único pseudogap, como representamos 
en (d). En dicha figura, el ángulo de incidencia alcanza $ = 7/5 radianes, 
manteniendo el resto de los parámetros anteriores, con los cuales venimos 
trabajando. Para este único pseudogap, el extremo inferior está asociado a 
la barrera de menor energía potencial, y el extremo superior a la barrera de 
mayor energía potencial, ya que el ángulo de incidencia es el mismo. Si sola- 
mente hubieramos considerado una barrera de potencial de altura 0.6 eV y 
$ = 1/5 - en vez de tres barreras de potencial de diferentes alturas -, el resul- 
tado hubiera sido un pseudogap con el mismo extremo energético superior, y 
distinto extremo energético inferior. Para que los extremos de un pseudogap 
se puedan ajustar independientemente, son necesarias al menos, dos barreras 
de potencial diferentes. 

En la figura Fig. 4.8 representamos el coeficiente de transmisión para 
tres barreras de potencial, dos de ellas a distintas alturas U¡ = 0.2 eV, U» 
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Figura 4.8: Gráfica de contorno representando el coeficiente de transmisión 
para tres barreras de potencial de idénticas dimensiones físicas L = 20 nm, 
D = 30 nm, siendo U, = 0.2 eV, Uz = 0.4 eV, y E = 0.1 eV. Para ángulos 
de incidencia f próximos a O radianes el coeficiente de transmisión es máxi- 
mo, mientras que en el extremo contrario es mínimo. Existe un intervalo 
(37/32, 107/32] radianes donde se alternan valores máximos y mínimos para 
el coeficiente de transmisión a medida que Uz aumenta. 


= 0.4 eV, y una tercera Uz sintonizable en un rango comprendido entre 0 
y 1 eV. Asumimos fija la energía del electrón incidente E = 0.1 eV, y las 
dimensiones físicas L = 20 nm y D = 30 nm idénticas para todas las barreras 
de potencial. En nuestra gráfica, las variables independientes son el valor de 
la altura de la tercera barrera de potencial así como el ángulo de incidencia 
del portador incidente (entre 0 y 7/2 radianes), obteniendo diferentes valores 
del coeficiente de transmisión que serán máximos (TP = 1) en torno a Y = 
O radianes, y nulos al aproximarnos al extremo contrario f = 7/2 radianes. 
Entre ambas zonas, aproximadamente en el intervalo de ángulos de incidencia 
comprendido entre 37/32 y 107/32 radianes observamos que existe cierta 
periodicidad en la gráfica asociada a los valores de U3, alternándose valores 
elevados y pequeños para el coeficiente de transmisión. Todo esto sucede, una 
vez que Uz ha superado la energía del electrón incidente. 
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4.3. Barrera de potencial escalonada 





Co C1 C2  C3 CnN-2 CN-1 CN 


Figura 4.9: Estructura general para múltiples barreras contiguas de potencial 
a lo largo de la dirección x. 


Consideremos ahora múltiples barreras (N) aplicadas a una lámina de 
grafeno monocapa. La altura de la barrera en una región entre posiciones 
contiguas c;_1 y cj es U, , según se ilustra en la figura Fig. 4.9. La región de 
entrada de los portadores es la región (0), y la de salida la región (N + 1), 
donde la altura de las barreras son nulas. El portador, que entra en la lámina 
por la zona (0) con energía E y momento (k,[0], k, + desplazándose hacia la 
derecha, viene descrito por una función de onda e**vY/;(x) cuya componente 
x es la dada por: 


If [KoL5] =0,1, Melide, /4_ 71] 
sI [ko[5] = 0,4, estat, /%, [7] 

IE [koi] = 0, eya, e=ikebilo,/¿T7] 
s¡Hf [Ko[i] =0,k,x — 1, —etelile, JE | 


(4.25) 
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siendo k+[j] = k. [5] + ik, y 











LS 
ky = a sin(9), (4.26) 
A ME-UN ,, 
ki] = y ( Em ) 2, (4.27) 
ss =stgn(B= Uy): Do = Un +1 =0. (4.28) 


Se ha introducido la condición en las funciones con el fin de evitar errores 
en el cálculo numérico para los valores de energía incidente para los que 
k« 13] = 0, en cuyo caso la E.D.O. que determina la solución se hace singular. 
En este caso, la solución es una función lineal en lugar de una exponencial [14], 
que resulta necesaria en el caso mencionado. Al atravesar los portadores la 
superficie que separa la región j de la ¿ + 1 se debe cumplir: 


Pilej) = Vjr1rles) (4.29) 


lo cual se traduce en las siguientes expresiones, que establecen la relación 
entre los coeficientes de la función de entrada sobre el límite x= = c; y los de 
salida: 





Pili] Prol] 
Aa za pea 0 ) ens o) Cta 0 ) Aj 
Bir1 0 ¿oli+ilez ) T£[ko [5+1]=0,1,2k [+1] o  e“belile; B, 
(4.30) 
siendo 
Pauli] = 1£ (ko(5]  =0,1, Vk- 07] 1 [x. pad di=te, +1] + 


A [Kol5] =0,1,-¿Vk+U 1] If [eli +1] =0,k,cj,¿V kl + 1] 
Po. TE [als =0, kyes, VR] If [els +1] =0,1-— k,c;, +1] se 

de [els] =0,kyc; — 1,1 11] If [ko li +1] =0,kycjp iv kaly + 1] 
Paljl = IE ol; ] =0,1, VA] 16 ko iio Rs [5 +1] m 

= asa lt [KoL5] E ik If [oli +1)=0,L iv ki + 1] 
Pal] = 1H [eli =0, k,C5, Vit If [eli E Vi +1] = 

e esté [els] O E ol If [eli 0 A EE y] 

(4.31) 
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La matriz de transmisión completa, que proporciona a través de sus coe- 
ficientes el estado de salida en la zona (N + 1) para un estado cualquiera de 
entrada en la zona (0), se obtiene a partir de las matrices anteriores como 


Án+1 ¿—ikp[N+1]e Ao 
( Bn+1 ) Ñ ( 0 él seed) ME Bo (4.32) 








con 
M = My: My-1 : My-2--:M23- Mi - Mo (4.33) 
y 
M, = 1 Pili Piel A | ¿olile[s=0,03.(e5=0j-10] 0 
J — Tflk2[¡+1]=0,1,2kz [41] | Por[i] Pozlj] d ¿—oliM[i=0,03,(ej=ej-1)] 
(4.34) 
con Det|[M] = Sit 


Las ecuaciones anteriores Eq[4.32-4.34] relacionan los coeficientes del es- 
tado de salida con los del estado de entrada. Estamos interesados en describir 
la situación en que entra un portador en un estado puro propagándose hacia 
la derecha, parte del cual se refleja y la otra parte se va transmitiendo por 
las regiones intermedias, donde sufre reflexión y transmisión en las distintas 
superficies de separación, hasta que emerge transmitiéndose exclusivamente 
hacia la derecha. Entonces, los valores que deben tomar los coeficientes en 
las regiones extremas son: 


Ao NN 1, Bo =", Án+1 = E, Bn+1 =0 (4.35) 


siendo r y t los coeficientes de reflexión y transmisión. 


t a 1 == Mii Mio 1 
al a 0 (4.36) 
y, como en las regiones extremas no existen barreras de potencial, resulta 


obviamente que Det|M] = 1 y obtenemos el coeficiente de transmisión total 
T' de la heteroestructura como 


l 


t=>—, 
Ma 


PIAR (4.37) 
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4.3.1. Barrera gaussiana 


Como aplicación del coeficiente de transmisión desarrollado en el punto 
anterior, consideramos a continuación una barrera de potencial que varíe 
suavemente, como es el caso de una gaussiana, que se ilustra en la figura 
Fig. 4.10. Además podemos considerar que una gaussiana sería la versión 
suavizada mediante líneas curvas de una barrera de potencial cuadrada, y en 
ese sentido quizás se pueda interpretar como más parecida al caso real. 







































































0 10 20 30 40 50 60 
x(nm) 


Figura 4.10: Barrera de potencial gaussiana (línea discontinua) cuyos paráme- 
tros son Uy =0.5 eV, a= 10 nm, zo = 30 nm. Hemos aproximado la barrera de 
potencial por una escalera (línea continua) con 20 escalones, para representar 
de forma visual la aproximación realizada. 


En particular, nuestra barrera de potencial gaussiana viene caracterizada 
por la expresión: 
U(x) = Us : exp E (4.38) 
donde Uy es la altura máxima de la barrera de potencial, xy es la abscisa 
en la cual está centrada, y a está relacionada con la anchura de la campana 
mediante la expresión FWHM = 242 1n 2 a, donde FWHM es la anchura total 
a la mitad del pico máximo. En nuestro ejemplo Uy = 0.5 eV, xy = 30 nm y a 
= 10 nm, además los portadores inciden sobre dicha barrera con un ángulo de 
10. Aproximamos la gaussiana por una escalera de 300 escalones, cuyo valor 
en cada escalón coincide con la semisuma en los extremos de dicho intervalo 
de la función gaussiana. A continuación representamos para la aproximación 
anterior el coeficiente de transmisión en la figura Fig. 4.11. En la cual obser- 
vamos un mínimo amplio cuando la energía del portador incidente coincide 
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con el valor máximo de la barrera de potencial original Uy. Para valores de 
energía inferiores a Up, tendrá un comportamiento oscilante asociado a las 
diferentes alturas de los tramos de barrera de potencial contiguos, mientras 
que para valores superiores a Uy no encontrará ningun obstáculo en la trans- 
misión, y alcanzará su valor máximo T'=1, manteniéndose indefinidamente 
en él. 























0.0 0.2 0.4 0.6 0.8 1.0 
E(eV) 


Figura 4.11: Coeficiente de transmisión frente a la energía del portador inci- 
dente para una barrera de potencial gaussiana caracterizada por la ecuación 
Eq[4.38], cuyos parámetros son Uy = 0.5 eV, a= 10 nm, zo = 30 nm y el 
ángulo de incidencia de los portadores es 4=10". Hemos aproximado la barre- 
ra de potencial por una escalera con 300 peldaños. Observamos un mínimo 
amplio cuando E coincide con Up. Para valores de energía inferiores a Up 
tendrá un comportamiento oscilante asociado a las diferentes alturas de las 
barreras de potencial contiguas, mientras que para valores superiores a Up 
no encontrará ningún obstáculo en la transmisión y tomará el valor máximo 
I=l, 


4.4. Resumen y conclusiones 


Después del estudio de dos barreras de potencial cuadradas aplicadas 
a una lámina de grafeno en un capítulo anterior, pretendemos seguir con 
el desarrollo ampliándolo a múltiples barreras de potencial. Éstas pueden 
realizarse experimentalmente mediante un dispositivo con múltiples puertas 
metálicas, mientras que la barrera de potencial escalonada puede utilizarse 
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para aproximar cualquier potencial que varíe suavemente. 


Para múltiples barreras de potencial con la misma altura, un aumento en 
el número de barreras no afecta a la anchura del pseudogap que aparece en la 
representación del coeficiente de transmisión frente a la energía, sin embargo 
aparecen resonancias muy selectivas en energía en el interior del pseudogap, 
que quizás pudieran utilizarse en la realización de filtros muy selectivos. 


En el caso de múltiples barreras de potencial de diferente altura, obser- 
vamos que cuando el ordenamiento desde la izquierda y la derecha de las 
barreras coincide, es decir tenemos simetría especular, el resultado para el 
coeficiente de transmisión es el mismo. También que al representar el coe- 
ficiente de transmisión frente a la energía, aparece un mínimo asociado a 
cada barrera de potencial, el cual puede llegar a convertirse, dependiendo de 
los parámetros utilizados, en un pseudogap de energía. El pseudogap puede 
ajustarse mediante las siguientes variables: ángulo de incidencia, energía po- 
tencial de las barreras y número de barreras de potencial. Dependiendo de los 
parámetros utilizados, el extremo superior de un pseudogap puede alcanzar 
o superar el extremo inferior del siguiente pseudogap en energía, y así llegar 
a colapsar en un único pseudogap. Para poder llegar a ajustar independiente- 
mente el extremo inferior y el extremo superior en energía de un pseudogap, 
necesitaríamos al menos dos barreras de potencial diferentes. La posibilidad 
de crear zonas de transmisión y pseudogaps a partir de voltajes externos aso- 
ciados a las barreras de potencial, ofrece la posibilidad de construir un filtro 
dependiente de la energía. 


La barrera de potencial escalonada brinda la oportunidad de simular cual- 
quier potencial que varíe suavemente mediante una aproximación por tramos 
escalonados. Hemos considerado interesante aplicarlo a una barrera gaus- 
siana, ya que es una versión redondeada, y por tanto quizás más real de 
una barrera cuadrada. Para dicho caso, en el coeficiente de transmisión no 
aparecen mínimos asociados a cada barrera de potencial, aunque sí aparece 
claramente un mínimo en energía asociado al valor máximo del conjunto de 
las barreras de potencial. A partir de dicho valor de energía, la transmisión 
es máxima T'=1. 
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Capítulo 5 


Barreras de potencial bajo 
campo eléctrico 


5.1. Introducción 


Continuando con nuestro estudio de transmisión [71] y transporte de co- 
rriente [86] en dispositivos GFET, incorporamos ahora al modelo considerado 
en capítulos anteriores un nuevo elemento. Hasta ahora hemos considerado 
al dispositivo como no polarizado, con sus portadores en equilibrio electro- 
químico con los contactos a un potencial fijo único. Existe gran interés en el 
efecto que produce un campo eléctrico sobre el grafeno [94], de modo que, 
persiguiendo este objetivo, vamos a introducir ahora una polarización del 
GFET aplicando un potencial débil entre drenador y fuente, lo que introdu- 
ce en la lámina de grafeno un campo eléctrico uniforme cuyo efecto sobre el 
transporte electrónico queremos considerar. Tratamos dicho campo por me- 
dio de un potencial adicional no constante, al que se encuentran sometidos 
los portadores de corriente, que varía linealmente con la posición entre fuente 
y drenador. Se sabe que las soluciones de la ecuación de Schródinger para 
un potencial lineal pueden construirse en términos de las funciones de Airy. 
Esta aproximación se utiliza en [95] para caracterizar la unión P-N de un 
nanoribbon de grafeno. Sin embargo, en nuestro caso este tratamiento no 
resulta apropiado, dado que los electrones en el grafeno no vienen descritos 
por una ecuación de segundo orden, como es la de Schródinger, sino por una 
ecuación de Dirac de primer orden cuyas soluciones son diferentes. 

Las soluciones de la ecuación de Dirac en un potencial lineal han sido 
consideradas directamente [96,97] o por medio de una transformación canóni- 
ca [98,99] que la convierte en una ecuación de Schródinger con un potencial 
de oscilador armónico. El sistema ha sido resuelto para un potencial de diente 
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de sierra [100], emergente de una estructura con múltiples pozo/barrera so- 
bre grafeno y para una barrera trapezoidal sobre grafeno [101], aproximando 
las zonas lineales por escalones. Soluciones analíticas exactas son muy esca- 
sas. Una deformación de esta barrera, que se convierte en triangular, sobre 
nanocintas de grafeno es estudiada en [102]. 


Para cumplir con nuestro objetivo, necesitamos discutir los diferentes mo- 
delos que caracterizan el transporte eletrónico en un GFET polarizado, con- 
trolado por diversas puertas de potencial. El modelo más simple consiste en 
considerar la energía potencial de los portadores constante en las diferentes 
zonas del GFET, tanto en los contactos extremos, como en zonas interme- 
dias o en las zonas donde aplicamos barreras de potencial constante mediante 
puertas [27, 103], es decir se aproxima el potencial lineal por diferentes ba- 
rreras de potencial constante, que promedian de algún modo el perfil del 
potencial existente. Es posible también [104] aproximar el potencial lineal 
por una serie de tramos rectos consecutivos de valor constante, mediante un 
procedimiento análogo al que hemos utilizado para la barrera gaussiana. Un 
modelo intermedio consiste en considerar tramos de energía potencial cons- 
tante en los extremos y que el campo eléctrico únicamente afecte a las zonas 
bajo la influencia de las barreras de potencial, como por ejemplo sucede para 
semiconductores tradicionales en [105]. 


En la figura Fig. 5.1 representamos nuestro modelo de GFET, con la 
lámina de grafeno encapsulada entre delgadas láminas de nitruro de boro he- 
xagonal (h-BN), que es un extraordinario dieléctrico que mantiene al grafeno 
plano y aislado del entorno. Inicialmente Dean lo utilizó para que al deposi- 
tar el grafeno sobre h-BN fuese plano. La rugosidad superficial del h-BN es 
del orden de 1/3 de la rugosidad del SiO, y aproximadamente toma el valor 
de 0.33 nm (106). Los dispositivos de grafeno sobre h-BN presentan respecto 
a los dispositivos con sustrato SiO una mejor movilidad, una reducción en 
la inhomogeneidad de portadores, y una reducción en el dopado intrínseco. 
Posteriormente, también se puso h-BN sobre grafeno para aislarlo del exte- 
rior. El h-BN no afecta a los potenciales que pretendemos aplicar al grafeno. 


Los portadores, cuya concentración en la lámina de grafeno se contro- 
la por la back gate, son impulsados por la diferencia entre los potenciales 
químicos de fuente y drenador us — Up = eVps que es proporcional a la ten- 
sión de polarización del GFET y, por tanto, al campo eléctrico longitudinal 
existente en la lámina. Los portadores se encuentran a potenciales fijos Vs y 
Vp a la entrada y salida del dispositivo, y tienen que atravesar una barrera 
de potencial determinada por las puertas top gate y back gate, a lo largo 
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Top gate 


Back gate 





Figura 5.1: Modelo del GFET donde el grafeno se encuentra encapsulado 
entre capas de h-BN, lo cual le aporta que sea plano, no corrugado, y le aisla 
del exterior. 


de cuya región de influencia el potencial eléctrico cae Vps. Es el efecto de 
este campo eléctrico longitudinal sobre el coeficiente de transmisión de los 
portadores en el dispositivo lo que vamos a estudiar en este capítulo. Para 
ello comenzaremos considerando campos eléctricos débiles, y posteriormen- 
te consideraremos una aproximación a lo que podría suceder para campos 
eléctricos elevados. 


5.2. FElectrón bidimensional de Dirac en un 
potencial lineal 


Buscamos la función de onda exacta para describir a un electrón bidi- 
mensional de Dirac sometido a un potencial que varía linealmente con su 
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Vos 





Figura 5.2: Vista superior del GFET, en la cual visualizamos el campo eléctri- 
co mediante flechas verdes entre drenador (D) y surtidor (S). La lámina de 
grafeno (azul) está sometida a una barrera de potencial Vr aplicada me- 
diante una puerta superior (naranja). 


posición. La energía potencial de dicho electrón se expresa como: 


U(x) = —eV (2) = —e(V, + Va 1) = W, - Waz (8.1) 


que deben ser solución del hamiltoniano: 


H = Up(0gPg + 09Py) + (W¡ - Wa ola (5.2) 


siendo 0, y 0y las correspondientes matrices de Dirac, e l¿ la matriz 
identidad. 

Cada uno de estos electrones del grafeno con energía E y pseudomomento 
próximo a los puntos de Dirac (k.,, k,), se describe mediante una función de 
onda espinorial de dos componentes 


¿yy 0) pi(x) 
( palz) ) 
cuy E0 (hy ¿ED p(a) +9 (2) ) 63 


V (zx, y, t) 


en la cual la función p(+) debe ser solución de la ecuación diferencial de 
segundo orden: 
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o"(x) 4 [€ - a =ki+ ¡A la) =0 (5.4) 


hur hup 
que al aplicarlo al potencial lineal que tenemos en la expresión Eq][5.1], 
obtenemos: 


o" (a) y (E — E sea = (13 + 7) )eta) =0 (5.5) 


Esta ecuación Eq[5.5] puede ser interpretada matemáticamente como la 
ecuación de Schródinger para un potencial cuadrático de un oscilador [99]. 
Como novedad, nosotros aplicaremos esta interpretación para resolver el po- 
tencial lineal en nuestro dispositivo de grafeno. Obtenemos así nuestra solu- 
ción exacta para el sistema de Dirac que estamos considerando, para lo cual 
introducimos 





ef ELO 
ela) =p (n(2)) (5.6) 
donde la nueva coordenada n(x) es 
da E-— W,; + Wo) x 
vV hupW,) 


F (1) se identifica como una solución de la ecuación diferencial de Hermite 
respecto a la variable 7: 


F"(n) — 2nF'(n) + 2iAF(n) =0 (5.8) 


ne) = (=1) 





(5.7) 


con parámetro característico 


z kjhu F 
2 Wa 
La solución general a esta ecuación puede construirse a partir de diversas 
funciones, hipergeométricas, de Hermite, de Weber, etc. que están todas rela- 
cionadas. Nosotros elegimos para construir la solución, la función de Hermite 


H[—1A, mn] y la función hipergeométrica ¡F (12, 2,177] como: 





(5.9) 


. A 1 
F(n) =C1 H|[-1A, nm] +0 1F1 (5, | (5.10) 


De este modo, el estado de los electrones del grafeno sometido a un poten- 
cial lineal puede expresarse como combinación lineal de dos pseudoespinores 


V(x, y, t) =C;i Wi (x,y,t) + Ca Wa(x, y, t) (5.11) 


101 


Potencial bajo campo eléctrico 





siendo 








-1)1/4 hw 
Hino] ALTE py HL m(2)] 
layt)=AS LO (5.12) 
H|[—iAn(a))]+ JW kyH [-4A-1,n(x)] 


7 = 1/4 /ñ : 
141 [+3 4007 IL po) Fr [+13 Eto) 
yv Wa 222, 
V)(z, Y, t) E A(z, Y, t) 











sá [+3 2a07 O Ena [1+:3, 207 
(5.13) 
y 
] ¡EW y Wo 
A(z, Y, t) == uy RH un j Ms sd (5.14) 


5.3. Una barrera de potencial sometida a cam- 
po eléctrico 


El estudio del potencial lineal no basta para caracterizar un GFET pola- 
rizado. La presencia de la top gate obliga a añadir una barrera de potencial 
constante, que superpuesta a la energía eléctrica establecida por la tensión de 
polarización Vps, proporcione la energía potencial total de los portadores de 
corriente al atravesar el GFET. Consideramos que la energía potencial entre 
los contactos extremos (surtidor y drenador) cae principalmente en la parte 
del grafeno que se encuentre debajo de la puerta metálica superior a la cual 
aplicamos una energía potencial constante U¡. Por tanto, los portadores en 
el grafeno conectado a cualquiera de los dos extremos, surtidor y drenador, 
tendrán aproximadamente la energía correspondiente a éstos Us y Up res- 
pectivamente. Este modelo se justifica a partir de mediciones experimentales 
que demuestran que el grafeno situado bajo una puerta metálica superior es 
más resistivo que el resto de grafeno que no está bajo la influencia de ninguna 
puerta superior. 


Así, para construir la energía potencial que influye en los portadores, 
hemos de considerar tres tramos. Las zonas extremas están caracterizadas 
por una energía potencial constante Us y Up, asociadas respectivamente a 
los terminales de surtidor y drenador, y bajo la puerta, se ha de añadir una 
barrera de potencial de altura U¡ y anchura de puerta lyate = 2 — 21, a la 
energía potencial lineal asociada al campo eléctrico. 
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Á U, 

Us 
v X1 Xx | Up Xx 








> 


Figura 5.3: Vista superior del dispositivo (a). Energía potencial U frente a la 
distancia cuando aplicamos una barrera de potencial U, dentro de un campo 
eléctrico asociado a la diferencia de potencial entre drenador y surtidor (b). 


De este modo, la energía potencial a la cual está sometido el electrón, 
representada en la Fig. 5.3, se expresa como: 


Us si z<o 
U(x)=3 —Fle—-xw1)+U54+ 01 si 11 << 22 (5.15) 
Up si Za<z 


siendo /" la fuerza asociada al campo eléctrico (en módulo) que aparece entre 
surtidor y drenador 


= -V 
pe (5.16) 


Lo — Ei Ta — L1 





En cuanto a los casos extremos para el campo eléctrico (en módulo), 
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haremos la siguiente distinción: 


F lgate <<U1 —> campo eléctrico muy débil (5.17) 
F -lgate >> U1 —> campo eléctrico muy elevado (5.18) 


e inicialmente estudiaremos el campo eléctrico bajo, que además es el que 
tiene más sentido tanto desde el punto de vista teórico como experimental. 

Los valores de pico en los extremos de la barrera de potencial, según la 
expresion Eq[5.15] son: 


Ups = U(x1) = Us ul Us (5.19) 
Unearka = U(z>) = Up +— Ur (5.20) 





La resolución de la ecuación de Dirac para este sistema, requiere conside- 
rar estados distintos para el electrón en las 3 zonas, como en casos anteriores. 
Este nuevo potencial no cambia la simetría del hamiltoniano y siguen con- 
servándose, como en capítulos anteriores, la energía E y el momento k,,. 
Llamamos ky = (kz, ky) al momento del electrón en la zona extrema l, y 
krrr = (q,,k,y) al momento del electrón en la zona extrema III. Ambos se 
encuentran relacionados a través de la energía, como observamos en las ex- 
presiones Eq[5.21-5.23]: 





| E—Ug| 
A .21 
k no cos Q (5.21) 
| E—Ug| 
SW, E RA .22 
19 e sen Q (5.22) 
de = Le (E -—Up)?-—(E-—Us)? sinó (5.23) 
hup 


De acuerdo con estas consideraciones, podemos construir la solución a 
nuestro potencial como: 








rr ta ada) File] Fs[] 
Vir(2,y,t) = E e 
nx, y, t) 2LW ; ( F,[x] ) ] ( Falx] ) 
e cb, 1 4 
? t ES AS 1Qx Y 5.2 
AA Conte) dd 
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siendo 
A | 1/4 ED)Vhop | 1/4 
Ple] = Hi», (1) [x] E NS E Aza), 


1): V hop 0% 
Folx] = H[-4A, (Dx 4 E ati, /42[2]), 























A v hu o E A 
Fslx] = 1 (5, 3 ica]? + ma tela] ¡F [1 OL —izla]? A 
e vV hu 0 A 
Fall =5ñA [5 3 —izla]? =- Am dla] ¡F [1 O —izla]?| j 
Es k?hv 
da es E AA (5.25) 
hupWoa 2W, 
Los parámetros auxiliares de las funciones anteriores son: 
Us — U, 
W, = Fx+U,+0s, W=F=-2 2 (5.26) 
Lo — T1 


donde hemos denotado sy = sign(E — Us) y sp = sign(E — Up) en las 
zonas extremas. Además hemos considerado que las partículas se propagan 
inicialmente de izquierda a derecha, pudiendo reflejarse o transmitirse en el 
límite entre zonas, pero en la región III, únicamente quedan partículas pro- 
pagándose hacia la derecha. 


También ahora se debe verificar la condición de contorno de que la función 
de onda debe ser continua en la frontera entre zonas, por lo que debe cumplir 
la igualdad en los dos límites x= 11 y 1= 22, 


W ¡(01 y,t) = ViJ(21, y, Pa Vi¡(x9, y,t) > Viri(22, Y, tl (5.27) 


Este procedimiento estándar nos proporciona el sistema de ecuaciones a re- 
solver. Cuatro ecuaciones con cuatro incógnitas que nos dan los coeficientes 
r,C1, C2 y t que omitimos. Una vez obtenido el coeficiente t de la tercera zona, 
buscamos el coeficiente de transmisión 7”. Teniendo en cuenta que ahora, en 
presencia de un campo eléctrico, las regiones de energía potencial constante 
en las zonas extremas toman diferente valor de modo que debe construirse 


como T = sssp E/2 | £ |2. Así encontramos 
lok 





105 


Potencial bajo campo eléctrico 





N 





Ts 7 (5.28) 
N =sssp16kqk, q, sinh pa Al (5.29) 
D = Abs[(sp q+Filx2] — qFa[co] ss k-Fs[x1] + kFalxs)) 

(ss k-Fi [1] + kPola1]))(ss q+Fs[xe2] — qFa[2217, (5.30) 


que, junto a las expresiones Eq[5.25] expresa el coeficiente de transmisión 
de los electrones (directo e inverso, pues ambos son iguales como se puede 
comprobar por cálculo directo) en términos de los parámetros que gobiernan 
el sistema E, Ug, Up, U,, 11,12 y P. Analizamos ahora gráficamente esta de- 
pendencia. Salvo que se indique lo contrario, el ángulo de incidencia de los 
portadores es de 10%: 


= Dependencia del potencial fuente-drenador. 


En la figura Fig. 5.4 representamos el coeficientes de transmisión frente 
a la energía para una barrera de potencial de 0.4 V y 100 nm para dis- 
tintos valores del campo eléctrico : 0, 0.5, 1 y 2 mV/nm asociados a Us 
= 0 y Up =0, - 50, -100, -200 meV, respectivamente . Así, observamos 
el efecto sobre T del aumento del campo eléctrico asociado a un au- 
mento en Ug— Up. A campo eléctrico nulo observamos pendientes muy 
verticales en las paredes laterales con muchas y amplias oscilaciones en 
T. Un aumento del campo eléctrico aumenta la anchura del pseudogap 
energético, atenúa la caída del pseudogap de energía, y disminuye la 
amplitud de las oscilaciones. 


= Dependencia de la anchura de la puerta superior. 


En la figura Fig. 5.5 se representa el coeficiente de transmisión de una 
barrera de potencial de 0.4 V con Us = 0 meV y Up =- 50 meV. Para 
anchuras de 100 nm (azul), 50 nm (verde) y 25 nm (rojo) se generan 
campos eléctricos de 0.5, 1 y 2 mV/nm. Observamos en dicha figura que 
cuánto más ancha es la barrera de potencial más profundo y definido 
es el mínimo, pero la anchura del pseudogap en energía básicamente se 
mantiene constante. 


Además al comparar la gráficas Fig. 5.4 y Fig. 5.5, con los mismos 
valores de campo eléctrico, se observa que un aumento del campo debido 
a un aumento de la energía de polarización produce un aumento en la 
anchura del pseudogap energético y una dismunución del número de 
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Figura 5.4: Coeficientes de transmisión para una barrera de potencial de 0.4 V 
y 100 nm de anchura dentro de un campo eléctrico nulo (lila), de 0.5 (azul), 
1 (verde) y 2 (rojo) mV/nm. El aumento del campo eléctrico provoca un 
aumento en la anchura del pseudogap, así como pendientes laterales menos 
verticales y oscilaciones de menor amplitud. 


oscilaciones en T', mientras que una disminución de la anchura de la 
barrera no modifica en la práctica la anchura del pseudogap y aumentan 
las oscilaciones. 


= Dependencia del potencial de puerta superior. 


En la figura Fig. 5.6 tenemos dentro de un campo eléctrico de 1 mV /nm 
una barrera de potencial de anchura 100 nm, y con distintas alturas 
0.2 (rojo), 0.4 (verde) y 0.6 (azul) eV. Al representar el coeficiente de 
transmisión frente a la energía, observamos que a medida que aumenta 
la altura de la barrera de potencial aplicada, el pseudogap energético 
se va desplazando hacia energías superiores en torno a dichos valores, a 
la vez que se ensancha y define más, es decir con bordes más verticales. 
El desplazamiento y ensanchamiento del pseudogap para barreras de 
potencial cada vez mayores es semejante a lo que habíamos obtenido 
para una barrea de potencial sin campo eléctrico aplicado. 


= Dependencia del ángulo de incidencia de los portadores. 


En la figura Fig. 5.7 tenemos dentro de un campo eléctrico de 1 mV /nm 
una barrera de potencial de 0.5 V y 100 nm, con distintos ángulos de 
incidencia del portador 1%, 5% y 10. A medida que aumenta el ángulo 
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Figura 5.5: Coeficientes de transmisión para una barrera de potencial de 0.4 
V, y anchuras 100 (rojo), 50 (verde) y 25 (azul) nm siendo Us = 0.2 eV y Up 
= O eV. Cuánto más ancha es la barrera de potencial más profundo y definido 
es el mínimo, pero la anchura del pseudogap se mantiene aproximadamente 
constante. 
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Figura 5.6: Coeficientes de transmisión para una barrera de potencial de 100 
nm en un campo eléctrico de 1 mV /nm para barreras de potencial de alturas 
0.2 (rojo), 0.4 (verde) y 0.6 (azul) eV. Al aumentar la altura de la barrera 
de potencial, el pseudogap se va desplazando hacia energías superiores, a la 
vez que se ensancha y se hace más vertical. 
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Figura 5.7: Coeficientes de transmisión para una barrera de potencial de 0.5 
V y 100 nm dentro de un campo eléctrico de 1 mV/nm para ángulos de 
incidencia de los portadores ¿ de 1? (rojo), 5? (verde) y 10? (azul). A medida 
que aumenta el ángulo de incidencia, los mínimos son cada vez más profundos 
y más anchos. 

















Figura 5.8: Coeficiente de transmisión para una barrera de potencial de 0.5 
V, y anchura 10 nm con campo eléctrico de 2mV/nm (rojo discontinuo), sin 
campo eléctrico (azul discontinuo), y anchura 100 nm con campo eléctrico de 
2 mV /nm (lila), sin campo eléctrico (añil). 
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de incidencia, los mínimos son cada vez más profundos y definidos, 
es decir se ensanchan. Sabemos que para incidencia perpendicular, es 
decir f = 0? el coeficiente de transmisión debe ser 1, lo cual enlaza muy 
bien con la forma de la gráfica obtenida para fp = 1”. 


Podemos interpretar de todo lo que se ha expuesto hasta ahora, que es la 
altura de la barrera el factor que más influye en el coeficiente de transmisión 
de una barrera de potencial sometida a campo eléctrico y éste tiene un papel 
menor. Si se compara el comportamiento de barreras con o sin campo eléctri- 
co como en la Fig. 5.8, se obtiene que, a medida que aumenta la anchura, los 
resultado con campo y sin él se van diferenciando cada vez más. Es decir, el 
campo eléctrico tendrá una mayor influencia cuánto mayor sea la anchura 
de la barrera de potencial. 


5.4. Dos barreras de potencial sometidas a 
campo eléctrico 


Consideramos a continuación el caso de dos barreras de potencial Uy y Uz 
genéricas, en principio distintas tanto en altura como en anchura, sometidas 
a un campo eléctrico debido a la diferencia de potencial entre drenador y 
surtidor. Al añadir ahora una barrera de potencial con respecto al caso ante- 
rior, aparecen dos nuevas zonas. Caracterizamos la energía potencial en este 
caso por medio de: 





Us si TI<LI 
—Fl(x—x)+U5+U, si 11<x<t>z 

Ult = <> Uso si Z2<1I<LI3 (5.31) 
—F (1-11) +Up+U32 si 13<U< IA 
Up si Ta < 


siendo F la fuerza asociada al campo eléctrico (en módulo) que se esta- 
blece entre surtidor y drenador, y que únicamente se presenta en la zona de 
grafeno bajo las puertas superiores, que introducen sendas barreras de po- 
tencial de alturas U¡ y Uz. En el modelo utilizado se considera que el campo 
eléctrico etablecido por la polarización del GFET es el mismo bajo las dos 
puertas 


F= Us Up =e- Ys Vo (5.32) 
(12 — 11) + (xa — 23) (12 — 11) + (24 — 23) 
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Figura 5.9: Vista superior del dispositivo (a). Energía potencial U frente a la 
distancia cuando tenemos dos barreras de potencial U, y Uz sometidas a un 
campo eléctrico (b). 


El valor de U;,n' asociado a la energía potencial constante que aparece en 
la zona intermedia entre las dos barreras, tal y como aparece representado 
en la figura Fig. 5.9, resulta: 


Us(ta4 — 23) + Uplx2 — 11) 


Oia = 
j (13 — 21) + (24 — L3) 


(5.33) 





de modo que siempre se cumple Us > Uin; > Up. 


De la figura Fig. 5.9 observamos que existen 3 zonas (las zonas extremas y 
la intermedia) en las cuales la energía potencial es constante, y otras 2 zonas 
(donde aplicamos las barreras de potencial) en las cuales la energía potencial 
varía linealmente con la misma pendiente en ambas, pero a diferentes alturas. 
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Las funciones de onda en cada zona son: 


$ e Ukyy= Ft) " 1 E 1 
Lt == 4 004 + pe ter E A 
0 ens) Corts ) 
(kyy-Ft+ A 4 y 2?) 
e Up a F, [2] F3|x] 
y , 56 == C + C: > 
n(x y ) 2LW 1 ( Faja] ) 2 ( Fla] 


e Ukyy= Ft) . 1 —ikgz 1 
ente Gp tas Jl) y 














¿(kyy Eb, O as Wa 2?) 
e h hu p 2hv p Gi; [2] G3 [x] 
y aos a C, , 
Iv(z, Y, ) 92LW 5 ( Go») [2] ) A E ( Ga[x] 
y t) a 1quT 1 (5 34) 


Para definir las dos nuevas zonas hemos tenido que introducir las siguien- 
tes funciones adicionales a las que se introdujeron en las expresiones Eq[5.25] 


(—D)4yRop 
VW 
(1) 44h 
Wa 


Gil] = Hit (1) yla] - ky HA 1, (1) [2] 


Golx H|[-iA, (1D) ylx]] y 


kE [A 1, (1) ye) 























Gale] = 1F; (5, 5 —iy[x] | + a Fue 1 [1 O —iy[x] | 
Al... vVhoup ¿TAS ES 
Gale] = 1F [3 > —iy[x] | =- a Fu) 1F; [1 a —iy[x] | 


(5.35) 


siendo los parámetros 


A E 
- Un Ñ pza E h2, (5.36) 
UF 





E-— 3 + Waz . 
E in = sig E — Us), (5.37 
ula 7, Sine =sign(E— Um), (537) 
Wi = Fo +0 + Us Wi Pr (68) 


De nuevo igualamos las funciones de onda en las fronteras entre zonas 
contiguas para obtener el coeficiente t en la zona V. 


Vi(x21,y,t) = Vyr(u1, y, t), Vir(x2,y,t) = Vir (xo, y,t), (5.39) 
Vin(23,y,t) = Viy (23, y, t), Viy(za,y,t) = Vylza, y,t). (5.40) 
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Y volvemos a calcular el coeficiente de transmisión Y” = sssp 2 Abs[t]?, 
obteniéndose explícitamente 


s558p 256 k q? ko qu Abs[x,]? sinh*[A3] 


Abs[eirolr3=22) (Ry = Ra)(Rz = Ra) = eino(u322) (Sy = S2)(S3 = Sa)]? 
(5.41) 





TT = 


siendo y, = y [xi], 21 = 2 [21] y 

















Ri = (ssk_Fj[21] + kPolza1D) (Sins 5-P3[22] + 5F1122)), 
Ra = (ssk_Fsla] + kFalar]) Sins 5-Fi[22] + 543[20)), 
Rg = (Sin k+Gilys] — kGalya))(sp a+Gslya] — aGalya)), 
Ra = EA + G3ly3 3 Ga [y3])(sp q+61 [Ya] - qG) [ya)), 
Si = (ssk-Fr [2] + kPola)) sin: 5+F3[z0] — Fa[z2)), 
S2 = (ssk-PE3l[2] + kErla1))sin: 5+F1[20] — kFo[z2)), 
S3 = (e Gi [yz + KG) yal) (sp q, -Galya] NN aGaly AD 
Sa =  (Sini-Galys] + «Galyal)(sp q+Gilya] — qGalyal). (5.42) 








Aparecen en este sistema nuevos efectos, como son: 


= Debido a la pérdida de simetría del sistema, la energía potencial total 
depende de la posición relativa entre sí, que ocupen las barreras de 
potencial (primera o segunda), lo cual está asociado al valor máximo 
de energía potencial que puedan alcanzar cada una de estas barreras 
de potencial UPew;, o bien Dent: Ya que para el pico más alto peak1, la 
barrera de potencial en primera posición se sumará con Ug, mientras 
que si aparece en segunda posición se sumará con Ui,+, teniendo en 
cuenta que el potencial en la zona intermedia U;, siempre es inferior 
al potencial en la parte izquierda Us, ya que hemos tomado Uy > Up. 





Uler!, =U (21) =Us+Ún, 
Da = U (12) = Uint + Un, 
bar2 (5.43) 
Uneak1 = U(x3) = Uint TT Oz, 
Upeaza = U (24) = Up + Uo. 


= La segunda barrera de potencial puede llegar a quedar enmascarada, es 
decir que su efecto puede llegar a no ser detectable. Eso sucederá cuando 
se verifique Eq[5.44]: 
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Figura 5.10: Representación de energía en dos barreras de potencial, cuando 
la segunda queda enmascarada al verificarse Eq[5.44]. 

















dls al el 


0.0 0.2 0.4 0.6 0.8 1.0 
E(eV) 





Figura 5.11: Coeficientes de transmisión para una barrera de potencial de 
0.4 eV y 75 nm separada 50 nm de la segunda. La segunda tiene una altura 
de 0.1 eV (rojo) o 0.9 eV (azul) con 25 nm de anchura. Campo eléctrico 2 
mV/nm y ángulo de incidencia 10?. Es patente el efecto de enmascaramiento 
en la gráfica roja al verificarse Eq[5.44]. 


En la figura Fig. 5.10 representamos el diagrama de energía para dos 
barreras de potencial, siendo el valor de Uz lo suficientemente pequeño 
como para que se verfique la expresión Eq[5.44], y por tanto esta se- 
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gunda barrera de potencial quede enmascarada. Este efecto se ilustra 
en la Fig. 5.11 donde se representa el coeficiente de transmisión frente 
a la energía del portador incidente para 2 barreras de potencial en un 
campo eléctrico de 2 mV/nm y ángulo de incidencia 10. La primera 
barrera de potencial, de 0.4 eV y 75 nm, está separada 50 nm de la 
segunda, que tiene 25 nm de anchura, y una altura de 0.1 eV (rojo) 
o 0.9 eV (azul). Con estos parámetros tenemos que, para Up = - 0.2 
eV, es Uin¿= - 0.15 eV, y Uin: + U2= - 0.05 eV en el primer caso y 0.75 
eV en el segundo. Con estas condiciones, en el primer caso la segunda 
barrera está enmascarada. Debe esperarse que el coeficiente de trans- 
misión presente un primer mínimo en ambos casos, tras lo cual, en el 
segundo caso (azul) vuelva a aparecer una segunda zona de extinción 
correspondiente a la segunda barrera que no existe en el primero. Este 
comportaiento es precisamente lo que aparece en la figura Fig. 5.11. 


I TI TI IV v 


























Figura 5.12: Dos barreras de potencial efectivo igual: alcanzan los mismos 
z barl __ prbar2 barl  _ bar2 . 
valores de pico Uneax1 = Upeak1 Y Upeaka = Upeako y son de la misma anchura 


L¡ = L3 = L. Son idénticas por encima de Us = 0 eV. 


= Efecto de la separación entre dos barreras de potencial sometidas a 
campo eléctrico. 


Denomino barreras de potencial efectivo igual, a las que alcanzan los 


a - barl __ prbar2 _ — Trbarn barl  __ prbar2 
mismos valores de pico Ur = Upeaja =*** = Upeater Y Upeaxa = Upeaico 


Pp 
darxa, Para diferentes alturas de las barreras de potencial 
U;¡, Uz, --- U,,, estando sometidas al mismo campo eléctrico, y además 


tienen la misma anchura L, = La =-+++ = L, = L. Estas se comportan 
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como barreras idénticas por encima de U, = 0 eV, como representamos 
en la figura Fig. 5.12. 


En la figura Fig. 5.13 representamos frente a la energía el coeficiente 
de transmisión de dos barreras con igualdad efectiva, siendo U¡ = 0.6 
eV y Uz = 0.8 eV dentro de un campo eléctrico de 4 meV/nm. (Puesto 
que Ug =0€eV y Uin: =- 0.2 eV, ambas barreras de potencial alcanzan 
los mismos valores de pico). Ambas son de la misma anchura de 50 
nm, separadas 1 nm (a) o 50 nm (b), y el ángulo de incidencia 10". 
Obsevamos que únicamente existe una zona donde el coeficiente de 
transmisión frente a la energía del portador incidente es mínimo, y que 
el número de picos aumenta a medida que aumenta la distancia de 
separación entre las dos barreras de potencial. 


Para dos barreras de potencial con igualdad efectiva, existen valores 
particulares de los parámetros para los cuales aparecen múltiples reso- 
nancias claramente definidas en el pseudogap energético que presenta 
el coeficiente de transmisión, de forma similar a las resonancias que 
habíamos obtenido para dos barreras idénticas (sin campo eléctrico 
aplicado). Para visualizarlo, representamos en la figura Fig. 5.14 el 
coeficiente de transmisión para dos barreras de potencial que cumplen 
esta condición, siendo U¡ = 0.3 eV y Uz = 0.5 eV, que debido al campo 
eléctrico dentro del cual están inmersas de 4 meV/nm, ambas alcanzan 


a . / > bart __ prbar2 _ 
los mismos valores de pico en energía potencial Ucax1 = Upeak1= 0.3 eV 
y Ubarno = Ucaja = 0.1 eV según la expresón Eq[5.43]. La anchuras de 


ambas es de 50 nm, y están separadas entre sí 1 nm, siendo el ángulo de 
incidencia 10%. Si consideramos la forma de la energía potencial frente a 
la distancia, por encima de O eV, ambas son completamente idénticas. 


Campo eléctrico elevado 


En la Sección 5.3, estudiamos la barrera de potencial sometida a un cam- 
po eléctrico débil, asociado a la diferencia de potencial entre drenador y 
surtidor. A continuación haremos lo propio para un campo eléctrico fuerte. 
Somos conscientes de que experimentalmente habría que tratarlo con sumo 
cuidado, debido a que hay un conjunto de fenómenos que no se han conside- 
rado, como el efecto de no-linealidad en la banda de dispersión para valores 
energéticos elevados, ni la posible ruptura del dieléctrico debido a un campo 
eléctrico elevado. Sin ambargo, se continúa con el desarrollo teórico, por si 
acaso pudiese resultar de interés, al ser una aproximación a dicho caso. 
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Figura 5.13: Coeficientes de transmisión para dos barreras de potencial efec- 
tivo igual, siendo U¡ = 0.6 eV y Uz = 0.8 eV en un campo eléctrico de 4 
meV /nm con 50 nm de anchura cada una y separadas 1 nm (a) o 50 nm (b) 
y ángulo de incidencia 10. Al alcanzar los mismos valores de pico y tener la 
misma anchura, dan lugar a un único pseudogap. Al aumentar las distancia 
de separación entre las barreras de potencial aumenta el número de picos que 


aparecen. 
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Figura 5.14: Coeficiente de transmisión para dos barreras de potencial efec- 
tivo igual, siendo U¡ = 0.3 eV y Uz = 0.5 eV, en un campo eléctrico de 4 
meV/nm. La anchura de cada una es 50 nm y separación 1 nm, siendo el 
ángulo de incidencia 10%. Aparecen múltiples resonancias en el pseudogap 
energético. 


Cuando se verifique la condición Eql[5.45], nuestra barrera rectangular 
asimétrica, con la parte superior lineal, puede llegar a convertirse de manera 
efectiva en una barrera triangular [102], ya que únicamente tendrá sentido 
representar energías superiores a Us. En la figura Fig. 5.15 representamos 
tres barreras de potencial U, > U2z > U3z de tal forma que se verifica la 
condición Eq[5.45]. Cuanto menor sea la altura de la barrera de potencial, 
menor será la altura y anchura del triangulo efectivo de energía potencial 


Us > U¡+Ubp (5.45) 


Sin embargo, estas barreras triangulares suponen utilizar campos eléctri- 
cos extremadamente altos, por lo que tienen sentido teórico, pero como co- 
mentamos anteriormente poco sentido desde el punto de vista de una posible 
realización práctica. 


En la figura Fig. 5.16 siendo $ = 10” representamos el coeficiente de 
transmisión para una barrera de potencial de 100 nm en un campo eléctrico 
de 4 mV /nm, para distintas alturas de la barrera de potencial que verifican 
la condición Eq[5.45]. Así, representaremos el caso límite 0.4 eV (lila). A me- 
dida que disminuimos la altura de la barrera 0.3 (azul) , 0.2 (verde), 0.1 eV 
(rojo), observamos que el mínimo es cada vez menos profundo, como debe ser 
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Figura 5.15: Distintas barrreras de potencial U,, U2 y Uz aplicadas sobre el 
campo eléctrico asociado a Us y Up. 


porque la barrera de potencial lineal por encima de Us es cada vez de menor 
altura y de menor anchura efectiva. Se observa además, que la envolvente 
de la pendiente descendente es semejante en todas las gráficas con amplias 
oscilaciones, mientras que la pendiente ascendente comienza a distintos va- 
lores energéticos, relacionados con U; y presenta oscilaciones mucho menores. 


A continuación relacionamos nuestro potencial lineal Eq[5.15] con otro 
aproximado mediante tramos constantes Eq[5.46]: 


Us si T<X] 
Dimwlr= 4 Us 3 Ei (5.46) 
Up si To<Z 


donde Un tomará los siguientes valores posibles: el máximo Uf2? = Ug+U,, 


el intermedio Uned = (Us + Up)/2+U,, y el mínimo Un" = Up + U,, para 
ver en qué medida los valores extremos y el intermedio se ajustan al caso de 
nuestro potencial lineal Eq[5.15]. 

En las figuras Fig. 5.17 representamos el coeficiente de transmisión en las 
expresiones Eq[5.28-5.30] para Ug = 0 eV, Up = - 0.2 eV, para un campo 
eléctrico de 20 mV/nm y una barrera de potencial U; de 0.5 eV (lila), con 
una anchura de 10 nm (a) y 100 nm (b). Para utilizar la expresión aproxi- 
mada Eq[5.46] tomaremos Un = 0.5 (rojo), 0.4 (verde) y 0.3 (azul) eV y 
el coeficiente de transmisión Eq|2.16]. En todos los casos el ángulo de inci- 
dencia es de 10. Obsérvese que en la Fig. 5.17 (a) no se alcanza el mínimo 
de transmisión. Ello es debido a que la barrera es muy estrecha (10 nm). Se 
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Figura 5.16: Coeficientes de transmisión para una barrera de potencial de 
100 nm en un campo eléctrico de 4 mV/nm para distintas alturas de la 
barrera de potencial 0.1 (rojo), 0.2 (verde), 0.3 (azul), 0.4 (lila) eV. Siendo 
pd = 10”. La envolvente de la pendiente descendente es similar en todas las 
gráficas, mientras que la pendiente ascendente comienza en distintos valores 
energéticos asociados a U,. 


observa, no obstante, que de las tres aproximaciones realizadas siguiendo el 
desarrollo anterior, el caso aproximado superior 0.5 eV reproduce bastante 
bien los resultados obtenidos con el potencial lineal Eq[5.15], lo cual es lógico 
puesto que la anchura es muy pequeña. Para analizar el efecto de la anchu- 
ra, se modifica únicamente ésta, que pasa de 10 nm a 100 nm, en la figura 
Fig. 5.17 (b). Ahora se observa que el caso intermedio superior no coincide 
con el caso lineal. Las gráficas para los potenciales aproximados tienen las 
paredes laterales muy verticales, frente al caso del potencial lineal en el cual 
aparecen pendientes finitas, con oscilaciones suaves. La anchura del pseudo- 
gap energético en el caso del potencial lineal es superior a cualquiera de los 
casos asociados a los potenciales aproximados. De hecho, salvando la pen- 
diente finita, su valor mínimo en energía parece coincidir con el valor más 
bajo para UM", y su valor máximo con el valor más alto para Uf"27. Es de- 
cir, los valores energéticos extremos del pseudogap en energía para el caso 
del potencial lineal coinciden básicamente (sin tener en cuenta la pendiente 
finita) con los casos aproximados mínimo y máximo en ambos extremos. 
Cuando se verifiquen las condiciones Eq[5.47 - 5.48], nuestras barreras 
rectangulares asimétricas, con la parte superior lineal, pueden llegar a con- 
vertirse de manera efectiva en dos barreras triangulares, de forma semejante 
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Figura 5.17: Coeficiente de transmisión para Us =0€eV, Up =- 0.2V, con una 
anchura de 10 nm (a) o 100 nm (b), para un campo eléctrico de 20mV/nm y 


una barrera de potencial Uy de 0.5 eV (lila). Casos aproximados Uaproz = 0.5 
(rojo), 0.4 (verde), 0.3 (azul) eV. Angulo de incidencia 10". 
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a lo que sucedía para una barrera triangular en la figura Fig. 5.15. Puesto que 
únicamente tiene sentido en nuestro modelo representar energías superiores 
a Us, para obtener barreras triangulares deberán verificarse las siguientes 
desigualdades: 


Us + Uint 
Us + Up 


< Us (5.47) 
< Us (5.48) 

En la figura Fig. 5.18 representamos en (a) la energía potencial para dos 
barreras de potencial de 0.2 eV la primera y 0.4 eV la segunda representado 
en verde, y de 0.1 eV la primera y 0.3 eV la segunda en azul. En ambos casos 
la anchura de la barrera es de 50 nm y están separadas 20 nm entre sí, en 


un campo eléctrico de 4 meV/nm. Así, para el primer caso (verde) Uer,, 


= Oneaj= 0.2 eV, Uporio= Uñeayo = 0 eV, mientras que para el segundo 
bart __ bar2 __ barl __ bar2  __ 
caso (azul) Uca = Upeaj= 0.1 eV, Unearo= Unpeajo = - 0.1 eV. Con estos 


parámetros ambas barreras cumplen Eq[5.47 - 5.48] y se trata por tanto 
de dos barreras con potencial efectivo igual. En particular, en el segundo 
caso, el sistema inicial de barreras por encima de O eV es equivalente a un 
sistema efectivo formado por dos barreras triangulares idénticas de 0.1 eV 
de altura y 25 nm de anchura separadas 45 nm. En (b) representamos el 
coeficiente de transmisión para esa energía potencial frente a la energía del 
portador incidente con ángulo de incidencia 10”. El coeficiente de transmisión 
presenta gran número de oscilaciones que originan resonancias, lo cual no 
había sucedido con las oscilaciones que aparecían para una única barrera 
triangular Fig. 5.16. Al igual que entonces, los mínimos comienzan a caer 
con la misma envolvente de la pendiente descendente, y son las diferentes 
alturas en las barreras de potencial las que marcan el comienzo del ascenso 
en la envolvente de estas gráficas. 


5.6. Resumen y conclusiones 


Después del estudio de una, dos y varias barreras de potencial, el siguien- 
te paso que hemos abordado es tratar de mejorar y refinar nuestro modelo. 
Para ello consideramos que sobre la lámina de grafeno no sólo influyen las 
barreras de potencial asociadas a las puertas superiores Vi, sino también 
el campo eléctrico asociado a la diferencia de potencial entre surtidor y dre- 
nador Vps, dando lugar a un potencial efectivo lineal en ciertas regiones del 
egrafeno. Enlazando con la medición experimental, la caída de potencial tiene 
lugar en la zona del grafeno que se sitúa bajo las puertas superiores, donde 
es más resistivo el grafeno que en otras zonas donde no existen puertas de 
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Figura 5.18: Dos barreras de potencial de U¡ = 0.2 (0.1) eV y Uz = 0.4 (0.3) 
eV verde (azul). En ambos casos con anchuras de 50 nm cada barrera de 
potencial, separadas 20 nm entre sí en un campo eléctrico de 4 meV/nm 
y ángulo de incidencia 10? (a). Coeficiente de transmisión para las energías 
potenciales representadas en (a) frente a la energía del portador incidente con 
ángulo de incidencia 10? (b). Observamos gran número de oscilaciones que se 
convierten en resonancias. Las diferentes alturas alcanzadas en las barreras 
de potencial marcan el inicio del ascenso en la envolvente de estas gráficas. 
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potencial por encima de él. Lo anterior nos ha llevado a definir la energía 
potencial en esas zonas de la lámina de grafeno bajo la influencia de puertas 
superiores de potencial mediante funciones lineales. Para su resolución es ne- 
cesario introducir en la definición de las funciones de onda en dichas zonas, 
tanto funciones de Hermite como funciones hipergeométricas. Lógicamente, 
en los terminales de drenador y surtidor tendremos una energía potencial 
constante, así como en las zonas del grafeno donde no estemos bajo la in- 
fluencia de una puerta superior de potencial. 


Distinguiremos entre campos eléctricos débiles cuando Fl yate << Utop gate, 
o bien campos eléctricos muy elevados cuando suceda lo contrario. Comen- 
zaremos estudiando lo que sucede para campos eléctricos débiles, que son 
los más razonables para su estudio. El principal efecto del campo eléctrico 
es favorecer el transporte en la dirección de dicho campo aplicado. El cam- 
po eléctrico provoca la pérdida de simetría en cuanto a la energía potencial, 
para un sistema compuesto por una lámina de grafeno sometida a una o dos 
barreras de potencial. Esto se traduce en la aparición de pendientes finitas 
en las paredes laterales de los mínimos o pseudogaps energéticos, con suaves 
ondulaciones, frente al caso en el cual no existía campo eléctrico que daba 
lugar a pendientes verticales con fuertes ondulaciones, cuando representamos 
el coeficiente de transmisión frente a la energía. 


Podemos actuar sobre el valor del campo eléctrico bien variando la energía 
potencial entre los contactos de fuente y drenador Us —Ubp o bien la anchura 
de la puerta superior. Aplicando uno u otro método obtendremos resulta- 
dos diferentes al representar el coeficiente de transmisión frente a la energía 
del portador incidente. En el primer caso, al aumentar el campo eléctrico 
mediante el aumento de Us — Up se provoca un aumento en la anchura del 
pseudogap energético, mientras que en el segundo caso, el aumento del campo 
eléctrico mediante la disminución de la anchura de la puerta superior provoca 
que disminuya la pendiente del pseudogap, pero éste se mantiene básicamen- 
te constante. 


Para dos barreras de potencial sometidas a campo eléctrico, el orden que 
ocupan las distintas barreras de potencial va a ser un factor importante, ya 
que influirá en el valor máximo que lleguen a alcanzar cada una de ellas. 
Sumaremos Ug a la que aparece en primer lugar y Uin+ a la que aparece en 
segundo lugar, siendo Un; < Us, y sabiendo que Uy, disminuye tanto más 
cuánto mayor sea el campo eléctrico y mayor sea la anchura de la primera 
barrera de potencial. Consecuencia de lo anterior es que la segunda barrera 
de potencial puede llegar a quedar enmascarada, es decir que sea indetectable. 
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Dos barreras de potencial sometidas a campo eléctrico con la misma 
energía potencial efectiva dan lugar únicamente a una zona donde el coe- 
ficiente de transmisión frente a la energía del portador incidente es mínimo. 
El aumento en la distancia de separación entre las dos barreras da lugar a 
un aumento en los picos que aparecen. Para ciertos parámetros particulares, 
dentro del pseudogap energético que obtenemos al representar el coeficiente 
de transmisión, aparecen claras resonancias. 


En el límite de campo eléctrico elevado, las aproximaciones realizadas de- 
jan de tener sentido, tanto por la relación de dispersión no-lineal a elevadas 
energías, como por la posible ruptura del dieléctrico utilizado. No obstan- 
te, continuaremos con el desarrollo, como una primera aproximación. Una 
barrera de potencial sometida a un campo eléctrico elevado, puede llegar 
a degradarse convirtiéndose de manera efectiva en una barrera de potencial 
triangular. En la representación del coeficiente de transmisión frente a la 
energía, para un mismo campo eléctrico y distintos valores de potencial apli- 
cados a la top gate, observamos que la envolvente de la pendiente descendente 
es semejante en todos los casos, mientras que la pendiente ascendente comien- 
za a distintos valores energéticos, relacionados con las diferentes alturas de 
las barreras de potencial. Para dos barreras con potenciales efectivos igua- 
les, observamos algo semejante aunque ahora aparecen muchas resonancias 
dentro del mínimo que no existían para una única barrera triangular. 
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Capítulo 6 


Constricciones en nanocintas de 
grafeno 


6.1. Introducción 


Este capítulo recoge el intento por nuestra parte de aplicar el formalis- 
mo de Dirac para obtener la conductancia en una constricción, inicialmente 
rectangular, realizada en una nanocinta de grafeno armchair. Para ello, a 
diferencia de capítulos anteriores, trataremos con nanocintas, las cuales pos- 
teriormente modificaremos para realizar un estrechamiento, y luego como 
extensión una constricción. En este proceso iremos viendo hasta qué punto 
es razonable la utilización de dicho formalismo. 


Una nanocinta de grafeno (GNR) es una cinta de grafeno de anchura na- 
nométrica ( típicamente inferior a 500 nm) y constante. La longitud de dicha 
cinta puede ser comparativamente mucho mayor que la anchura, por lo que 
se le puede considerar un sistema quasiunidimensional [7]. En la práctica, 
basta con que esta anchura sea inferior a la longitud de scattering de los 
electrones del grafeno. Esta pequeña anchura en la cinta da lugar a efectos 
de confinamiento cuántico, ya que los electrones básicamente se propagan a 
lo largo de la cinta, mientras que oscilan como en una onda estacionaria en la 
dirección de confinamiento (transversal). Esto difiere bastante del movimien- 
to de los electrones en una lámina de grafeno ilimitada, donde los fenómenos 
estacionarios no se manifiestan. 


La forma de los bordes de la nanocinta, principalmente zig-zag (ZGNR) 
o armchair (AGNR), [7, 14, 107-110], interviene de modo determinante en 


las propiedades de la cinta. La estructura electrónica de las nanocintas de 
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grafeno es muy sensible a la geometría de los bordes y a la anchura de la 
cinta [111]. La anchura finita de las nanocintas rompe el espectro energético 
en un conjunto infinito de bandas [109,112], en vez de una única banda como 
sucedía en el caso de una lámina de grafeno. Una cinta ZGNR es siempre 
metálica. Sin embargo, una cinta AGNR dependiendo de su anchura puede 
tener un comportamiento metálico o semiconductor [14]. Las cintas con bor- 
des zigzag tienen estados de borde (funciones de onda) localizados en el nivel 
de Fermi, mientras que esto no sucede en cintas con bordes armchair [113]. 
Los estados de borde son muy importantes cuando la anchura de la cinta es 
nanómetrica. En las cintas semiconductoras, la energía del gap E, es inver- 
samente proporcional a la anchura de dichas cintas [107,110], lo cual provoca 
que para anchuras superiores a 100 nm, E, tienda a cero, y la estructura de 
bandas tendería a ser la misma que existe en una lámina de grafeno [7]. En 
la actualidad se conocen técnicas precisas para controlar la anchura de una 
cinta [108, 112, 114, 115] y las características de sus bordes [116-118], hasta 
el punto que es posible fabricar a gran escala cintas tan estrechas como 10 
nm [119] con una muy buena movilidad de los portadores de carga, donde se 
hace necesario considerar los efectos de borde. 


Los antecedentes sobre el estudio de las constricciones en las cintas de gra- 
feno, son las constricciones para gases bidimensionales [120] que se habían 
estudiado años atrás. En 1988 tuvo lugar el descubrimiento de cuantización 
en la conductancia por van Wees [121], y Wharam [122,123] asociada a la 
constricción de un gas de electrones bidimensional por punto de contacto, que 
consiste en una constricción corta y estrecha entre dos regiones eléctricamen- 
te conductoras de mucha mayor anchura. La constricción tiene una anchura 
del orden de la longitud de onda de Fermi, y es mucho menor que la longitud 
de camino libre medio de los portadores de carga [124]. La conductancia de 
los puntos de contacto cuántico está cuantizada en unidades de 2e?/h. La 
transición entre las regiones ancha y estrecha en una constricción puede ser 
gradual (adiabática) [125] o abrupta [2, 126]. 


Las constricciones sobre nanocintas de grafeno (GNC) han sido estudiadas 
teóricamente por Muñoz-Rojas [127], Darancet [128] y Yin [129], mediante la 
teoría de enlace fuerte (tight-binding), y también usando su semejanza con 
una guía de onda de grafeno formada conectando múltiples segmentos de 
AGNR de diferentes anchuras [130,131]. Los estudios teóricos sobre constric- 
ciones dieron paso a su construcción física y al desarrollo de experimentos 
tendentes a clarificar su comportamiento [132-138]. Experimentalmente se 
han llegado a obtener constricciones con un canal muy estrecho [139] y, en 
este proceso de reducción progresiva, se ha llegado incluso hasta tamaños de 
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sencillas cadenas de átomos de carbono [140, 141]. 


Nos proponemos analizar en esta sección estos fenómenos, centrándonos 
en las nanocintas con terminación en armchair, que, debido a su versatilidad, 
parecen las más adecuadas en el diseño de dispositivos electrónicos [142]. 
El estudio teórico de las nanocintas se ha hecho tradicionalmente mediante 
métodos de tight-binding [107] o mediante la teoría de densidad funcional 
[108], mientras que trabajos analíticos basados en el formalismo de Dirac 
como el nuestro son más bien escasos. 


6.2. Estados electrónicos en una nanocinta de 
grafeno 


Para la obtención de los estados electrónicos, seguimos la metodología 
utilizada por Castro Neto [14], Brey [109] y Wurm [143]. Ahora se ha de 
tener en cuenta la existencia de confinamiento en una dirección (orientamos 
la cinta de modo que la dirección de confinamiento sea el eje x), y hemos 
de considerar la ecuación efectiva de Dirac incluyendo la contribución de los 
dos valles. El estado electrónico en la cinta correspondiente a una energía E 
debe escribirse como 


V(r,t) = al E ) 4 0 ( o ) (6.1) 


siendo 
D(r) = [La(r), Og(r), 9, (1), Pg (1)]* (6.2) 
Aplicamos la ecuación de Dirac a una nanocinta armchair de grafeno 
infinita en la dirección y, y de anchura finita W en la dirección x. Puesto 
que los bordes de la cinta que consideramos es armchair, contienen átomos 
de los dos tipos: A y B. Así, se debe exigir que se anulen a la vez en los dos 


extremos tanto la contribución a la función de onda de la subred A como de 
la subred B. Es decir 


v(0) =W(W)=0 (6.3) 


La anchura efectiva considera la primera línea de posiciones atómicas 
vacías en los bordes de la cinta, normalmente ocupada por átomos de hidrógeno, 
es decir W.;¿ = W +a [130,143], aunque la denotaremos de aquí en adelante 
como W. En la Fig. 6.1 representamos una cinta AGNR cuya dirección lon- 
gitudinal de transmisión es el eje y, con una anchura efectiva finita W. Los 
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Figura 6.1: AGNR donde representamos átomos de la subred A (círculos ro- 
jos), átomos de la subred B (círculos azules) y átomos de hidrógeno (círculos 
verdes). W es la anchura efectiva de la cinta, y suponemos longitud infinita 
a lo largo del eje y. 


bordes de la cinta están constituidos por átomos de hidrógeno representados 
por círculos verdes, siendo los círculos rojos y azules, átomos de la subred A 
y B, respectivamente. 

La función de onda electrónica queda reducida a: 


ikyy E ] 
as Es ] iO 

sin| 7 e 6.4 
lira] (A; (6.4) 
donde la componente x del momento en el cristal k, = q,, se encuentra 
cuantizada, con 


Wir) = 





nT dr 
dn = W pa K e 34 (6.5) 
siendo K la componente x del vector de onda en el punto de Dirac K. Existen 
ciertas anchuras de la cinta para las cuales existe una banda n cuyo vector 
de onda cuantizado q, es nulo. El espectro de energía para esa banda n no 
presenta gap, y por tanto el comportamiento de la cinta es metálico. En cual- 
quier otro caso será semiconductor. 
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Para una cinta dada de anchura (efectiva) W, y a una energía dada e = 
E/hvp, el número de estados propagantes viene determinado por la condición 
e >| q. |, es decir, cuando 


A O (6.6) 

TT T 

= Para una energía e dada, el número del modo más bajo posible en 
la representación de E frente a k, que contribuye a la conducción no 
es n = 0, sino nin; = Int[W(K — e)/1] + 1, siendo Int[A4] la parte 
entera de A, y el número del modo más alto de conducción €s Nsup = 
IntW(K +e)/1] . 


= El número de modos propagantes N, que aumenta con la anchura de 
la cinta W y con la energía e es: 


W W 
N= Int —(K de 0)| — Int ¡—(K E 0)] (6.7) 
T T 
= La separación entre dos modos consecutivos (Eq[6.5]) es inversamente 
proporcional a la anchura de la cinta 


T 
4n41 — Un = W' (6.8) 


Así, a medida que aumenta la anchura de una cinta, aumenta el número 
de modos que contribuyen a la conducción, a la vez que disminuye la sepa- 
ración entre dichos modos. En la figura Fig. 6.2 representamos el número de 
modos frente a la energía de los portadores para cintas AGNR de anchuras 
W = 13-a (rojo) y 26-a (azul). Ambas son semiconductoras, puesto que para 
energías próximas a O eV no existe ningún modo propagante. Para valores 
de energía no muy pequeños, el número de modos conductores de la cinta 
ancha es mayor que en la cinta estrecha. Además el intervalo de energías 
entre la aparición de dos modos sucesivos es menor en la cinta ancha que en 
la estrecha. 


En la figura Fig. 6.3 (a) se representa la energía frente a k, para una cinta 
AGNR de anchura W = 13 - a. Los modos que contribuyen a la conducción 
para valores de energía hasta 1000 meV son n = 16 (rojo), 17 (verde) y 18 
(azul). Observamos que no aparecen de manera correlativa, es decir a medida 
que aumenta n no aumenta el valor mínimo de la energía a la cual aparece 
dicho modo. En la figura (b) hacemos lo mismo para una cinta AGNR de an- 
chura W = 26-a. Al ser la cinta de mayor anchura, existen ahora más modos 
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Figura 6.2: Representación del número de modos propagantes frente a la 
energía para AGNR de anchuras W = 13-a (rojo) y 26 - a (azul). Para 
energías no demasiado bajas, en una cinta ancha aparecen más modos que 
en una cinta estrecha. 


que en el caso anterior. Los modos conductores de menor energía son n = 
32 (rojo), 33 (verde), 34 (azul), 35 (rojo discontinuo), 36 (verde discontinuo) 
y 37 (azul discontinuo). Se observa claramente que a medida que aumenta 
el número n asociado al modo (dentro de los que contribuyen), los valores 
de energía (para k, = 0) asociados a dichos modos comienzan a disminuir, 
hasta aproximarnos al origen de energías. A continuación, y siguiendo con el 
aumento de n, hacen lo contrario; es decir empiezan a aumentar, alejándose 
así del origen. 


6.3. Estrechamiento en una nanocinta 


Realizamos un estrechamiento conectando dos nanocintas AGNR de di- 
ferente anchura. En la zona 1 tenemos una AGNR ancha, de anchura efectiva 
W¡, y en la zona Il una AGNR estrecha, de anchura efectiva Wa como en la 
figura Fig. 6.4. Ambas son idealmente de longitud semiinfinita a lo largo del 
eje y, y la interface tiene lugar en y = 0. Aunque en sentido estricto hay una 
muy pequeña separación de ap/2 entre la columna con átomos de la subred 
A y la columna con átomos de la subred B, el efecto de este desplazamiento 
de red para cintas no demasiado estrechas es despreciable. En cualquier ca- 
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Figura 6.3: Representación de la energía frente a k, para los modos propa- 
gantes de una AGNR en (a) con W = 13-a para n = 16 (rojo), 17 (verde) y 
18 (azul), y en (b) con W = 26 - a para n = 32 (rojo), 33 (verde), 34 (azul), 
35 (rojo discontinuo), 36 (verde discontinuo) y 37 (azul discontinuo). No hay 
correlación entre el orden n y la energía mínima en que aparecen, sino un 
ordenamiento con efecto rebote. 


so, no suele ser tenido en cuenta en un modelo de tratamiento continuo del 
problema como es el que usamos aquí. 

Construimos las funciones de onda siguiendo la expresión general de los 
estados de la cinta Eq[6.4]. Un electrón en la zona I con energía E = shupe y 
modo «a, procedente del surtidor idealmente situado en —oo, incide sobre la 
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Figura 6.4: Unión ancho - estrecho (W; - W2) en una AGNR. Cada zona se 
supone semiinfinita a lo largo del eje y. 


interface ancho-estrecho. En dicha interface se produce en la zona 1 una onda 
reflejada con todos los momentos posibles, incluidos los correspondientes a 
ondas evanescentes, compatibles con el valor de energía de la onda incidente. 
En la interface también se produce en la zona Il una onda transmitida en 
sentido +00 con todos los momentos posibles en dicha zona. 

La conexión entre los estados en las dos zonas, ancha y estrecha, se realiza 
imponiendo la continuidad de la función de onda en el límite entre zonas eje 
y =0. Las condiciones de contorno naturales (suponemos que la cinta ancha 
termina en átomos de la red A) son: 





VA (r,0) si 0<r<W 
A a IN) SS Wa, 
OS Í 0 si W¿<zx<Wi, 000 
Vé(x,0) si 0<r<W 
B == Ie) SS Wa, 
AS l si Wa<zx<W;. (6.10) 


Dado que no hay átomos de B en la interfaz, no hay condición para 
vBP (7,0) en Wa < x < Wj, que puede tomar cualquier valor, lo cual hemos 
representado por una línea discontinua. Con respecto a las expresiones ante- 
riores Eq|[6.9-6.10], si en lugar de finalizar la cinta ancha en átomos de la red 
A, hubiese sido en átomos de la red B, intercambiaríamos en las ecuaciones 
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anteriores A por B. 

Podemos extraer las amplitudes de scattering proyectando estas ecuacio- 
nes, primero en la zona ancha de la cinta y después en la estrecha obteniendo 
relaciones suficientes para determinar los coeficientes de reflexión y transmi- 
sión para cada estado incidente. 








ar | y Wa 2 gr, Jt[5, a] 








Pel W, 01. 
(6.11) 
Wa k1y [0 ] = Wa eel , 
—21 W, 0:10] mis pra e ni Sl n]S!j, m)] — Sn )t[m, 0] 
(6.12) 


siendo 





Wa bsin lar 
S/a, b] = | sin ES sin ÉS de = SN ; — | (6.13) 
0 T LN 
Wi 


En la práctica, para cada cinta en concreto y para cada valor de energía, 
el sistema anterior Eq[6.11-6.13], que constituye un sistema lineal de ecua- 
ciones algrebraicas, tiene que resolverse numéricamente. 


Las condiciones de contorno impuestas sobre las funciones de onda y su 
continuidad en el límite entre zonas, hace que la densidad de corriente sea 
también continua. Puesto que la densidad de corriente y la corriente son 
independientes de la coordenada y, podemos elegir las expresiones para la 
misma en cualquiera de las dos regiones, puesto que son idénticas. Elegimos 
la expresión correspondiente a la corriente transmitida, como se acostumbra. 
La intensidad de corriente transportada por un estado incidente q: procedente 
del surtidor y la transmitida al drenador, transportada por todos los estados 
propagantes, resultan ahora, para este nuevo sistema contable de estados que 
caracteriza la nanocinta de anchura finita, 


ev 
cl l= — sky la al, Tesoa [01] A y [tl $ 0 al kay | Jl, (6.14) 


jprop 
y el coeficiente de transmisión para el modo incidente a: 


A (6.15 
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Figura 6.5: Conductancia frente a la energía para un estrechamiento de una 
AGNR, siendo W, = 26 -a y Wa = 13 - a (rojo). Conductancia ideal de la 
cinta estrecha W = 13-a (lila). La conductancia a través del estrechamiento 
viene determinada básicamente a partir del número de modos propagantes 
de la cinta estrecha. 


donde la suma recorre solamente los estados propagantes (y no los evanescen- 
tes), que son los que contribuyen a la corriente. Para obtener el coeficiente 
de transmisión total hay que hacer la suma sobre todos los modos entrantes 
O. que se propagan por la cinta, siendo la conductancia proporcional a ella. 


e > Tlo] (6.16) 


En la figura Fig. 6.5 representamos en rojo la conductancia frente a la 
energía de los portadores para una cinta AGNR que presenta un estrecha- 
miento, pudiendo considerarse como dos nanocintas semiinfinitas con an- 
churas W, = 26 -a y Wa = 13- a, mientras que en lila representamos la 
conductancia ¿deal asociada al número de modos propagantes de la cinta se- 
miconductora estrecha W = 13-a que contribuyen a la conducción, teniendo 
en cuenta que ¿dealmente cada modo va a contribuir con una conductancia de 
2e” /h puesto que suponemos longitud infinita y por tanto la transición entre 
escalones es abrupta al no existir scattering [144]. Sin embargo, si hubiese 
sido una cinta armchair metálica o zigzag los plateaus estarían separados 
4e?/h ya que entonces se preservarían las simetrías de spin y de valle [88]. 
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Figura 6.6: Constricción en una AGNR : W; - Wa - W3(D), siendo D la dis- 
tancia entre las zonas extremas, semiinfinitas a lo largo del eje y. La relación 


entre D y el número de hexágonos h es : D = oa . y 


Aunque el número de modos de la cinta ancha y de la estrecha es diferente, 
la conducción a través del estrechamiento va a venir definida principalmente 
mediante el número de modos propagantes en la cinta estrecha. 


6.4. Nanoconstricción rectangular en una na- 
nocinta 


Una nanoconstricción de grafeno (GNC) es un estrechamiento localizado, 
de una longitud nanométrica, en una nanocinta de grafeno en el camino de 
la corriente de portadores. Cuando las dimensiones del estrechamiento, largo 
y ancho, son semejantes, a la nanoconstricción se le llama punto cuántico de 
grafeno (GQPC) [145]. A continuación consideramos una nanoconstriccción 
con forma rectangular en una cinta AGNR, que denotamos Wi - Wa - W3(D) 
y medimos en hexágonos. Básicamente, este sistema constituye una amplia- 
ción del caso anterior, donde ahora tenemos tres zonas (en vez de dos) como 
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se muestra en la figura Fig. 6.6, claramente diferenciadas por sus anchuras en 
la AGNR. Las zonas l, 11 y III de anchuras respectivas W,, Wa y W3, deben 
verificar las condiciones Wa < W, y Wa < W3 para que efectivamente sea 
una constricción. La distancia entre la zona 1 y la zona II es D, por lo que 
la constricción viene definida a partir de su anchura Wa y longitud D, siendo 
D = [(3h-1)/2] -(a/v3), siendo h el número de hexágonos. Los límites entre 
las tres zonas están situados por tanto en y = y; y en y = ya = y1 + D. Se 
pueden considerar nanoconstricciones tanto asimétricas con W, 4 W3, como 
simétricas con W, = W3. Para este trabajo se ha elaborado el cálculo general 
aplicable a ambos tipos de nanoconstricciones, pero las simulaciones se han 
realizado sólo para el caso simétrico. 


Para aplicar las condiciones de contorno, consideramos el acabado de la 
cinta. Elegimos que el borde de la nanoconstricción en y = y termine en 
átomos de la subred A, y debido a la geometría de la red hexagonal del gra- 
feno, tendremos que el borde en y = yz terminará en átomos de la subred 
B. Si hubieramos hecho la suposición contraria, intercambiaríamos A por B. 
Como consecuencia de la simetría del hamiltoniano, el coeficiente de transmi- 
sión no cambiaría para dicho caso. Despreciamos, de nuevo la muy pequeña 
separación de ap/2 entre la columna con átomos de la subred A y la columna 
con átomos de la subred B. 


Las condiciones de contorno que hemos de imponer son: 








v7 (1, y1) = Í nte 1) E AN (6.17) 
VE (xy) = l (ay) E a a (6.18) 
VA (1, Y2) = Í unto y2) a A (6.19) 
v?, (2, y») = l Vil, 92) E A (6.20) 


Aplicamos estas condiciones de contorno a las funciones de onda de las 
tres zonas, y proyectamos el estado electrónico en las zonas anchas y estrecha 
de la cinta, obteniendo el siguiente conjunto de ecuaciones: 
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rlp,a] = 
Qu» e le a e 
(Q1+[p] Opa + 2 me (cel Le Q1+ lp] Caló, ) : 
(6.21) 
2 po 51ALi pltl;, a] = 
[pe PRkIC, [p, a] + Q2+ [pj PRPIC_ lp, a], (6.22) 


W. 
Em Sls,pr[i, o 7 Sla,p] +Cy[p,0]+0C_[p,0], (6.23) 


3: 


00 W. ta : a 
al=Y_ y y, lp, 5) (Publ, [j, a] + e PrlilC_[,0]), (6.24) 


siendo Sla,b] como en la expresión Eq[6.13] y Rla, b] 











T (2 


Las ecuaciones Eq[6.21-6.24] resuelven el problema de determinar, para 
un estado entrante dado, todos los coeficientes por medio de un sistema lineal 
de infinitas ecuaciones algebraicas. En la práctica, para unos valores deter- 
minados de los parámetros de la cinta, se resuelve el sistema lineal anterior 
truncándolo en un valor conveniente que determina la precisión del cálculo. 
Se calcula el coeficiente de transmisión de cada modo conductor por medio 
de la expresión Eq[6.15], y luego la conductancia de la cinta con Eq[6.16]. 


Enfocamos nuestro análisis en la dependencia del transporte con factores 
geométricos tales como la relación entre las anchuras de las zonas extremas 
y de la zona intermedia, así como la longitud de ésta. Vamos a representar 
la conductancia de la cinta Eq[6.16] en función de la energía de los porta- 
dores para cuatro diferentes combinaciones: zona intermedia metálica (M) 
o semiconductora (S) con extremos metálicos (M) o semiconductores (S). 
Comparamos en la misma gráfica nuestros resultados con los obtenidos por 
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Figura 6.7: (a) Conductancia para las constricciones 6 - 3 - 6 (9) mediante 
las ecuaciones de Dirac en azul y tight-binding en rojo. (b) Idem para 7 - 3 
- 7 (9). Obsevamos un buen acuerdo entre los métodos de las ecuaciones de 


Dirac y tight-binding. 


el grupo Quantum Nanosystem Group de la Universidad Complutense en la 
simulación de estas mismas cintas mediante el método tight-binding. 


En la figura Fig. 6.7 (a) representamos el caso 6 - 3 - 6 (9) correspon- 
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Figura 6.8: (c) Conductancia para las constricciones 6 - 4 - 6 (9) mediante 
las ecuaciones de Dirac en azul y tight-binding en rojo. (d) Idem para 7 - 4 
- 7 (9). Obsevamos un cierto acuerdo entre el método de las ecuaciones de 


Dirac y el tight-binding. 


diente a una cinta de tipo M-M-M, y en (b) el caso 7 - 3 - 7 (9) de tipo 
S-M-5S. En ambos casos la zona intermedia es metálica, y el acuerdo entre el 
método utilizando las ecuaciones de Dirac desarrolladas y el método tight- 
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binding es bueno. En la figura (b) la conductancia comienza a crecer cuando 
se alcanza la energía a la cual empieza a contribuir el primer modo de la 
cinta más ancha. En la figura Fig. 6.8 (c) representamos el caso 6 - 4 - 6 
(9) correspondiente a una cinta de tipo M-S-M, y en (d) el caso 7 - 4 -7 (9) 
de tipo S-S-5. Es decir, en ambos casos la zona intermedia es semiconduc- 
tora, y el acuerdo entre ambos métodos se aleja sobre todo en los mínimos, 
aunque en el máximo local coinciden. En todas estas figuras Fig. 6.7 - 6.8, 
la conductancia está condicionada principalmente por las características de 
la cinta más estrecha. En las cintas AGNR de anchura 3 o 4 hexágonos (o 
constantes de red) que se están simulando, existe como máximo un modo 
activo de conducción, y por ello en las representaciones, la conductancia no 
supera el valor unidad en unidades de 2e?/h. 

Inicialmente habíamos pensado utilizar el método de la ecuación de Dirac 
para resolver una nanoconstricción porque el método tight-binding a partir 
de ciertas dimensiones resultaba impracticable. Sin embargo, nuestra sorpre- 
sa surgió cuando observamos que nuestro método presentaba problemas de 
convergencia, incluso para anchuras de la cinta relativamente pequeñas. Así, 
el método tight-binding resultó más ventajoso. 


Buscábamos investigar las propiedades de transporte de los electrones 
y los efectos de cuantización asociados al tamaño en una GNC. Para ello 
nuestro grupo de investigación Nanotechnology Group de la Universidad de 
Salamanca ha utilizado un método de fabricación de grafeno que permite un 
buen procesado litográfico con elevada calidad cristalina, al mismo tiempo 
que mantiene una elevada movilidad de los portadores [146,147]. La rugo- 
sidad del borde en la muestra experimental es pequeña, es decir del orden 
de 1ó 2 nm, lo cual es conveniente ya que un aumento en la rugosidad con- 
lleva un aumento en el scattering. En la figura Fig. 6.9 en (a) se muestra 
una imagen SEM de nuestro dispositivo, que consiste en una barra de gra- 
feno de 1x2 ¡um?, y una constricción 85x85 nm?, y en (c) se representan 
las dimensiones geométicas. Las muestras utilizadas tienen un pico de Di- 
rac próximo a cero voltios debido a la baja interacción del grafeno con el 
substrato. GNC similares basadas en heteroestructuras de grafeno y h-BN 
eran estudiadas por Terrés [135], que aunque también consiguieron eleva- 
da movilidad para los portadores, la anchura de la constrición superaba los 
200 nm debido al proceso de grabado utilizado en la heteroestructura. Otro 
tipo de GNC para anchuras menores de 200 nm era realizado con grafeno 
encapsulado por Kumar [138], sin embargo la rugosidad en el borde no era 
despreciable. Anteriormente Tombrós [134] había realizado pequeñas GNC 
en grafeno suspendido con elevada movilidad, sin embargo la técnica de fa- 
bricación utilizada no estaba bien definida. Otro tipo considerado era GNC 
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Figura 6.9: (a) Imagen SEM de una cinta de 1x2 ym? y una GNC de 85x85 


nm”. (c) Imagen geométrica de la constricción. 


sobre obleas de Si/SiO>, que presenta baja rugosidad de borde y un buen 
control del tamaño, sin embargo la movilidad de los electrones era muy ba- 
ja debido al contacto directo entre grafeno y óxido [132,137], lo cual evita 
una buena observación de los efectos asociados a la cuantización por tamaño. 


En la figura Fig. 6.10 se representa la conductancia frente a la magnitud 
adimensional W - kp, donde W es la anchura de la constricción y kp = yrn 
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Figura 6.10: Conductancia frente a la magnitud adimensional W -kp. La línea 
negra corresponde a la medida a baja temperatura, mostrando la tendencia 
lineal esperada, cuyo ajuste se representa mediante una línea discontinua. La 
línea azul muestra el resultado de la simulación teórica. 


es el número de onda de Fermi. La densidad de portadores n estaba relacio- 
nada con la tensión back-gate V¿y mediante n = (Coz/e)(Vog — Von) donde 
Cor es la capacidad del SiO2 por unidad de área, e es la carga elemental, y 
Vonp es el voltaje en el punto de Dirac. La línea azul muestra el resultado 
de la simulación teórica realizada mediante la aproximación tight-binding a 
los vecinos más próximos por colaboradores pertenecientes al grupo Quan- 
tum Nanosystem Group de la Universidad Complutense de de la GNC con 
bordes rectos, donde se ha supuesto que la cinta de partida es armchair. La 
línea negra corresponde a la conductancia medida a baja temperatura (5 K), 
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mostrando la tendencia lineal esperada, asociada al régimen balístico, cuyo 
ajuste se representa mediante una línea discontinua. Esto es cierto salvo para 
valores bajos de W - kp donde aparece una pequeña contribución adicional 
debida a un dopado residual en la muestra de grafeno, lo cual evita poder 
realizar experimentalmente una densidad de portadores nula. Los resultados 
de la figura Fig. 6.10 son semejantes a los obtenidos por Yannouleas [131] que 
encontró picos en la conductancia de sus nanoconstricciones y los atribuyó a 
la cuantización longitudinal debida a la longitud finita de la GNC. Las pen- 
dientes de las líneas negra y azul son muy semejantes, y la pequeña diferencia 
entre ambas procede de imperfecciones como bordes más redondeados en las 
esquinas de la GNC real. En general podemos decir que los datos obtenidos 
experimentalmente y las simulaciones numéricas muestran un buen acuerdo. 


6.5. Resumen y conclusiones 


El propósito de este capítulo fue inicialmente utilizar el formalismo de Di- 
rac, para aplicarlo en una nanoconstricción realizada sobre una cinta AGNR, 
que se iba a realizar en el laboratorio. Para ello se comenzó aplicando el 
formalismo a nanocintas AGNR, en las cuales debido a la anchura finita, el 
vector de onda transversal resultaba cuantizado, dando lugar a un conjunto 
finito de modos de propagación; lo cual constituye la diferencia fundamental 
con el comportamiento de la lámina de grafeno ilimitada. El número de mo- 
dos propagantes NV en una AGNR aumentaba con la anchura de la cinta W 
y con la energía e. Posteriormente realizamos un estrechamiento de una cinta 
AGNR, obteniendo que la conductancia viene dada a partir del número de 
modos que se propagan en la parte estrecha de la cinta. Para comprobarlo, 
comparamos la conductancia del estrechamiento de una AGNR con la con- 
ductancia ideal asociada a la cinta estrecha. 


Siguiendo un proceso semejante al utilizado para el estrechamiento, es- 
tudiamos una nanoconstricción rectangular para una cinta AGNR, con el 
propósito de poder utilizarlo cuando en el laboratorio obtuviesemos la nano- 
constricción de grafeno. Una vez obtenido el desarrollo matemático, lo apli- 
camos a una cinta estrecha para comparar los resultados con los obtenidos 
mediante el método tight-binding. Encontramos que cuando la constricción 
era metálica los resultados mostraban un acuerdo relativamente bueno, que se 
separaba claramente cuando la zona intermedia era semiconductora. Además 
al aumentar la anchura de la cinta, nuestro método presentaba problemas de 
convergencia. Por todo ello, junto con las ventajas que presentaba el método 
tight-binding de ser escalable y poder definir la forma de los bordes, hizo 
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aconsejable utilizar este último método para realizar la comparación con los 
resultados experimentales. En nuestro grupo se realizó una GNC mediante 
un proceso litográfico que estaba bien definido, con una rugosidad de bor- 
des muy pequeña y que además presentaba una elevada movilidad para los 
portadores. Posteriormente se obtuvo una dependencia lineal entre la con- 
ductancia y kp, como cabía esperar para un régimen balístico, tanto desde el 
punto de vista experimental como a partir de la aproximación tight-binding 
utilizada. 
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Apéndice A 


Densidad de portadores 


A.l. Expresión general para la densidad de 
portadores 


Para calcular la densidad de portadores en el grafeno, necesitamos conocer 
la densidad de estados g(H) - el número de estados en un intervalo unitario 
de energía para una superficie unitaria de grafeno-, así como la probabilidad 
de que éstos estén ocupados f(L). Si en una superficie unitaria de grafeno 
tenemos dZ estados cuánticos en el intervalo de energías entre E y E+dbE, 
entonces 
_dZ 
dE 


La probabilidad de llenado del estado con energía E es igual a f(L), por 
lo que el número de electrones es: 


g(E) (A.1) 


dn= f(EJdZ = (E)(EJdE (A.2) 


El número total de electrones cuyas energías están comprendidas entre 
Es y E, es: 


n= |" HEEE (A.3) 
E 
En dos dimensiones [7], el número total de estados en el intervalo de 
energías entre E y E+dE es el área diferencial en el espacio recíproco K 
dividido entre el área que ocupa un estado k de este espacio 


2rTkdk 
mM . 


g(EME = (2m) 





(A.4) 
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Puesto que E = hupk podemos escribir la expresión Eq|A.4] en función de 
la energía 
1. EdE 2 EdE 
EME =m:- . SH A5 
E) 2 (hup)?? Tr  (hup)? ES) 
puesto que m =4=2 x 2, un 2 asociado a la degeneración de spín, y el 
otro 2 asociado a la degeneración de valle, por lo que la densidad de estados 
es: 





inst (A.6) 


Tr (hup)? 


La densidad de estados es un número positivo, por lo que para huecos 
(con energía negativa) es necesario tomar el módulo en la expresión anterior. 

La probabilidad de ocupación para electrones está dada por la función de 
distribución de Fermi-Dirac 


1 


A 14 e(E-Ep)/kaT (A.7) 





La densidad de electrones n en equilibrio, dentro de todo el rango posible 
de energías, es: 





2 ds E 
0 e / rie (A.8) 
A continuación se hace el cambio de variable y = E/kgT y np = Ep/kgT 
2 kgT d E n 
=-=| == ————d A.9 
Ñ 2) 0 14 ermano 0 


donde esta integral representa la integral de Fermi-Dirac de orden 1, ya 
que en general la integral de Fermi-Dirac de orden s se expresa como: 


O 00 n 
1 / Ea dn (A.10) 


Por lo que la expresión exacta para la densidad de electrones [148] es: 


2 (kgT Y' Ep 
=-=| ==] Fil —Ñ A.11 
ñ T (E ) : (E ) 
Se pueden analizar los casos extremos: caso no degenerado (np << 0) y 
el caso degenerado (np >> 0). 
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Para el caso no degenerado Ep << kgT' , tenemos la aproximación 
Fi(np) > e”, si además la energía de Fermi es nula obtendríamos una 
dependencia cuadrática con la temperatura 


2 EEN" 
=p A.12 
8 (E) ) 


T 


Para el caso degenerado, cuando Ep >> kgT' entonces Fi(np) > np”, 
por lo que la densidad de portadores [149] es: 


2 (kg T Y” 1(ErY' 1(EpY 
n 2 == E] = - (A.13) 
TA hup 2 A kgT TA hvup 
Debido a la relación de dispersión simétrica para electrones y huecos, 
obtendríamos para huecos una expresión semejante a Eq[A.11]. Además, te- 


niendo en cuenta que la función de distribución de huecos es 1 — f(B), la 
expresión exacta para la densidad de huecos es: 


DARE? HE 
E e A.14 
p T (E) e ( Es) ( ) 
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Hemos calculado analíticamente el coeficiente de transmisión para una 
lámina de grafeno de anchura indefinida sometida a una o dos barreras de 
potencial unidimensional. Para construir los estados electrónicos, hemos te- 
nido en cuenta las simetrías del hamiltoniano de Dirac que describe a los 
electrones en la lámina. Se concluye que tanto la energía de los portadores 
como la componente perpendicular de su vector de onda, deben mantenerse 
constantes. La relación entre el ángulo de incidencia y el de transmisión en 
una barrera de potencial depende de la energía del electrón incidente y de la 
altura de la barrera de potencial Uy. Bajo ciertas condiciones (E > 0.5 Uy) 
aparece un ángulo límite, de tal forma que para ángulos superiores, la trans- 
misión es mediante estados evanescentes en lugar de estados propagantes. 

El coeficiente de transmisión resulta unitario para incidencia perpendicu- 
lar, y también en otros ángulos que verifiquen ciertas condiciones, ello provoca 
que no pueda existir confinamiento de portadores de carga mediante una ba- 
rrera de potencial eléctrico. El coeficiente de transmisión es mínimo, próximo 
a cero, para ciertos valores de los parámetros, de manera que podemos llegar 
a obtener un pseudogap de energía. 

El mismo modelo ha podido aplicarse al caso de dos barreras de poten- 
cial independientes. El coeficiente de transmisión resultante es independiente 
del orden en el cual coloquemos las barreras de potencial. Asociada a cada 
barrera de potencial puede aparecer un pseudogap, que dependiendo de los 
parámetros considerados, llegue a colapsar en un único pseudogap. Al con- 
siderar dos barreras de potencial simétricas aparecen unos picos (tipo delta) 
de resonancia equiespaciados dentro del pseudogap. En la medida en la cual 
las dos barreras dejen de ser simétricas, esos picos de resonancia se van de- 
teriorando, disminuyendo de amplitud, desplazándose o ensanchándose. 


A continuación, hemos construido un modelo para el transistor de efecto 
de campo de grafeno (GFET) que relaciona los potenciales exteriores apli- 
cados en las puertas del dispositivo con los parámetros utilizados en nuestro 
modelo. Las variables experimentales que gobiernan el GFET están asocia- 
das a los potenciales aplicados en los terminales de fuente, drenador, puerta 
superior y puerta inferior. Éstas están asociadas al perfil de la barrera de 
potencial considerada y a la energía de Fermi de los portadores. Obtenemos 
una expresión para la conductancia de los GFET bajo ciertas condiciones ex- 
perimentales, en particular para temperaturas próximas a cero kelvin, y baja 
polarización. En principio, consideramos incidencia no-perpendicular, para 
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evitar la paradoja de Klein. Además dependiendo del ángulo de incidencia 
de los portadores sobre la barrera de potencial, la conductancia obtenida es 
diferente. Inicialmente consideramos que los portadores pueden proceder de 
todas las direcciones posibles y además llegar a impactar en la barrera de 
potencial con la misma probabilidad. Posteriormente, tratamos el caso de 
direcciones privilegiadas, es decir que tienen más probabilidad de proceder 
de esas direcciones que de otras, lo cual simulamos mediante funciones de 
ponderación tipo campana de Gauss o tipo delta. Físicamente, esta disper- 
sión que puede suceder al considerar direcciones privilegiadas, la asociamos 
a muestras con pocas o bastantes inhomogeneidades. 


Aplicamos múltiples barreras de potencial a una lámina de grafeno, y en 
vez de obtener una expresión analítica cerrada (como había sucedido para dos 
barreras de potencial), utilizamos el Método de la Matriz de Transferencia 
(MMT). En el caso de múltiples barreras de potencial con la misma altura, 
realizamos representaciones para el coeficiente de transmisión frente a varios 
parámetros, como el ángulo de incidencia que da lugar a la aparición de reso- 
nancias, la energía con la aparición de minibandas y un pseudogap semejante 
al obtenido para una única barrera de potencial, pero con la aparición de re- 
sonancias muy selectivas en energía dentro de ese pseudogap. En el caso en 
que varias barreras de potencial sean de diferente altura, observamos que el 
resultado depende del ordenamiento de dichas barreras, de tal forma que si 
existe simetría especular el resultado es el mismo. Asociado a cada barrera 
de potencial aparece un mínimo en energía que dependiendo de los paráme- 
tros usados se puede convertir en un pseudogap en energía, de tal forma que 
varios consecutivos podrían llegar a colapsar en un único pseudogap. Para di- 
señar los extremos inferior y superior de un pseudogap necesitamos al menos 
dos barreras de potencial de diferentes alturas. La posibilidad de crear zonas 
de transmisión y de pseudogaps mediante voltajes aplicados externamente 
ofrece la posibilidad de construir un filtro dependiente de la energía. 

Por otra parte, la barrera de potencial escalonada permite simular cual- 
quier potencial que varíe suavemente. Nos ha parecido interesante aplicarlo 
a una función gaussiana. En la representación de su coeficiente de transmi- 
sión aparecen varios mínimos en energía, siendo el menor áquel asociado a la 
barrera de potencial más elevada, a partir de dicha energía es unitario. 


Con el fin de mejorar y refinar nuestro modelo, hemos considerado que 
sobre la lámina de grafeno además de influir las barreras de potencial asocia- 
das a las puertas superiores e inferior Vra — Vga , también lo hace el campo 
eléctrico asociado a la diferencia de potencial entre surtidor y drenador Vhs, 
dando lugar a un potencial efectivo lineal en ciertas regiones de la lámina 
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de grafeno. Lo anterior, cambia el potencial y por tanto el Hamiltoniano, lo 
que nos lleva a construir de nuevo los estados electrónicos, que ahora son 
diferentes. Ya no resultan de una superposición adecuada de ondas planas, 
sino que es necesario introducir funciones de Hermite e hipergeométricas. 


La consideración del campo eléctrico frente al caso en el cual básicamente 
lo habíamos considerado despreciable, provoca la pérdida de simetría para 
la energía potencial en un sistema compuesto por una lámina de grafeno 
sometida a una o dos barreras de potencial. Así, en la representación del 
coeficiente de transmisión frente a la energía aparecen pendientes finitas en 
las paredes laterales de los mínimos o pseudogaps energéticos, con suaves 
ondulaciones, frente al caso en el cual no existía campo eléctrico que daba 
lugar a pendientes verticales con fuertes ondulaciones. Para variar el valor 
del campo eléctrico podemos hacerlo bien variando la energía potencial entre 
los contactos Ug — Up, o bien la anchura de la puerta superior, dando lugar 
a resultados diferentes. En el primer caso el aumento de Ug— Up provoca un 
aumento en la anchura del pseudogap energético, mientras que en el segundo 
caso, la disminución de la anchura de la puerta superior provoca una dis- 
minución en la pendiente del pseudogap, pero éste se mantiene básicamente 
constante. 


Al someter dos barreras de potencial a campo eléctrico aparecen nuevos 
efectos que no tenían lugar cuando no existía campo eléctrico. El orden que 
ocupan las distintas barrreras de potencial influye en el valor máximo que 
lleguen a alcanzar cada una de ellas. En consecuencia, la segunda barrera de 
potencial puede llegar a quedar enmascarada, es decir que sea indetectable. 
También estudiamos dos barreras de potencial sometidas a campo eléctrico 
que se reducen de modo efectivo a barreras idénticas y dan lugar únicamente 
a un mínimo al representar el coeficiente de transmisión frente a la energía 
del portador incidente. El aumento en la distancia de separación entre las dos 
barreras proporciona un aumento en el número de picos que aparecen, y para 
ciertos parámetros particulares, aparecen claras resonancias en el pseudogap 
energético. 


En el límite de campo eléctrico elevado, la barrera de potencial puede 
llegar a degradarse, convirtiéndose de manera efectiva en una barrera de po- 
tencial triangular, que estudiamos - aunque las condiciones utilizadas no sean 
suficientemente adecuadas - como una primera aproximación. En la represen- 
tación del coeficiente de transmisión frente a la energía, para un mismo campo 
eléctrico y distintos valores de potencial aplicados a la top gate, observamos 
que la envolvente de descenso en el mínimo energético es aproximadamen- 
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te igual, mientras que la pendiente ascendente comienza a distintos valores 
energéticos, relacionados con las diferentes potenciales aplicados a la top 
gate. Para dos barreras con potencial efectivo igual, observamos algo seme- 
jante aunque ahora aparecen muchas resonancias dentro del mínimo que no 
existían para una única barrera triangular. 


Finalmente, hemos considerado que nuestra lámina de grafeno se estre- 
cha a lo largo de una dirección, hasta obtener una nanocinta de grafeno, que 
al tener una anchura finita se caracteriza por un vector de onda transversal 
cuantizado, dando lugar a un conjunto finito de modos de propagación, a di- 
ferencia de lo que sucedía en una lámina de grafeno indefinida. La nanocinta 
idealmente infinita, sufre a continuación un estrechamiento, de tal forma que 
se comporta como la unión de una cinta ancha y otra estrecha. La conductan- 
cia básicamente viene dada a partir del número de modos que se propagan 
en la parte estrecha de la cinta, para lo cual comparamos la conductancia del 
estrechamiento con la de la cinta más estrecha. 


El siguiente paso en esta reducción progresiva de las dimensiones, ha 
sido considerar una nanoconstricción rectangular en una nanocinta, con el 
propósito de poder utilizar los resultados para realizar una comparación con 
una nanoconstricción que ibamos a realizar en el laboratorio. Comparamos 
para una nanoconstricción estrecha nuestros resultados con los obtenidos me- 
diante tight-binding, obteniendo que cuando la zona de la constricción era 
metálica los resultados mostraban acuerdo, que se separaba cuando la zona 
intermedia era semiconductora. Por lo cual, junto con las ventajas que presen- 
taba el método tight-binding de ser escalable y poder definir la forma de los 
bordes, hizo aconsejable utilizar este último método para realizar la compa- 
ración con el caso experimental. Nuestro grupo realizó una nanoconstricción 
de grafeno mediante un proceso litográfico que estaba bien definido, con una 
rugosidad de bordes muy pequeña y elevada movilidad para los portadores. 
Obtuvimos una dependencia lineal entre la conductancia y kp, como cabía 
esperar para un régimen balístico, tanto desde el punto de vista experimental 
como a partir de la aproximación tight-binding utilizada. 
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